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I. SUMMARY OF PRINCIPAL SCIENTIFIC RESULTS 
FROM THE LUNAR HEAT-FLOW EXPERIMENT 


The Apollo Lunar Heat-Flow Experiment was conceived to make direct 
observations of the rate of heat loss through the surface of the moon at 
the Apollo landing sites. Two of these experiments were emplaced success- 
fully during the Apollo program: one at Hadley Rille (Apollo 15) and the 

other at Taurus Littrow (Apollo 17) . A third experiment was flown on 
Apollo 13, but when the mission was aborted en route it never reached the 
moon. A fourth instrument was installed at Apollo 16, but a broken cable 
between the central transmitting station and the heat-flow experiment ren- 
dered it useless. The principal components of the experiment were probes, 
each with twelve thermometers of exceptional accuracy and stability, chat 
recorded temperature variations at the surface and in the regolith down to 
2.5 m. The Apollo 15 experiment recorded temperatures for a period of four 
and a half years and the Apollo 17 probes, which were still returning accurate 
results when the instrument was turned off, recorded four years and ten 
months of lunar surface and subsurface temperatures. These data provided a 
unique end valuable history of the interaction of solar energy with lunar 
surface and the effects of heat flowing from the deep interior out through 
the surface of the moon. The interpretation of these data resulted in a 
clearer definition of the thermal and mechanical properties of the upper 
two meters of lunar regolith, direcc measurements of the gradient in mean 
temperature due to heat flow from the interior and a determination of the 
heat flow at the Apollo 15 and Apollo 17 sices (see Table 1 for values 
of significant new observations that resulted from the experiment.) 

S urf ace Temperatures and the Nacure_of the Shallow Reg o lit h 

Surface thermometers of the heat-flow experiment provided a complete 
history of the surface variation of temperature during the observation period. 
Two important discoveries stemmed from the interpretation of these results: 

1) A 30 to 35® difference in mean cemperature occurs across the upoer 
few centimeters of the regolith because, in this relatively fluffy layer, 
radiative heat transfer predominates. At lunar noon, nearly 70% of the 
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heat transfer In this layer occurs by particle to particle radiation. 

2) The thermal response of the regolith to variations In solar radia- 
tion, due to eclipses, rotation of the moon on its axis, and the eccentricity 
of its orbit about the sun, is weLLdescribed by a simple two layer model of 
the upper 2.5 m of the regolith, which is comprised of a thin (2-3 cm) surface 
layer of low density (1.1 - 1.2 gm/cc) very fine material of extremely low 
conductivity (0.9-1. 6 mW/m 9 K) . This layer overlies a higher density (1.75 - 
2.10 gm/cc) tightly compacted layer of fines which have a quite uniform 
conductivity of 10 to 15 mW/m 9 X. The thermal conductivity of the deeper 
layer was determined by careful analysis of the downward propagation of long- 
term surface variations of temperature, such as the annual variations and 
disturbances caused by the visits of the astronauts to the sites. The 
value of conductivity inferred for this deeper layer is considerably higher 
than can be achieved by compaction of lunar sells in the laberatory, sugges- 
ting that on the moon soils below a few centimeters have been compacted into 
a mechanical configuration that has evolved over millions of years which 
cannot be duplicated in the laboratory. This characterization of the lunar 
soil layer will be important when the time comes to establish bases or 
scientific stations on the lunar surface. 

T he Temperature Gradient and Hea t Flow 

The increase of mean temperature with depth, due to heat flow from the 
interior, varied between the four sites, due to local disturbances and regional 
variations in surface heat flow. The average gradient observed was 1.6 9 K/m. 
This value is in agreement with the vertical temperature gradients deduced 
from radiotelescope observations of long wavelength radiation from the whole 
disk of the moon. The observed increase in brightness temperature with 
wavelength when combined with the electromagnetic properties measured on 
Apollo soil samples. Indicate a gradient of 1.5 to 2.0 9 K/m. The significance 
cf this agreement is that it indicates that the mean of the four Apollo 
measurements is probably representative of the nearside of the moon. 

The observed mean gradient of 1.6*K/m and the thermal conductivity of 

2 

the regolith below 3 cm combine to yield a heat flow of 17 mW/m with an 
estimated uncertainty of 20 %. How representative this value is of the true 



average heat loss from the moon is difficult to say with so few data that 
vary so widely. There are some lines of evidence that suggest it may be 
quite representative. The large variations of heat flow which are observed 
over the surface of the earth are mainly associated with the boundaries cf 
converging and diverging lithospheric plates. The lithosphere of the moon, 
on the other hand, is rigid and static so that variations associated with 
tectonics are unlikely. Regional variations in heat flow on the mocn are 
more likely associated with variations in abundance of long-lived radio- 
active isotopes of potassium, thorium and uranium. These isotopes were 
mobilized by the igneous activity that formed the lunar crust and flooding 
cf the mare basins. Orbital gamma-ray observations, global photomapping 
and Apllo samples allow seme estimate of the variability of heat flow. 

Our analysis suggests that this variability is small. It is described 

in Section VI-4. Lastly, a critical analysis of crust and mantle tem- 

2 

peratures that would be implied by a heat flow between 12 and 18 mW/m 
gives values that are in good agreement with seismic and magnetic data that 
depend on mantle temperatures. The heat-flow result implies temperatures 
within a few hundred degrees of melting at depths of about 300 km. 

For a planet as small as the moon, much of the heat associated with 
accretion will have escaped during its 4.5 b.y. history; consequently, most 
of the heat escaping through the surface at present is produced bj radioactive 
isotopes in the interior. Thus, the moon is close to steady state. If we 
assume that this is the case for the moon, then the observed heat-flow 

_9 

value would require a uranium content cf approximately 30-45 x 10 g/g which 
is consistent with the results of more detailed thermal models. This amount 
is similar to that estimated for the mantle of the earth. 
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II. EARLY ORGANIZATION OF THE 
APOLLO LUNAR HEAT-FLOW EXPERIMENT 


Even chough. In concept, Che heac-flow experiment was simple, it in fa.L 
applied physical principles chac have been well understood for more than 
a century, the development of any Apollo experiment was a complex under- 
taking, usually involving hundreds of persons, a dozen or more corporations 
and several branches of NASA. Because of this complexity, this history will 
not attempt to cover every event or person involved. There were many acti- 
vities about which I am totally unaware and that only a long research effort 
could uncover. The purpose of this section is rather to record the interes- 
ting and sometimes instructive story of the twisting path the development of 
our experiment followed. It will, in some cases, describe the results of 
efforts which were eventually discarded for one reason or another, and conse- 
quently will not appear in any scientific report. 

The importance of knowing the total loss of heat energy from the surface 
of a planet has long been recognized as an important constraint in thermal 
history and total average composition of chac planet. Important papers by 
Urey (1957), Levin (1962), MacDonald (1959) and others during the 1950’s 
and early 60’ s used mathematical thermal models of evolution of the moon to 
explore the consequences of certain origins (hot or cold accretion) , age 
(the same as the solar system, or younger) or composition (chondricic or 
otherwise) on the present state of the moon. These models predict a present 
day heat flow from the moon against which measurements could be compared, 
which provide a direct test of any of the hypotheses of evolution and com- 
position. Thus, the value of surface heat flow is a fundamental planetary 
quantity like mean density or moment of inertia, chat provide clues to the 
composition and internal distribution of properties in a planet. 

Determining the rate of heat loss from the moon did not have to await 
lunar landings. Theoretically, the average heat flow can be measured from 
earth, using radiotelescopes. Telescopic IR measurements from earth and 
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measurements of thermally genaracad microwave emission spacers, particularly 
as a function of solar phase, were ene of the r* 'aary sources of information 
on the detailed nature of the lunar surface, befere close-up photography became 
available from spacecraft. Except for the very near surface optical charac- 
teristics of the lunar surface, the thermal properties, temperature, thermal 
conductivity, density and heat capacity of the lunar surface layer were 
better known chan any other physical properties. Telescopic thermal obser- 
vations on the moon were the only source of information on the nature of the 
subsurface. 

Interpretation of microwave and infrared observations are strongly limited, 
by the uncertainty of the electrical properties of the lunar surf.ee, by 
the low level signals at long wavelengths, and by the uncertainty of the 
distribution of thermal properties with depth. Consequently, the relia- 
bility of estimates of heat flow based on these observations is not good, 
and the possibility of making in- situ determinations of surface heat flow 
on a landing mission would help resolve many of these uncertainties. 

It is not known exactly who first proposed that a measurement of the 
rate of heat loss from the moon should be made as part of the scientific 
effort of the Apollo Lunar Landing Program. During the early planning 
phase of Apollo, NASA sponsored many studies, some by scientific groups 
working in NASA, some by Industrial groups and some by academic scisrtists to 
recommend science programs that sight be carried out when the Apollo lunar 
landings were achieved. Reports of these studies usually contained. a 
'shopping list' of experiments and measurements that could be performed on 
the lunar surface. Some thermal measurements were usually Included on such 
lists, usually surface temperature and surface thermal properties. A few 
proposed trying to measure heat flow from the moon. 

My first exposure to the concept of lunar heat-flow measurements came 
via a phone call from Dr. Paul Lowman at the Goddard Space Center in Mary- 
land. Paul asked simply how difficult would' it be to measure heat flow 
through the surface of the moon. 

Soon afterward, in 1964, I was asked to join a small ad hoc committee tc 
address the problems of in situ thermal measurements on the lunar surface 
The membership of this cotmnlttee consisted of Arthur Lachenbruch 
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of Che U.S. G.S. in Menlo Pack, California; Sydney Clark of Yale University; 

Richard Von Herzen, then with Scripps Institution of Oceanography; and myself. 

It was a good mix; the two first mentioned were experienced with continental 
heat flow determinations, while Von Herzen and I were actively making heat 
flow measurements in the seafloor. The first meeting of this group took 
place in Washington. Sydney Clark’s first comment was that the only worth- 
while thermal measurement would be that of heat flow from the lunar interior. 

We all concurred and the remainder of this meeting was devoted to learning 
more about how and when the first Apollo missions would take place, the nature 
of the lunar surface, and the most feasible approach to measuring heat flow. 

Richard Von Herzen's term on the committee was short-lived as he accepted 
a position with UNESCO in Paris. Later in 1964, a conference in Houston at the 
Manned Spacecraft Center brought all the ad hoc committees representing many 
other scientific disciplines together. Lachenbruch, Clark, and I spent most f 

of our time at this meeting considering techniques for measuring lurar heat 
flow. The most straightforward approach to determine the flux through the 

t 

lunar surface layer or lunar soil is by direct measurements of the tempera- 
ture gradients and the thermal conductivity of the soil. Because ...e mea- 
surements would have to be made very close to the lunar surface, it was early 
recognized that the heat flow might be seriously disturbed by near-surface 
effects and, consequently, it was proposed that numerous (ten or more probes) 
be used to measure the flux. 

An alternative technique we considered was to set up heat flux meters on 
the lunar surface. These meters would consist of thin circular disks (or 
blankets) made cf material of very low but known thermal conductivity. If 
the ratio of the disk radius to thickness was SO to 1 or greater, the mean 
temperature dlffererce across the disk would provide a measure of mean heat 
flow. As in the case of the probes, several 'heat-flow meters' would have 
to be set up to measure and, it was to be hoped, average the local near 
surface variability. 

Sydney Clark was the principal author of a joint proposal to study methods 
for making measurements of surface heat flow or. the lunar surface. NASA gave 
a grant to Yale and Lamont-Doherty Geological Observatory of Columbia University 
to carry out these studies in response to this proposal. Arthur Lachenbruch 
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did aoc wish to coacinus active participc jn in Che proposed program 
because of che responsibilities he already had at U.S.G.S. However, 

* subsequent to the meeting at the Manned Spacecraft Center, he contributed 
a short report on che effects of surface temperature variations on measu- 
rements of heat flow. 

Under this grant, Sydney Clark carried out two analytical studies, the 
first addressed the problem of how to measure the very low level heat flow 
escaping from the moon in the midst of the large periodic variations of 
temperature at the surface produced by the moon's diurnal cycle, and the 
feasibility of using the heat-flow meter or blanket approach to measuring 
the surface heat flow. The results of this study were summarized in the 
final report of MSG-400 of Yale University (1965) * One principal conclusion 
was that the 'blanket' method of measuring lunar heat is not feasible because 
of the necessity of closely matching the local albedo and because a blanket' 
which would effectively alter periodic fluctuations, would take a prohibitively 
long time to come to equilibrium with steady-state surface, conditions. 

It is interesting to note that the blanket technique had a far more 
serious problem which was not fully appreciated until later. Llnsky (1966), 
Watson (1964) and others had pointed out the importance of radiative heat 
transfer in evacuated powders. The radiative component depends strongly on 
temperature and, as a result, behaves nonlinearly. A rectifying effect 
occurs that produces a net transfer of heat energy into the subsurface during 
the lunar day, "hich must be balanced by a steady mean outward flux. This 
effect, which is confined to the upper few cm, can induce steady-state fluxes 
that are ten or more times greater than the flux from the interior. Thus, a 
direct measurement of the flux at the surface would not have been represen- 
tative of the heat flow from the interior. Our subsequent measurements on 
the moon showed that this heat flux Induced by Che nonlinearity of radiative 
heat transfer in che lunar soil is very large. For example, a temperature dif- 
ferential of 40°C was Induced over a few cm of near-surface material at the 
Apollo 15 site. 

A second study carried out by S. Clark investigated techniques of mea-? 
surlng thermal conductivity with a probe inserted into Che lunar soil. The 
results are summarized in an Interim Report under Grant MGR -07-004-039 at 
Yale University. This report formed che basis of che technique that was 
eventually used to make in situ conductivity measurements on the moon. The 
method proposed consisted cf a cylindrical ring source of heat and a temperature 
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sensor placed within the heater, or at some distance from the heater, along 
the probe. A finite difference thermal model was used by Clark to investi- 
gate the problems interpreting the temperature change with time, in terms of 
the thermal conductivity of the surrounding material. This work showed that 
a reasonable determination of the thermal conductivity could be obtained by 
curve fitting, although it was clear that the thermal linkages and properties 
cf the probe and its contact with the lunar medium had to be well known. 

During the same period, scientists at Arthur D. Little, Inc. (ADL) of 
Cambridge, Massachusetts, in particular, Peter Glaser and Alfred Wechsler, 
had proposed a number of experiments to measure radiative and thermal pro- 
perties of the lunar surface as well as surface temperature as part of the 
Apollo program. ADL has been engaged in a study of the characteristics of 
the lunar surface layer under contract to NASA. In the spring of 1965, 

Peter Glaser visited Lamont-Doherty Geological Observatory to discuss ADL's 
Interest in working with us on lunar thermal experiments. This was Che 
beginning of a cooperation that survived through the entire program. 

1 965 Summer Studies & the LSEP Co nce pt’. During the summer of 1965, 

NASA and the National Academy of Sciences (NAS) sponsored tvo formative 
workshops, one at Falmouth, Massachusetts and the other at Woods Hole, 
Massachusetts. The purpose of both studies was to help NASA formulate a 
firmer set of scientific objectives and priorities for the Apcllo lunar 
program. S. Clark and I participated in the Falmouth study. 

Among the geophysical experiments proposed to be carried out on the moon, 
a direct measurement of heat flow had a high priority. One of the mo3t 
important concepts to emerge from the study was that of the Lunar Science 
Experiments Package (LSEP) . This concept recommended a single transmitter 
and power supply to be set on the lunar surface, to which various arrays of 
experiments could be connected. This concept, was accepted by the Apollo 
program planners and, with the addition of an 'A' for Apollo, became the 
ALSEP science stations that were eventually set at five locations on the 
lunar surface (see NASA Special Report, NASA SP-88, 1965). 

At these meetings, NASA management also selected experiment teams for 
those experiments that had beer given a high priority. The heat-flow 
experiment team named Included Marcus G. Langseth as Principal Investigator 
and Sydney Clark as Coinvestlgator. During the meeting, Gene Simmons, 
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a scientist interested in geothermal problems at Southern Methodist University, 
asked to be included in the team and was accepted. It was recommended by the 
NAS panel meeting at Woods Hole, that this team work closely with ADL scientists 
in developing the experiment. These decisions were the basis of a contract be- 
tween NASA and L-DGO to help in the development of a lunar heat-flow experiment 
for Apollo. 

One further outcome of the Falmouth meeting had an important effect on the 
subsequent technique to measure heat flow. It was recogrized that, to avoid the 
large surface fluctuations, it would be necessary to implant the probes to mea- 
sure the temperature gradient and thermal conductivity deep in the lunar soil. 
Because of the uncertain penetrability of the lunar surface layer, it was re- 
commended that a portable, astronaut-operated drill be developed to assure deep 
burial of the probes even in hard rock. Thus, the development of a hard rock 
drill (soon to be designated as the ALSD, Apollo Lunar Surface Drill), became 
a part of the Lunar Heat Flow Experiment. 

Organization for Flight Experiment Development 

By 1966 NASA was firmly committed to the ALSEF scientific stations as part 
of the Apollo lunar landing sites. An ALSE? Management Office was established 
at the Manned Spacedraft Center; experiment teams, each headed by a principal 
Investigator were established, and means of communication and coordination be- 
tween the growing number of elements in the program were developed. 

The first task was to select an ALSE? system contractor who would be res- 
ponsible for the overall management of the development, test ng, manufacture, 
flight qualification and flight support for the ALSEP system, from among pro- 
posers from major aerospace corporations. The field was soon narrowed to three; 
namely, the Bendix Aerospace Corporation of Ann Arbor, Michigan; Space General 
in Los Angeles, California; and, TRW, also in the Los Angeles area. It seemed 
as though there were an Interminable series of meetings at which ve listened 
politely to long engineering and management presentations of their development 
programs. In the process we learned a whole new language and we offered 
criticism and evaluation where we could. The eventual winner of the competition 
was the Bendix Corporation, Aerospace Systems Division. 

At the experimental level competition was equally vigorous and diverse. A 
subcontractor to develop each ALSEP experiment had to be chosen. Four companies 
were being considered to develop the Heat Flow Experiment, Arthur D. Little, 

Texas Instruments, Bendix and Teledyne Corporation. It was difficult at the time 
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co evaluate che proposals submitted by these companies, because the heat- 
flow experiment design was at che conceptual stage. There was no precedent 
for such an Instrument and It had to operate in che lunar subsurface envi- 
ronment which was poorly defined. Consequently, It was decided co go with 
the group chat had che greatest experience in original thermal design and 
knowledge of the characteristics of che lunar surface. The obvious choice 
was Arthur 0. Little, Inc. (ADD of Cambridge, Massachusetts. There was a 
protracted negotiation between che National Aeronautical Space Administra- 
tion, the Lamont-Doherty Geological Observatory of Columbia University, 

Arthur D. Little Co. Inc. and the Bendlx Corporation over the responsibility 
for the overall integration of the heat-flow instrument. I personally 
favored having Arthur D. Little Inc. assume that role, which meant chat the 
Arthur D. Little Inc. would deliver a complete instrument subsystem to be 
'plugged' into che ALSEP. I did not have my way. The final decision made 
che Bendix Corporation the principal contractor and Arthur D. Little a sub- 
contractor to Bendix, with a primary responsibility for instrument develop- 
ment and the manufacture of che heat-flow subsurface heat-flow probes, which 
were the primary components of the experiment. Gulton Industries Incorporated 
of Albuquerque, New Mexico was chosen to develop the electronic circuitry. , 
This award was given on che basis of a critical circuit design which automa- 
tically compensated for drifts in the offset and gain in the necessary am- 
plifiers and the A-D converters needed to measure temperature. The circuit 
design and its application to the measurement of temperature difference was 
primarily due to Jay Cox of Gulton Industries. 

The Dril l 

For most of the conceptual designs of che experiment a drill, capable of 
penetrating the lunar surface, would be needed. NASA decided to manage this 
development outside of the ALSEP contract. There were two companies vying 
for the contract; The Martin Marietta Company of Baltimore, Maryland and 
the Northrop Corporation of California. The Martin Marietta Company had 
earlier experience developing space tools for NASA, inducing a surface drill, 
and che Northrop Corporation was involved in a deep lunar drill development 
(deeper than 100 feet). The Martin Marietta Company won the contract. 

Nearly all of 1966 was spent on che process of reviewing proposals, 
making site visits and contract negotiations; in summary, selecting the 



industrial members of the heat -flow experiment development team. The 
complete team was defined by the end of the year. 
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Ill . CONSTRAINTS OF THE CONCEPTUAL DESIGN OF THE 
HEAT-FLOW EXPERIMENT 


The principal constraints on the conceptual design of the experiment were: 

1) The mechanical properties of the lunar surface layer; 

2) the thermal properties of the lunar surface layer and its heat budget; 

3) the weight and dimensional constraints of the spacecraft; and, 

4) the limitations of astronaut activity on the lunar surface. 

Since our knowledge of the subsurface was limited, we decided to make the 
experiment in such a way that reliable results could be obtained over as wide 
a spectrum of mechanical and thermal properties as pcsslble. What follows is 
a stannary of the status of information about these constraints when the design 
was started. 

Numerical Properties 

Some indications of the physical, characteristics of the lunar soli has been 
inferred from the thermal and optical respor.se of the lurar surface to solar 
radiation decades before the Apollo program. The facts that, at noon the moon's 
surface temperature was nearly that required for radiative steady state, and 
the rapid cooldown after sunset, or during an eclipse, indicated, that the surface 
layer had an extremely low thermal inertia (Kpc) 1/2 (Wesselink, 1948 and 
Jaeger, 1953). Microwave observations throughout a lunar day indicated that 
subsurface materials to a depth of several centimeters had a very low con- 
ductivity as well (Piddington and Minnett, 1949) . The highly insulating 
thermal properties of the lunar surface could be best explained if the surface 
material was finely powdered material in a vacuum. Photometric observations of 
low albedo, peaking of brightness at zero phase and low polarization also 
Indicated that the lunar surface was covered nearly everywhere by extremely fine 
powders (e.g. Hapke, 1971). 

Further verification of the dusty nature of the lunar surface came with the 
Panger missions. High resolution photographs shoved that topography at the 
smallest scale was rounded and smoothed by a blanket of materia] that pre- 
vailed everywhere except over fresh craters. The ranger photographs also 
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revealed that the lunar mare surfaces were peppered everywhere by craters up 
to 30 meters in size (see. Ranger VIII and IX, Technical Report # 32-800, p. 
161-170). This observation suggested the mare surface; had been hit every- 
where at least once by objects large enough to create craters twenty to thirty 
meters in diameter. Such craters would be from three to six meters deep, 
based on experiments in materials simulating these thought to exist on the 
lunar surface. If craters this deep, or deeper, had excavated material during 
a hypervelocity impact, then the surface layer of the moon would consist of 
fine fragments and rock debris to at least these depths. 

The Surveyor Missions provided the first eye level view of the lunar sur- 
face. These views showed several meter-sized boulders and numerous smaller 
rock fragments strewn over the lunar surface. (Surveyor Project Final Report, 
Part II. Technical Report 32-1265 JPL) . 

The Surveyor Missions also provided the first direct information or. the 
distrlbutior. of grain sizes in the near surface layer of the regolith. These 
distributions were found to vary from site to site; but, in general, only 11 
of the volume was comprized of fragments ten certimeters' or larger in 
diameter. (Figure III-l). 

The soil mechanics experiments which were flovn on seme of the Surveyor 
Missions cculd perform simple mechanical tests or. the soil to a depth of 
several, centimeters. The tests showed the soil to have little cohesion and 
2 . density of about 1.5 gm/cm^. However, trench digging experiments shewed the 
resistance of the soil to penetration increased at depths of five to ten 
centimeters . 

Thermal Properties 

Infrared and microwave observations of the lunar surface had provided much 
information on the thermal properties and temperatures of the regolith. 
Observations of the temperature drop of the lunar surface during a toral 
eclipse indicated the thermal inertia of the very near surface material is about 
7.7 x 10~^cal cnT^sec” 1 ^®!^ . At 200®K, this corresponds to a conductivity 
of 4 x 10”^cal cm”^sec"’^°K, or 1.5 x 10 cm ^®K \ assuming 0c * 0.16 cal 

_3 

cm ®K. The surface temperature decay curve during the lur.ar night hewever 
was best fit by a thermal inertia of 1.25 x 10 which corresponds to a 
conductivity at 200®K of 7.7 x 10”^cal cm ^sec ^®K or 3.2 x 10 \‘att 

cm _ ^*K~^. These results are extracted from a summary given b> Shorthill (1972). 
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The difference in the thermal inertia between the two types of observations; 
the eclipse vith an uaibral period of about three hours and the lunar night- 
time of about 350 hours, resulted from an increase in thermal conductivity 
vith depth. 

High resolution infrared scans of the lunar surface showed certain areas 
that remained relatively hot during an eclipse. These became known as hot 
spots and generally coincide with the location of large fresh craters. It 
is now generally agreed that these anomalies are due to higher abundance of 
solid rock fragments around these craters. The solid rocks have much greater 
thermal inertias (about 4.5 x 10” cal cm sec~' t, K”‘ L ) than the typical regolith. 
However, before the Surveyor program the nature cf these hot spots wes 
debatable. 

T he heat budget ; The thermal properties from remote measurements permitted 
calculations of the transient and steady state heat budget in the regolith. 
Bather conservative limits on physical properties were assumed for such cal- 
culations based on the high variability of the properties observed over the 
surface. 

Periodic variations, due to the lunation cycles, were of primary concern. 

At the surface they had a peak to peak amplitude of nearly 30C>*K, but vith the 

generally low thermal conductivities the resulting thermal wave propagating 

downward into the subsurface would be rapidly attenuated. An example computed 

by Clark (1965) is shown in Figure IIX-2). The model is a 20 cm layer with 

-3 -2 -1 -1/2 

a thermal inertia of 2.24 x 10 cal cm °K sec over solid rock, of basaltic 

composition. It is clear that the amplitude is reduced to orly a few degrees, 

20 cm below the surface. 

There is else an annual variatior in mean temperature due to the eccentri- 
city of the earth-moon systems' orbit about the sen. The amplitude of this va- 
riation is abcut 3*K peak ro peak; however, the attenuation factor is 3.5 

times lover, so that some variation due to this annual cycle might be observed 

at depths of a meter or more. 

Steady State Gradient Due to Heat Flow : The bases for estimating the value 

of the gradient in mean temperature vith depth resulting from heat flow from 
the lunar interior were: 
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1) Microwave observations of brightness temperature gradient and models cf 
the thermal evolution of. the moon, based on ccsmologically ressorahle abun- 
dances of heat-generating long-lived radioactive 1 sc topes. 

2) The range of estimates of thermal conductivity which takes into account 
the observations described above and also the possibility that a significant 
Increase in the total amount of sclid rock fragmerts below depths of ters of 
centimeters could significantly increase the conductivity there. 

Microwave observations were Interpreted by Troitsky (1965) to Indicate that 

heat flow from the moor, was similar tc that observed for the earth (i.e. 6.3 yW 
_2 

cm ). We tock this value as an upper limit. For the lower limit, we used e chon- 

dritic model of the moon, with a uranium abundance of only 11 ppb and an initial 

temperature near zero *C (sec, for example, Reynolds et al., 1972 , who gives 

_2 

e value of about 1.0 yw cm for such a model). 

-«? -l 

The range of conductivities of the regolith was taken to be 2 x 10 " W cm 
-1 -4 -1 -1 

*R to 4.2 x 10 W cm °R (see, for example, Langseth et al., 1968). The 
high upper limit reflected uncertainties as tc the conductivity below 20 cm. 

When these two ranges are combined, the range of expected gradients resulting 
from heat flow from the interior of the moon is 2.4 x 10 ~ a K/cm.to 0.314®K/cm 
cr, in terms of the temperature difference expected over a one meter depth 
interval, 0.24 *K to 31.4 *K. 

Steady State gradient Associated with the Radiative Heat Transfer in the 
Upper Regolith : Linsky, 1966 pointed cut the possibility of a substantial 
temperature gradient in the upper few centimeters cf the lunar surface material 
due to the strong temperature dependence of finely divided particulate material 
in vacuum. This dependence results from the prominent role of radiative heat 
exchange between the individual particles. Evacuated powders as a consequence, 
transfer heat far more efficiently at high temperature than at low temperature. 

This nonlinear behavior results in a net flux of heat into lunar regolith during 
the lunation cycle w.ilch must be counterbalanced by a steady temperature gra- 
dient to maintain the overall energy balance. The consequences of this behavior 
of lunar surface material are important to the design of the heat-flow ex- i 

perimert. Firstly, the mean temperature gradient in the uppermost layer (0- 
10 cm) will be due to the nonlinearity of heat transfer and rot tc the heat < 

flux frem the interior. As we noted earlier, this phcncmeron ruler out the 
use of a surface blanket ot wefer-type heat flow meter . Secondly, any change 
in the radiative properties of the near surface layer will cause e change in the 
subeurface temperature because there is a reasonable chance that the ratio of 
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radiative to conductive beat transfer will have been changed. 

It should also be noted that this effect would have Important consequences 
for the Interpretation of microwave observations from any planetary body that 
has a powdery surface and a very tenuous atmosphere 

S pace Constraints 

One of the most significant constraints placed on the experiment's design 
was the size of the compartment In which the Instrument was to be stowed during 
the flight to the moon. The maximum length of an cb;ect to be placed in the 
compartmer.t was about 63 cm. This meant that any probe we wished to emplace to 
depths greater than 60 cm had to be folded, collapsed, or taken apart to be 
stowed. 

The large accelerations the .instrument would encounter during launching 
cf the Apollo rocket, due to vibration and shock, were another important 
consideration. Most accurate temperature sensors suffer shifts in tl elr 
calibrations when they are subjected to large accelerations. 

Less demanding constraints, as far as design was concerned, were the: 
maximum weight of about 4.5 kg, the maximum of peak pover allowed which was 
9 watts, and the range cf temperatures that the experiment would experience 
90 ®K to 400 ®K. 

A stronaut Surface Activi ty 

The disturbance cf a dust covered lunar surface by an astronaut as ha is 
setting up experiments on the surface car. have a profound effect on the subsur- 
face thermal regime. The footprints and the dust kicked up by walking charge 
the surface radiative properties (the surface albedc and the IR emissivlty) and 
the thermal Inertia. These changes in surface properties could cause a change 
In surface temperature and consequently, a charge in heat flow through the 
subsurface. These changes take years to propagate dowr to the depth of a 
meter or more. 

Other design considerations had to take into accourt that the astronaut would 
be greatly accumbered by his space.suit. In particular, the gloves greatly 
reduced the manual dexterity as compared to the use of the. bare hand. Addi- 
tional considerations were: that the deployment of the experiment on the sur- 

face had to be accomplished in the shortest possible time In order to maximize 
the number of tasks tb“ astronauts could accomplish during tht-.ir brief stay. 
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Apollo program plans called for the heat-flow experiment to be flown on 
Apollo missions 13, 15 and 16. Mission 13 never landed on the lunar surface 
because of the explosion of the oxygen tank in the CSM fuel cell. After many 
months of preparation with James Lovell and Fred Halse, the commander and 
lunar module pilot on the mission, it was very disappointing to have our first 
flight instrument burn up as the lunar module reentered the atmosphere over 
Australia. In retrospect, a heat-flow measurement at the Fra Mauro site 
(the Apollo 13 target site) would have been very valuable. Later orbital data 
showed this region to have the highest surface abundance of long-lived 
radioisotopes; thus, a measurement there would have provided a valuable end 
point for a heat flow and surface heat production correlation. Secondly, we 
would have had an earlier indication of the problems of drilling the regolith 
and could have been better prepared for the later missions. 

Apollo 15 - Hadley Ril.le 

.The time was too short between the Apollo 13 and 14 missions to reconfigure 
the Apollo 14 ALSEP to carry the heat-flow experiment. Consequently, our 
first opportunity to emplace the experiment on the lunar surface came with the 
Apollo 15 mission to Sima Hadley. 

Both astronauts worked to set up the ALSEP station and array of experiments, 
however. Drilling the holes and emplacing the heat-flow probes was the task of 
David Scott, the mission commander. The drilling proceeded badly. After the 
first meter of penetration, Commander Scott reported that the drilling was 
getting 'stiff' and downward progress stopped. He managed to drill two 
sections of the borestem down to a depth of about one meter, and to continue 
an additional section of borestem would have to be added. To disconnect the 
drill power head frcm the drillstring required a counterclockwise twi3t of 
the head relative to the drillstring to release a collet over the tapered 
upper joint. However, there was no sufficient counter tcrque on the stem 
to release the collet. With great, difficulty Dave Scott managed to put a wrench 
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on Che stem. The struggle vas painful to watch, the spacesult made the 
stooping and bending movements extremely difficult. It was or.ly the physi- 
cal strength and persistence of Dave Scott that allowed him to accomplish 
the task. After adding two more sections, Scott managed to continue pene- 
trating to a total depth of 170 cm. By this time he was applying his full 
lunar weight to the powerhead and the drill's motor was stalling. At this 
point, it was decided to remove the powerhead and insert the first heat-flow 
probe, even though the bit was still one meter above the target depth. 

Again, it was extremely difficult to get the collet to release its grip on the 
upper taper joint; but, again Scott succeeded after a cumbersome struggle. 

After the emplacement of the first probe, it was decided to go ahead and 
turn on the experiment. This would allow us to see if the experiment was 
working properly before expending, the valuable time and astronaut effort on 
completing the drilling for the second probe. It operated perfectly. 

The second drill hole was started, but immediately the same difficulties 
were encountered. The science team in the ALSEP room decided to call a stop 
to deployment of the heat-flow experiment until the second EVA. - The power- 
head was left on the drlllstem after penetrating about one meter. 

Between EVA's there were hasty conferences at mission control to try to 
diagnose the problems cf the drill based on performance. The concensus of 
opinion was that Scott was applying too much downward force on the drill and 
not allowing the drilled material to clear out of the hole. Before the 
second EVA Scott was advised to try to hold the drill back, tc actually pull 
up if necessary to hulp clear the cuttings. 

Soon after the second EVA began, Scott resumed the task of drilling. He 
tried to restrain the downward progress of the drill, but commented that the 
stem seemed to be augerlng deeper into the regolith. At one point in the 
drilling, there was a noticeable release of torque and Scott commented that 
it had freed. Later, analysis of the films and the photographs showed that 
the drillstem had actually parted at one of the joints below the surface and, 
therefore, the stem attached to the powerhead no longer had a bit. Even without 
the bit, the drill continued to penetrate twerty to thirty additional centi- 
meters and then stalled. We decided to halt the drilling at that point. 

Scott had the same difficulty removing the powerhead as before, but finally 
managed it. The second probe was inserted, but when it reached the bottom 
of the hole, the cop of the probe still protruded. Commander Scott felt chat 
there was some obstruction in the hole, but a rapid miscalculation on earth 
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Figure IV- 1 


Temperature histories of the sensors on probe i during the first 300 hours after emplacement. 
The subsurface geometry of the probe and temperature as a function of depth are shown on the 
left. Temperature as a function of time is shown on the right. 










30 


assured him that everything, was as expected. 

It* was abundantly clear that the drill that we had designed to penetrate 
the lunar regolith was not up to the task; nonetheless, we had thermometers 
buried in the subsurface. The deepest was buried to a depth of 137 cm. This 
depth was sufficient to assure an accurate measurement of the subsurface 
temperature gradient. In Figures IV-1 and IV- 2 we show the evolution of 
the temperature versus depth profiles measured by the probes, for the first 
300 hours after emplacement. Notice that, the initial temperature profiles 
cf both probes decreased with depth; but, as the probe slowly equilibrated 
with the surrounding regolith, a gradient of nearly 1. 7°C/m was established 
in the lower 70 cm of the hole. It took more than a month for the probe to 
fully equilibrate, which attests to the extremely low thermal conductivity of 
the regolith . A feature of the early temperature profiles of probe 1 was 
very pertinent to the drilling problem. The thermometers adjacert to the 
tapered joint show high temperature indicating a high frictional heating at 
this point. 

T he Deep Core Sample 

The lunar drill was also used on Apollo 15 to obtain a deep core about 
two and a half meters long. Scott again operated the drill and it penetrated 
the regolith quite readily. The difficulties came when the drill was to be 
extracted. Scott simply could not get the c*rill to pull up despite applying 
power; which, if the drill were held steady, would help clear cuttings out of 
the hole. Instead, the drill tended to twist Itself farther down. Finally 
with great effort and Jim Irwin's help the two pulled the drill up. 

Later analysis of the core samples showed that the material at depths 
greater than two meters had lain undisturbed for a half billion years. 

Diagnosing the Drill Problem 

The problem with the drill was simple. The subsurface of the regolith 
was much denser tban we had supposed. The regolith, at least to a depth of 
a few meters, is mainly composed of extremely fine dust of rocky and glassy 
fragments; yet it is packed so densely that it cannot be further compressed 
without great effort. This character of the lunar soil layer required that 
all of the drill cuttings pass up alongside the drill stem to the surface. 

The bulges in the stem at each joint prevented this; and, consequently, the 
drill was augerlng its way Into the moon as Dave Scott had noted. We had 



mistakenly assumed, based on models of the lunar surface layer, that the 
cuttings could be easily compressed into the adjacent material. 

The poor performance of the drill on the Apollo 15 mission could be 
completely explained by this hypothesis and It was demonstrated by running 
tests in a test bed of highly denslfled basalt powder. During these tests, 
the drill would penetrate beautifully until a joint entered the soil. Then 
the penetration slowed and eventually the motor stalled. 

The problem was solved by changing the design of the joint so that there 
was no bulge and so that the spiral fluting that moved the cutting up the hole 
continued through the joints. To do this, we had to bond titanium joints 
onto the ends of the composite borestems. At the same time, wc lengthened 
each section so that fewer joints were required. We adjusted the lengths 
so that no thermometer was Immediately adjacent to a titanium joint. The 
redesigned drill stem drilled the high density basalt powder model easily. 

A redesign of the wrench that decoupled the powerhead frem the borestem 
was also needed to' make it easier for an astronaut to use. To overcome 
the problem of retrieving the deep drill core an extractor, with a greater 
mechanical advantage, was designed. All of these modifications were made, 
tested and flight rated In the eight month Interval between Apollo 15 and 16. 
Apollo 16 - The Highland Site 

Charles Duke, che lunar module pilot on the Apollo 15 mission was active 
in the redesigning and retesting of the drill. Being taller than John Young, 
the commander of Apollo 16, Charlie was assigned the drilling task, since It 
was easier for him to manage the long first section. Charlie was a perceptive 
engineer and helped greatly with all the phases of the modifications to che 
drill. 

The drilling went perfectly for the first heat-flow hole. Each section 
sank slowly, but surely into the regolith. However, as Duke was inserting the 
first probe, John Young who was busily engaged in setting up ocher experiments 
of the ALSZP array, caught the flat cable that connects che heat-flow experiment 
to the ALSEP central station over his boot. Before he felt che tug of che 
cable, or before any word of warning could be relayed from the gi md, the 
cable snapped rather easily at the connector that plugged into the central 
station. 
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John Young apologized dejectedly when he discovered the heat-flow cable 
had become completely torn off. Charlie Duke grumbled about the time wasted 
in drilling the hole and inserting a probe that would never be used. Upon 
closer inspection, the astronauts found that the cable had torn off cleanly 
and squarely right at the outermost clamp of the astromate plug, the con- 
nector to the ALSEP transmitter (see Section III). Strain relief had not 
been provided for a tearing motion. The outer rim of the damp from which 
the cable emerged had a sharp square corner that undoubtedly helped in starting 
the tearing action. The tensile strength of the cable was sufficient to 
withstand almost any imaginable tug that an astronaut might apply. Because 
of the soft nature of the lunar surface and the poor downward visibility cf 
the astronaut that the irelmet afforded there was a good llkelyhood that they 
could trip ever one or more of the many wires that radiated out from the 
central station. The cables and plugs were designed to withstand this even- 
tuality. It was another design failure that completely negated all of the 
heat-flow team's preparation for the mission. 

There was a faint ray cf hope. The cable was torn so cleanly that there 
was a possibility it could be repaired. .Thus, between the first and the third 
EVA's we began a nonstop effort to discover a way to remake the twenty or 
sore electrical connections. If the 'Astromate Connector' could be disconnected 
from the central station, if the probe could be pulled out of the borcstem 
without damage, if the whole experiment could be bagged and carried back into 
the LEM, if the tools existed to open the clamp, if the insulation or the 
flat cable could be removed to expose the surface of the ribbons of copper 
conductor without damaging them, if the cable could be inserted back into 
the plug so that most of the contacts were good, and if the repaired experi- 
ment could be set top on the lunar surface, the heat flow experiment could 
be saved. 

6 P° 93 ^ble Repair of the Apollo 16 Accide nt 

Developing a fix was an all-out effort. NASA assigned its 'Tiger Team' 
to the problem. The principals were Dor. Arabian and Fred Haise. Haise donned 
his space suit to determine whether the 'Astronaut Connection* could be 
disconnected from the central station. It could be done if the station housing 
were tilted backward and the astronaut virtually got down on all fours. The 
rest of us worked on ways of stripping the Insulation away from the conductors 
on the ribbon cable. Any sharp edged tocl could be used as long as one were 
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patient enough. However, it was not certain what tools the astronauts had 
on board the lunar module. Each astronaut carried a discretionary packet of 
things not inventoried. There was some hope that either Duke's or Young's 
packet contained a sharp edged tool. Another possibility wss co use one of 
the lunar rocks. The harder breccias were an ideal abrasive, so that they 
could be used much like a coarse whetstone. By rolling ribbon cable over a 
hard mandrel and rubbing firmly, but carefully, with a flat side of a stone 
the insulation could be removed with only slight damage to the copper conductor 
inside. 

By 7 a.m. the next morning, a procedure had been written and Don Arabian 
was calking Fred Haise threugh a simulation cf the repair. After the Simula* 
tion the repaired plug was tested and most, but not all, of the contacts were 
good. However, enough were good to give useful temperature data thanks to 
considerable redundancy in the use of the conductors. The amazing thing to 
me was that, while the enthusiasm and energy of most of us seemed to flag 
through the wee hours of the morning, Arabian and Haise seemed to gain in 
strength. By 8 a.m. they were working with the vigor most of us had started 
with twelve hours earlier. I later learned that Arabian had been up the two*’ 
preceding nights working on the problem of the errant thruster controls in 
the command module that nearly cancelled the Apollo 16 lunar landing. 

At about 9 the next morning, an ad hoc beard consisting of top members 
of the Apollo Program Directorate convened. Arabian presented the procedure 
for carrying out the repair. There were many attractive reasons for trying it. 
Firstly, it would be a convincing demonstration of the valuable role men play 
in space missions. The repair would have strong public relations appeal, 
especially the possibility of using lunar material tc do it. The heat* flow 
program had the highest priority of all the ALSEP experiments, sc that one 
of the principal objectives of the mission would have been saved. This was 
tbe last chance to get a heat-flow measurement in the highland areas during 
the Apollo Program. 

But, the reasons for not trying the repair were stronger. To gather up 
the experiment at the end of EVA 2 would cut about a half hour of time from 
other parts of the science program. Another hour at the start of EVA 3 to set 
the repaired experiment back up would have precluded the trip to North Crater, 
a primary geological target of the mission. The probability of a successful 
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repair was not very high, so chac all this tine could conceivably be wasted. 

Most importantly however, was the fact that the repair to the cable had to 
be done inside the lunar module and it would cut into the rest time of Charles 
Duke and John Young between EVAS 2 and 3. After EVA? they would have to 
perform the critical launch from the lunar surface and rendezvous without 
rest. 

On balance then, it was clear that the correct decision was to forsake the 
heat-flow experiment. This decision Rocco Petrone made unequivccably after 
lengthy discussion and consideration of all sides. 

Apollo 17 - Taurus Littrow 

In contrast to the Apollo 15 and 16 problem the operations involved with 
setting up the experiment at the Apollo 17 site of Taurus Littrow went smoothly. 
Each borestem was drilled to its prescribed depth and the probes installed. The 
operation of the experiment system was perfect and it continues to work as of 
1976 «^ith virtually no failures. . The only degradation in performance is noise 
that appears on the tnermocouple channels during the lunar mcrning (see 
Section V) . 

By the end of the Apollo Program, there were two working hear-flow expe- 
riments operating on the lunar surface, one at Rlma Hadley with the probes 
burled to shallower than planned depths, and the other at Taurus Littrow, 
which is set up as planned. The shallow depth of burial at Apollo 15 proved 
eventually to be beneficial. The four and a half years reccrding lunar 
temperatures at relatively shallow depths has provided an excellent record 
of lon£,-term temperatures in the moon's subsurface. The first major failure 
in the Hadley Rille experiment occurred in December of 1975, four and a half 
years after installation. 
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V. EXPERIMENTAL RESULTS 

A. SURFACE TEMPERATURES AND THE NATURE OF THE SHALLOW REGOLITH 
1 . Surface Temperature Deductions 

Ac each of Che haat-flow stations, sections of the electronics" cable 
containing one or more thermocouples of + 0.5 # K accuracy were deployed 
above the lunar surface. These thermocouples are in radiative balance 
with the lunar surface, Che solar flux and space and hence provide a 
measurement of the surface brightness temperature. At the Apollo 15 
site, a total of five thermocouples Ctvo at the probe 1 station; three 
at probe 2) were exposed above the lunar surface between the protruding 
borestem and the heat-flow electronics' station. Because of poorly defined 
orientation and an unknown amount of lunar dust coverage, interpretation 
of the daytime thermocouple data in terms of surface temperature could not 
be accomplished with less than 20 e K uncertainty. However, since daytime 
surface temperatures are primarily controlled by the solar albedo, theo- 
retical models taking into account latitude and previous albedo estimates 
were used to supplement the nighttime data and calculate the mean surface 
temperature over a lunation. 

At the Apollo 17 site, each of the probes was deployed to the nominal 
depth of 2.3 m, resulting in the placement of two above surface thermo- 
couples (one at each station) , approximately 15 cm from the cop of each 
borestem. Improved determination of orientation allowed noontime tem- 
perature estimates of +7®K accuracy. 

Fortunately, more accurate measurements of surface temperature could be 
made during the lunar night, when cable orientation and surface solar 
albedo can be neglected. The cooldown response of the surface in the 
absence of solar flux is controlled primarily by the thermal properties 
of the upper 5-15 cm of the regolith and can, therefore, be utilized to 
deduce a thermal property profile for Che near surface layer. 

In the absence of solar flux, the radiative balance of each thermocouple 
can be expressed as: 
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were a ■ cable radius, 

pc ■ volumetric htac capacity of cable, 

c^ r « infrarad absorptivity of cabla, 

e - tmissivity of cabla, 
c 

e ffl ■ emissivlty of lunar surface, 

F * view factor of thermocouple element to the surrounding topography, 

ff * Stefan-Boltzman constant; and, 

T ,T • cable and lunar surface temperature respectively, 
c m 

Except for periods of rapid temperature change (such as during an eclipse, 
or immediately following lunar sunset), the finite time constant of the cable 
may be neglected and Equation V-A-l may'tre reduced to the following expression 
for T 
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The primary sources of error In calculation of surface brightness temperatures 
using Equation V-A-2 derive from uncertainties of the radiative properties 
of the cable and topographic effects. During eclipses and immediately fol- 
lowing lunation sunset, errors in the determination of 6T c /5t must also be 
taken into account. 

Radiative properties of the thermocouple cable, e and a . , as well as 

c cir 

the term accounting for the cable r s thermal mass, aPc, were determined ex- 
perimentally by Clifford Cremers of the University of Kentucky on a sample 
piece of thermocouple cable. Values of ape were also checked by an evaluation 
of cable component properties. Values of £ c *nd o cir were adjusted to account 

for partial cable covering by lunar fines as was indicated by the different 
temperatures reached by the five thermocouples at the Apollo IS sice at the 
end of lunar night. More detailed descriptions of the surface temperature 
analysis can be found in Kelhm et al. (1973) and Kelhm and Langseth (1973). 

2 . Temperature Dependence of >- he Surface Layer thermal Conductivity 

A major result of the surface temperature measurements at Hadley Rille and 
Taurus Littrow was the calculation of the mean lunation temperatures: 207 
+ 7*K at Apollo 15 and 216 + 5*K at Apollo 17. The lower mean temperature at 
the Apollo 15 site is due primarily to the early morning shadowing effect 
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produced by a topographic high to the east of that site. The mean surface 
temperatures observed are 35-40 *K lower than the equilibrium temperatures 
measured at meter depths. Only a small part of this mean temperature rise 
(no more than 5*K) can be accounted for by the temperature gradient produced 
by the measured heat flow. The mean temperature rise is due mainly to the 
non-linear effects of radiative heat transfer within the highly porous dust 
layer about 2-3 cm thick at the surface. During the warm lunar day heat -s 
transferred more effectively into the lunar surface than it can be trans- 
ferred out during the cold lunar night. To conserve net flux over a lunation, 
a mean-temperature gradient is established mainly confined to the porous 
surface layer. The size of this gradient is largely a function of the relative 
contributions of the radiative and conductive components of the effective « 
thermal conductivity. Because most of the total mean- temperature difference 
takes place in the upper few centimeters, the relative contribution of the 
radiative and conductive terms can be determined unambiguously only in this 
region. If one characterizes the conductivity temperature dependence by the 
Linsky (1966) parameter, 

>350 ' > r (350)/K c 

where is the radiative component, and K c is the conductive component, the 
observed mean temperature gradients at Hadley Rille and Taurus Littrow require 
H 35 O * 2.5 - 3.0 and R 350 * 1.7 - 2.2 respectively. These values indicate a 
somewhat larger radiative component then that inferred from measurements of 
returned lunar fines (Cremers and Birkebak, 1971 and Cremers and Hsla, 1973) 
and suggest that thermal properties in the low density surface layer are a 
strong function of in situ grain-contact configuration. 

3. Thermal Response During the Lunar Night 

As discussed earlier, the surface temperature response in the absence of 
solar radiation is primarily a function of the shallow regolith thermal 
properties: p, the density; c, the specific heat; and k, the thermal con- 

ductivity. Variations with depth and temperature can also be constrained 
by the nature of the surface cooldown. Matching of the postsunset surface 
temperatures was reduced to the selection of an appropriate conductivity 
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function, k(z,T), by utilizing ocher data to constrain the density profile and 

temperature dependence of the specific heat. Density profile selection 

appropriate to each sice was based on soil mechanics data provided by J.K. 

Mitchell (personal communication) for Apollo 15 and D. Carrier (personal 

communication) for Apollo 17. The specific heat as a function of temperature 

was modeled by a least square fit to the data of Robie et al. (1970) on 

Apollo 11 fine samples. Similar measurements on Apollo 15 soil sample 

15301,20 and Apollo 16 soil sample 60601,31 by Hemingway et al. (1973) yielded 

nearly Identical values of specific heat versus temperature as the Apollo 11 

sample. Since the ratio of the radiative to conductive component of the 

’•hermal conductivity was constrained by the observed mean temperature gradient 

discussed earlier, only the conductive component, K , remained to be varied 

as a function of depth to fit the cooldown surface data. One dimensional 

finite difference models were utilized to fit the nlghtime data. Figure 

V-A-l shows the fits to the nighttime surface cooldown data along with the 

derived mean conductivity profiles (solid curve in inset). In both the 

Apollo 15 and 17 models, a low conductivity layer of about 2 cm thickness, 

with 210°K conductivities between 0.9 and 1.5 x 10~^W/cm r °K was required to 

fit the steep drop in surface temperature for the first 40-60 hoi s following 

sunset. The Apollo 15 model then required a steep, but not discontinuous, 

-4 - 

rise in conductivity with depth, to a value of ^ 0.8 x 10 w/ cm ®K at a 
depth of 10 cm, to produce the increased flattening of the surface cool- 
down curve through the lunar night. The Apollo 17 mo .'el, however, required a 
much sharper jump in conductivity at a depth of 2-3 cm to produce the abrupt 
flattening of the cooldown curve at 10 s phase past sunset characteristic of 
the Taurus Littrow site. The required increase in conductivity with depth 

from 3-15 cm is slight enough that the Apollo 17 model may be considered 

-4 

essentially a two layer model. Conductivity values cf 1.2 - 1.3 x 10 are 
required below 5 cm. For both the Apollo 15 and 17 sites, the conductivities 
inferred at depths of 10-15 cm are in excellent agreement with the values 
deeuced from the annual wave data for the deeper resolith. 

The principal advantages of the heat-flow experiment surface temperature 
measurements was the ability to measure local temperatures immediately 
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Surface temperatures during the lunar night at the Apollo 17 (solid docs) 
and Apollo 15 (open circles) heat flow 3ltes. Vertical bars are estimated 
errors. The continuous curves are theoretical curves derived from the ther- 
mal property models shown in the inset. 


Figure V-A-l 
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following sunset. The conductivity models inferred from the data are in 
essential agreement with earlier models based on earth-based infrared 
brightness temperature measurements. However, the greater detail of thermal 
property determination could only be accomplished with the addition of precise 
postsunset data. Measurements of + 2*K accuracy within the first 10-20 hours 
after sunset are critical to an unambiguous modeling of the near surface 
thermal property profile. For example, as shown in Figure V-A-l, slight (2®K) 
changes in the nightime cooldown curve within 40 hours following sunset can 
lead to relatively large changes in the derived conductivity profile 
(compares curves 17 and 17a) . 

4. Thermal Response During Eclipse 

Because of the relative short duration, surface cooling during periods of 
partial or total solar eclipse are controlled primarily by the first centi- 
meter of the surface layer. At the Apollo 15 site, surface temperatures were 
computed from the thermocouple data for the complete eclipse of August 6, 1971. 
Umbral surface temperatures of + 5®K accuracy were obtained. (The decrease 
in accuracy relative to lunation nighttime temperature is due to uncertainties 

in the calculation of 6T /6T which must be included during the rapid tempe- 

c 

rature variations of an eclipse; see Equation V-A-l ) . The conductivity of the 
upper centimeter was varied in the one-dimensional models until agreement with 
the umbral surface temperatures was obtained. The preeclipse temperature pro- 
file was computed from the thermal property profile deduced from the lunation 
cooldown and has little effect on the eclipse deductions. Surface temperatures 
and feasible conductivity model fits for the August 6, 1971 eclipse at the 
Apollo 15 site are shown in Figure V-A-2. The resulting best fit conductivity 
function for the upper centimeter yielded: 

K (T) - 0.67 x 10' 5 + 0.42 x 10’ 12 T 3 

3 

indicating at 210® conductivity 1.05 x 10~ 3 w/cm~®K, in agreement with the 
lunation results for the near surface layer. The accuracy cf deduced umbral 
surface temperatures led to an uncertainty of + 20% for the conductivity fits. 

A similar analysis in a total eclipse at the Apollo 17 site also yielded a 
surface conductivity in agreement with the lunation cooldown data. 
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Figure V-a-2 


Surface temperature history during the umbral phase of the lunar eclipse of August 6, 1971. 
The vertical bars span a range of two standard deviations centered at the reduced surface 
temperature measurement. The solid lines represent computer generated models characterized 
by conductivity functions, ’•* (T) , for the t. ? centimeter. 



5 . C omparison With Laboratory Measurements 

Laboratory maasuramants of tha thermal properties of returned regolith 
samples provide a direct comparison with tha in situ measurements from tha 
heat-flow sites, the most striking result of these measurements both for small 
samples of lunar soil and intact core tube samples, is the rather narrow range 
of conductivity values obtained. Small soil and sample measurements (Cremers 
et al. , 1970; Cremers and Birkebak, 1971; Cremers, 1972; and Cremers and Hsia, 
1973), utilizing the line heat source technique, yielded conductivities in 
the range 0.5 -3.0 x 10~^ w/cm~*K over a temperature range between 100 and 
400°K. All samples were characterized by increasing conductivity with tem- 
perature. The relative contribution of the radiative term to the effective 
conductivity was characterized by R^q [- (350) /K] values between 0.5 and 

1.5, approximately a factor of two less than that deduced in situ. Most 
of the measurements were done at low bulk densities ( ''■•1.3 gm/cm^). 

However, measurements at densities as high as 1.9 (characteristic of the 
deeper regolith) did not yield significantly different results. 

Conductivity measurements of unopened core-tube regolith samples have been 
made using a new radiative flux technique (Horai et al. , 1977). By conducting 
the measurements with minimal disturbance to the cores, it was hoped that the 
thermal conductivity characteristic of in situ regolith could be obtained. 
However, for the two Apollo 17 drill cores 70002 and 70006, preliminary results 
yielded conductivities in the range 1.8 - 3.6 x 10~^ w/cm”*K, in close 
agreement with the line source laboratory results. The Apollo 17 drill- 
core samples contained material from 1-3 m depths below the lunar surface with 

3 

bulk densities of 1.74 - 1.30 gm/cm . The large discrepancy between the 
laboratory measurements and the in situ results at depths greater than 10 cm 
suggests that the deeper lunar regolith possesses an intragranular thermal 
contact that has been enhanced by mechanical stress over millions of years 
of metecritic impact. Only the near surface layer (0-2 cm), which has been 
overturned on a much more rapid time scale, exhibits the low conductivities 
obtained for the returned samples. At greater depths, whatever long-term 
mechanical process is involved, conductivities 3-5 times the laboratory results 
occur. The annual wave and transient deductions, as well as the short term 
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in situ conductivity experiments and surface temperature cooling histories, 
could not be accounted for by a regolith characterized by the laboratory 
conductivity results at depths greater than 10 cm. These results emphasize 
the importance of in situ measurements for reliably assessing the near- 
surface thermal regimes of other planetary bodies. 
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V-B. OBSERVED TEMPERATURE HISTORIES IN THE SUBSURFACE 
AND ON THE SURFACE 

Temperature sensor locations and final subsurface emplacement geometries 
are shown in Figure V-B-l. The accuracies and resolution of each of the sensors 
were given in Section III-B. Figure V-3-2 presents a three and one he. If year 
history of temperatures at some representative sensors. Fcur distinct temporal 
components are evident. 

1) We see initial cooldown of the probe and borestem after emplacement 

of the probes showed that the initial temperatures were 295 - 305 °K. Cooling 
to ambient regolith temperatures is rapid at the start; however, full equili- 
bration requires nearly three lunations as Figure V-B-2 shows. 

2) We see a diurnal component on sensors within 80 cm of the surface. 

Below 80 cm, the nearly 30C>°K peak-to-peak surface variation is attenuated 
to a level below the noi3e of the absolute temperature measurement. Close 
examination of the Apollo 17 sensors shows a small variation with a diurnal 
period. This variation is not due. to heat transfer downward 'in the borestem, 
but Induced by the electronic circuitry which is apparently slightly tempe- 
rature sersltive (see Section III-B) . Two notable features of the diurnal 
variation are the high frequency contert of the variations and the fact chat 
they are nearly in phase with the surface variation. These features indicate 
that the predominant mode of heat transfer is radiative transfer along the 
borestem. 

3) We see an annual temperature variation, which is most clearly seen by 
the Apollo 15 thermometers at 45 and 59 cm. As will be discussed in detail 
later, this variation is propagated downward through the regolith and its 
amplitude and phase lag are rot strcngly affected by the presence cf the 
borestem and probes. 

4) All sensors show a slow increase in temperature with time. The 
shallowest sensor begins to detect this increase after a few lunations, 
whereas the deepest sensors, for example those at Apollo 17 deeper Chan 150 cm, 
begin to detect the rise after a year. This component of the temporal 
variation is due to an increase of a few degrees K in surface temperature 
around the borestem caused by the astronaut’s footsteps. Notice the sub- 
stantial difference in this effect between probe 1 and probe 2 at the Apollo 

17 site. 
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Configuration of the heat-flow probes in the lunar regolith. 
Additional details about the heat-flow instrumentation may 
be found in Langseth et al. (1974) 
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Figure V-B-2 


Subsurface temperature histories covering a 3. 5 yr period of representative sepsors 
on the four heat-flow probes. Initial cooldown of the probes after Insertion can be 
seen at the start of each history. The small spikes in the first 60 days result from 
the short-term conductivity experiments. 
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In Figure V-B-(3-14) are shown temperature records of principal sensors 
Che first and last lunations observed on the surface. Figure V-B-l (a & b) gives 
the sensor designations used to Identify each curve. The records shew many 
interesting details. 

At the Apollo 15 site five of the thermocouples are in parts of the cables 
that are hanging above the lunar surface and are exposed to the full solar flux. 
The sixth, TC14, is inside the part of the borestem projecting above the surface. 
During the first lunation, shortly after the experiment is turned on, there is 
an abrupt change In temperature of all exposed thermocouples and TREF (which 
is the sensor on the thermal plate of the electronics) . This results from 
realignment of the electronics and the probe I cable during EVA 2 and the 
insertion of probe 2 which was done on EVA 2 (see Section IV) . 

An eclipse occurred during this lunation which is clearly seen by all 
thermocouples. The abrupt fall of temperatures is repeated at sundown. There 
is a pronounced decrease in cooling rate at about 140°K, which is due to radiation 
frem surrounding topography which is still sur.llt. The slow warming of the 
thermocouples before dawn is probably due to a similar effect. However, at 
Hadley Rille the topography west of the experiment is much lover than that east 
of tfce experiment which explains the difference in magnitude of the effect 
between post sunset and predawn parts of the record. 

Relatively low temperatures observed by TC12 and TC13 at probe 1 during the 
morning are due to partial shadowing by the very small scale surface roughness 
and the fact that the cable is nearly aligned with sclar radiation in the 
morning. Similar effects are seen to reduce temperatures at probe 2 thermo- 
couples during the evening. Large dips in temperature occur when the borestem 
shadow crosses the section of cable containing the thermocouple. 

The records of subsurface thermometers, during the first lunation at 
Apollo 15, show in more detail the initial cooldown of the probes and the 
diurnal variation of temperature. The topmost sensor TG 11A actually shows a 
drop in temperature during the eclipse. This is because the sensor is radia- 
tlvely coupled to the top of the borestem where the temperature drops nearly 
200*K during the eclipse. The affects of two thermal conductivity measure- 
ments made late in the lunar night of the first lunation are alsc noticeable. 
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The records frcm che 42nd lunation ac Apollo IS shew chat an eclipse also 
occurred. Notice that many of the small details noted three and a half years 
earlier are still repeated. 

The records frca the Apollo 17 probes show quite different temperature 
variations during a lunation because of the deeper probe depths. Only one 
thermocouple at each probe, TC12 and TC22, are exposed above the surface. 

The cables emerging from the top of the borestem were purposely aligned north 
and south, resulting in a far more symmetrical temperature history during the 
lunar day. 

TC 13 and 23 are buried about IS cm below the lunar surface (see Figure V-B- 
lb). A 50°K peak to peak diurnal variation in temperature is observed. 

This is much greater than would be observed if the periodic surface variation 
was propagating downward by conduction alone. These sensors have a relatively 
strong radiative link with che top of the borestem which is further evidenced 
by the sycchronousness of the variations with solar phase. Both temperature 
reccrds on thermocouples at the IS cm level are asymmetric relative to solar 
noon, but each in an opposite sense to the other. The reason for this is cot 
known for certain, but is believed to be related to the different orientations 
cf a reflecting shield placed on tep of the borestem. 

The thermocouple at 6S cm below the surface records a variation that is 
barely preceptible at the scale of che plots. In large part this variation is 
an artifact. 

The source of this artifact has bean traced to c.opper/Kover (Cu/KO) 
junctions in each thermocouple’s electronic circuit, which are meunted on a 
circuit beard in the electronic's housing. The areas of the circuit boards where 
these junctions are located have variable temperature during a lunation cycle 
and, consequently, introduce a varying EM? at the input to the data sensing 
circuitry which is extraneous. This source of error in the thermocouple 
measuring circuits has been thoroughly analysed and is reported in a Bendlx 
Corporation Aerospace Systems Division Technical Memorandum ASTIR TM30/31 by 
B.D. Smith. This analysis allows the thermocouple temperature to be corrected 
so that absolute temperature can be determined to + 0.26°K. The data shown in 
the figures has not been corrected. 

Another problem with the Apollo 17 thermocouple observations is seen 
in the 25th. lunation data. A large + 15®K noise appears only during the 
lunar morning, starting after the solar phase reaches 30°. This obviously is a 



thermally Induced noise but Its source has not been ascertained. This noise 
appeared during the 2nd lunation, and appears at a lower level during the 1st 
lunation. 

The records of subsurface temperatures at Apollo 17 are monotonous com- 
pared to the Apollo 15 observations. The initial cooldown is clearly seen 
during the first lunation; subsequently, the only discernible variation is 
due to electronically induced noise. 
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V-C. THERMAL PROPERTIES OF THE DEEPER REGOLITH 


Through fortuity as well as design, measurements of the surrounding 
regollth thermal properties were provided by analysis of four distinct 
components of transient temperature histories recorded by the subsurface 
sensors of the heat-flow experiments. Preliminary estimates were provided 
within the first sixty days following the experiment emplacement by analysis 
of the relatively short-term transients produced by probe regollth equili- 
bration and the activation of small probe heaters. Because of the small 
amount of heat dissipation Involved, the early measurements characterized 
the thermal peopertles of primarily the first few centimeters of material 
adjacent to each borestem. After two or more years of data were available 
at each probe site, analysis of long-term temperature variations, charac- 
teristic of a much larger volume of surrounding regollth material, were 
used to refine and revise the early thermal property estimates. The long 
period transients resulted from annual variations In the solar flux at the 
lunar surface and reequilibration of the local subsurface thermal regime 
following astronaut- Induced alteration of che thermal and radiative pro- 
perties of the surface layer. 

The present section is devoted to a description of the methods of analysis 
utilized for the regollth thermal property estimates based on four distinct 
components of transient temperature variations observed at the heat-flow 
stations. Particular emphasis is placed on details of che long-term analysis, 
upon which the final heat-flow values were based. At che end of che section, 
a comparison and discussion of the different thermal property deductions is 
presented. 

1 . Short-term Measurements 

a. Initial Probe-Borestem Equilibration 

When the borestem was drilled into che soil and che probe inserted, both 
were at temperatures considerably higher than their subsurface surroundings. 

The higher temperatures resulted from the temperature of these components above 
the surface before emplacement and from heat generated during drilling. Because 
of che complexity of the drilling process, it was difficult to estimate the 
amount of heat dissipated and che distribution of heat along che borestem. 
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Measurements made within minutes after the probes were Inserted showed that 
the probes were approximately 40 to 50 "K warmer than the surrounding un- 
disturbed lunar soil. Because of the very low thermal conductivity of the lunar 
soil, this heat was dissipated very slowly. Even after 45 days, some of 
the thermometers indicated continued cooling. 

In general, the rate of equilibration of each probv sensor to the sur- 
rounding lunar medium depends on several different parameters, including thermal 
and geometrical properties of the probe/borestem system, thenr.il properties 
of the adjacent lunar medium, the thermal contact between the borestem and the 
adjacent zone of disrupted lunar material, and the total heat energy excess 
of the system upon completion of drilling. For long times ( t>20 hr) after the 
probes were inserted, however, the probe sensor temperatures above equilibrium 
depended essentially on two quantities, the total initial energy excess per unit 
length of the system AE and the thermal conductivity of the adjacent regolith 


T (t) - T - AE/4trk t (C-l) 

p « O 

where T p (t) is a probe sensor temperature at time t since emplacement and T m 
is the equilibrium sensor temperature. 

The form of equation (C-l) has been verified by finite difference models 
and by analytical solution appropriate to simplified distribution of AE 
between the probe and borestem. The analytical solution, which cakes the form 
of equation (C-l) for long times, is described in detail in the appendix of 
Langseth et al. (1972). 

Estimates of the regolith conductivity k^ made from the long-time cooldown 
data using equation (C-l) once a value for AE was assumed. Because of the 
uncertainties of the contribution of the drilling process to the total excess 
heat energy, two different assumptions were made in regard to the AE term. 
First, it was assumed that drill heating effects were negligible and that all 
excess heat energy initially resided in the probe and borestem. The probe and 
borestem were assumed to be at the same initial temperature T q after drilling, 
as measured by the probe sensors. For this assumption, the initial energy 
excess can then be expressed simply as 

s 6 )-(t 0 - tj 


(0-2) 


J 


where S p and ere Che heac capacities per unic length of the probe and borestem, 
respectively. The resultant problem has an analytical solution that was used 
to determine regolith conductivity estimates from the cooldown data, lhese 
assumptions constituted the minimum Initial energy case. 

For a maximum initial energy estimate, it was assumed that a 2.2 am thickness 
of disrupted lunar material had received enough h :at during the drilling pro- 
cess to attain an initial temperature equal to that of the probe and borestem. 

The 2.2 m contact zone thickness was chosen to correlate with the difference 
in drill bit radius and borestem radius. The subsequent cooldown problem was 
Chen solved, assuming radial heat dissipation, by using a finite difference 
model. The curve matching procedure was performed for times >20 hr when all 
parameter effects, with the exception of k m> were negligible; that is, when 
equation (C-l) was valid. 

The cooling histories for all subsurface sensors chat were not affected 
by diurnal variations were analyzed to estimate the thermal conductivity of 
the surrounding regolith for the two limiting estimates of initial temperature. 
Results have been tabulated in Langseth et al., 1972, 1973 and are presented 
with the other conductivity measurements in Figure V-C-16 of this report. 

b. Beater-activated conductivity measurements 

Following substantial equilibration of the probe-bores tem systems to Che 
surrounding regolith thermal regimes, active conductivity measurements were 
made using the thousand Ohm heaters that surrounded each of the gradient sensor 
housings. The experiments could have been operated in either of two modes 
by energizing the heaters at 0.002 or 0.5 W, depending on the conductivity of 
the surrounding regolith. Because of the low conductivities Inferred from 
tie cooldown analysis only the lower power was used at both heat-flow sites. 

The discussion is therefore confined to this mode. 

After each heater was energized, the temperature rise at the gradient sensor 
enclosed by the heater depended in a complex way on the thermal properties 
of nwarby probe components, borestem, and lunar material, as well as the contact 
zones between these elements. A detailed finite difference model of the 
conductivity experiment at each heater location was used to interpret the 
temperature history of the gradient sensor in terms of the conductivity of 
material external to the borestem. Numerical model computations, laboratory 
experiments, and lunar experiments indicated that, for times >20 hr, the 
temperature rise AT(t) was well defined by the simple relationship 
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AT(t) « c^in (c) + c 2 (C-3) 

where c^ and c 2 are constants . The form of equation (C-3) is the same as that 
for a heated infinite cylinder in an infinite homogeneous mediuu after many 
time constants. Indicating that axial gradients at the heater location are 
quite small and that the principal heat dissipation at long times takes place 
in the adjacent regollth. 

For an Infinite cylindrical source 

4 

c. - Q/4irk (C-4) 

i in 

where Q is the power per unit length in W/cm. Thus, c^ depends solely on 
the heater power per unit length and on conductivity. The constant c^ can 
be determined- easily because it is the slope of thv* temperature rise curve when 
plotted versus In t; therefore, cylindrical sources are often used as a practical 
technique for measuring conductivity. 

Conductivity was determined from lunar experiments by comparing reserved 
slopes on a logarithmic time scale with values of c^ calculated wit;, the finite 
difference models. Parametric studies, in which certain thermal properties 
were varied singly in the numerical model, showed that for times >20 hr, c- 
was very nearly Insensitive to changes in pc of the surrounding medium, 
changes in borestem conductance, and changes in the thermal links between the 
probe and borestem and the borestem and lunar medium. However, c^ was sensitive 
to changes in conductance in the probe body, which could alter the flow of 
heat from the heater axially along the probe. Assumptions of thermal properties 
in the numerical models that influence axial heat transfer along the probe 
were probably the largest source of error in the heater activated conductivity 
determinations. 

Active conductivity measurements were carried out at ail subsurface gradient 
sensors free of large diurnal variations within the first sixty days following 
deployment of each heat-flow station. Heaters were sequentially energized for 
periods of 36 hours or more and local conductivity estimates were made by 
matching long-time slopes of the gradient sensor temperature rise with pre- 
dictions of the finite difference models. Small transients due to incomplete 
Initial cooldown equilibration or slight diurnal effects were accounted for by 
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Figure V-C-l 




extrapolation of treads observed immediately prior to heater turn-on. The 
absolute level of observed temperature rise at a given heater location was 
utilized to estimate the contact resistance between the borestem and adjacent 
regolith by comparison with the finite difference models. Examples of the 
curve matching exercise and subsequent thermal property results can be found 
in Langseth et air, 1972, 1973. The heater activated conductivity results are 
included in Figure V-C-16 of this report. 

2_. Long-Time Temperature History Analysis 

a. The Probe-Borestem-Regolith Numerical Model 

It became increasingly apparent with the return of temperature data from 
the Apollo 15 heat-flow probes that the exchange of heat within the probe- 
borestem regolith system would require careful analysis before much of the 
data could be utilized. Because of the shallow emplacement of both heat- 
flow probes only the lower section of probe 1 was deployed deep enough to be 
free of measurable diurnal variations. The presence of significant higher 
frequencies and the absence of a discernible phase lag, which characterized 
the probe diurnal variations, implied the importance of radiative transfer 
along the probe-borestem axis (significant annual temperature variations 
were observed at all sensors. In addition, thermocouple measurements within 
the borestem section protruding above the surface indicated that the diurnal 
variation at the top of the borestem differed significantly from the lunar 
surface diurnal response in amplitude, frequency composition and average 
temperature over one lunation. Thus, to interpret the time varying data i r 
terms of regolith thermal properties and the steady-state regolith gradient, 
detailed numerical models of both the probe 1 and probe 2 in situ geomecrles 
were developed. Including both conductive and radiative heat exchange. A 
schematic of the probe 2 model is shown in Figure V-C-l. Vertical dimensions 
are to scale. Radial dimensions are indicated along the bottom boundary. 
Radial symmetry is assumed. Within the regolith thermal properties are 
modeled as both temperature and depth dependent. Lateral variation of 
properties is limited to the contact zone adjacent to the borestem. A fixed 
finite value of conductivity is used to evaluate conductive exchange verti- 
cally along the borestem walls. Distinct conductivities are assigned to 
each of the probe nodes for evaluation of axial conductive exchange. The 
variation of conductivity along the probe depends on the number of electrical 
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connections within a particular probe zone. The zones corresponding to the 
large platinum resistance thermometer location^ at either end of each probe 
section, are nearly ten times as conductive as other probe zones. Heat 
capacitance variation is accounted for in a similar way. 

Blackbody radiative exchange takes place at the surfaces of all borestem, 
probe and cable zones. View factors between all radiating zones are calculated 
based on the formulas found in Hottel and Sarofim ( 19*67) for concentric 
cylinders. The surface boundary condition balances the outgoing flux, solar 
insolation and the temperature gradient at the surface. A similar flux balance 
equation governs the heat balance of the exposed borestem zones. At the 
bottom boundary, 15 cm below the deepest drill bit penetration, a constant 
heat flux is assumed. At the far radial boundary, lateral heat flux is set 
to zero. 

The computation scheme utilizes forward difference equations which essen- 
tially assume a steady flux across each node boundary over the time step of 
about 100 sec. The general formulation of the difference equation is outlined 
in Appendix A. Temperature dependent thermal properties are recomputed after 
each time step. Initial conditions were established by running one dimensional 
regolith models with the same vertical zone sizes, chermal properties and solar 

insolation function as the complete model until a periodic steady state solution 
« 

was achieved. In this way, predeployment regolith temperatures were established. 
The regolith thermal property profile used was based on the preliminary in situ 
conductivity measurements and near surface (0-15 cm) conductivity deductions 
inferred from the measured temperature histories (see Langseth et al., 1972; 
Keihm et al., 1973). Temperature dependence of regolith conductivity was 
derived from the measured mean temperature difference between the surface and 
the shallowest probe sensor at 35 cm. The density profile *ns based primary 
on the core tube measurements of Carrier et al. (1973). S -ic heat variation 

with temperature was modelled as a polynomial fit to the data of Robie et al. 
(1970). It is to be emphasized that the principal objective of the two dimen- 
sional models was an attempt to understand the complex thermal interaction of 
the probe-borestem system and surrounding regolith and subsequently to interpret 
sensor temperature variations with time and depth in terms of the regolith 
variations and thermal properties. For this reason, uncertainties in the 
choice of regolith thermal properties were of secondary importance. The 
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primary question was to what extent did the probe-borestem thermal linkages 
modify the temperature regime which it was designed to measure. 

the first models runs were carried through three lunations with no annual 
component present in the solar insolation function. The purpose of these 
initial runs were twofold: first, to test the model's capability of repro- 

ducing the diurnal characteristics of the real time data, and, second, to 
examine the feasibility of interpreting the diurnal variations and time averages 
of the actual data in terms of regolith thermal properties and gradients. 
Comparisons with the lunar data containing significant diurnal variations is 
shown in Figure V-C-2. Note the difference in shape between the probe 1 and 
probe 2 lunar data. The modelling results simulate the high frequency content 
of the probe 2 data but do not compare well with the smoother, more surface-like 
response of the probe 1 data (Differences in amplitude between the modeling 
results and the lunar data indicate that our original choice of regolith 
thermal properties based on active conductivity experiments required 
adjustment. However, the matching of diurnal amplitudes observed at various 
depths along the probe by varying the regolith conductivity is neither a 
straightforward nor unambiguous process. In fact, it will be shown in the 
next section that the regolith conductivity must be lowered to raise the model's 
diurnal amplitudes to the level observed at depth. This result agrees with the' 
later analysis of the long-term data.) The probe 1 and 2 geometries differed 
only in depth of probe emplacement and in the presence of a disk shaped radiation 
shield made of reflective material above probe 1 within the borehole (Figure V- 
C-3) . Both of these differences have been accounted for in the probe 1 and 2 
models, yet the thermal model of the probe 1 response is essentially similar to 
that at probe 2. Various model parameters were adjusted in an attempt to re- 
produce the shape of the probe 1 diurnal data. These included varying the 
reflectivity of the radiation shield, the borestem-regolith contact resistance 
and the borestem and regolith conductivity. Yet, the shape of the modelled 
diurnal response remained essentially unchanged - a steady rise in temperature 
from sunrise throughout the lunar day with a sharp transition to cooling im- 
mediately after sunset. The lack of a significant morning phase lag and the 
absence of afternoon cooling Indicate that the probe sensors are responding 
primarily to diurnal pulse propagation down the probe-borestem hole rather 
than the diurnal response of the adjacent regolith. In large measure this 
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Comparison of diurnal phase lags and amplitudes for undisturbed regolith (FRB) 
and probe sensors using two different regolith conductivity profiles in the 
Apollo 15 probe numerical model. 


Figure V-C-4 
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heat transfer is by radiation along the tube. Unlike the lunar surface, the 
protruding borestem maintains a nearly constant high temperature during the 
lunar day due to the orthogonality of the borestem-lunar surface orientation. 
Daytime flux into the protruding borestem derives both from direct solar 
heating and radiation from the lunar surface. The solar contribution maxi- 
mizes near sunrise and sunset while the heat emitted by the surface peaks 
near noon. The summing of these effects results in the approximate constancy 
of daytime protruding borestem temperatures, which serve as a hot source for 
the axial pumping of heat downward during the day. Radiative transfer explains 
the nearly immediate response of the probe sensors at sunrise. Both radiation 
and borestem wall conduction subsequently deliver heat to the probe sensors from 
the protruding borestem 'reservoir' throughout the lunar day; at sunset, the 
borestem cools almost immeiiately as the opposite process takes place. 

The question remains as to why the probe 1 sensor's diurnal response more 
closely resembles the regolith diurnal variation than the probe 2 data or the 
modelling results for either probe. One feasible explanation is that the re- 
flective plate which caps each borestem appears from photographs to be some- 
what tilted up at the probe 1 site, possibly resulting in radiation from the 
lunar surface being reflected directly Into the probe 1 borehole. This could 
explain the early afternoon peaking of probe 1 sensor temperatures since they 
would be more directly coupled with actual surface temperature variations than 
the probe 2 sensors or either model. 

b. Probe-Regolith Coupling; Diurnal Variations 

It is apparent from the data shown in Figure V-C-2 that the axial thermal 
linkages of the probe-borestem system contribute to the character of the probe 
sensors' diurnal variations. To assess the effects of regolith thermal pro- 
perties on the probe response, models were run with different regolith thermal 
properties for three lunations each. Although this is not normally enough time 
for a periodic steady-state solution to be obtained, it is sufficient to esti- 
mate the effects of regolith properties on probe diurnal response. Results are 
shown in Figure V-C-4. For each of the two models shown, the diurnal tempe- 
rature variations of all probe and far radial boundary (undisturbed) regolith 
nodes were Fourier analyzed, using the second and third lunations of the three 
lunation computer runs, to yield the amplitude and relative phase lag of the 
diurnal components. Models Cl and C2 differed only in regolith thermal 
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Comparison of mean temperature, annual phase lag and annual amplitude between 
undisturbed regolith (FRB) and probe as computed from 2.2 year run of the 
Apollo 15 probe 2 numerical model. 


Figure V-C-5 



76 


conductivity below 20 cm: GO ? 1*5 (x) for x > 20 (Figure V-C-4a) . 

Note in Figure V-C-4b the suppression of the diurnal phase lag with depth 
as sensed by the probe. This results directly from the radiative coupling 
of the subsurface probe and the above-surface borestem section. The in- 
creased conductivity model does suppress the phase lag even further; how- 
ever, the resolution of average regollth conductivity based on the probe 
phase lag data is much less than that implied by the undisturbed regollth 
zon6s • 

Figure V-C-4c illustrates most clearly the importance of the probe- 
borestem thermal linkages in determining the probe's diurnal response. At 
the lunation frequency, the regollth acts primarily as a medium into which 
the diurnal wave disperses by conduction effectively 'lowering the probe 
diurnal response for higher regollth conductivity. Thus, we see that the 
diurnal amplitude response of regollth and probe are essentially negatively 
coupled in that a higher regollth conductivity results in larger diurnal 
variations within the regollth (lower attenuation) but reduced variations 
within the probe (higher attenuation) . This effect dominates at all depths 
below 20 cm, but quantitatively depends In a complex manner on the amplitude 
and phase correlation between the probe and regollth variations at a given 
depth. For these reasons, the utilization of probe diurnal variations as 
an accurate means to estimate regollth thermal properties had to be abandoned. 

c. Probe-Regolith Coupling: Annual Variations in Mean Temperature Profile. 

Previous models, both analytic and numerical, had been utilized to examine 
the effect of the relatively high conducting probe-bores tem system upon the 
local regollth steady-state gradient. It was found that for depths at which 
diurnal variations were negligible, the gradient as measured by the probe 
differed from the undisturbed local regollth value by less than 51. However, 
the shallow emplacement at the Apollo 15 site did result in significant 
diurnal variations at four of the working probe sensors and the question of 
how reliably the time average probe temperature profiles reflected the regollth 
gradient required investigation. Of prime importance was the effect of the 
mean temperature difference between the lunar surface and the protruding 
borestem section and the resultant equilibrium thermal regime in the sub- 
surface. 

In addition to the shallow sensor diurnal components, annual temperature 
variations were detected in the long-term data at all but the deepest probe 1 



sensors. Once again, correlation between regollth annual variations and 
thermal properties demanded investigation. To examine the correlation of 
both the shallow sensor mean temperature gradient and the measured annual 
variations with the comparable temperature components in the undisturbed 
regollth, the probe 2 numerical model was run with the nominal thermal 
property profile for a simulated period of 2.2 years. The run of over 27 
lunations was required to achieve a satisfactory periodic steady-state 
solution for the probe-borestem-regolith thermal regime and to allow an 
accurate Fourier analysis to obtain the annual components. Diurnal and 
annual components were identified using the same finite Fourier transform 
technique as was employed for the actual lunar data. Both probe nodes and 
far radial boundary nodes (undisturbed by the probe presence) were analyzed. 
Probe and regollth comparison of mean temperatures, phase lag, and annual 
amplitude profiles are snown in Figure V-C-S. The large offset in mean 
temperature (Figure V-C-5a) at depths less than 50 cm results directly from 
the mean temperature difference between the protruding borestem and the 
lunar surface. Only at depths below significant diurnal variations ( < 0.05*K) 
does the regollth temperature regime totally constrain the probe gradient. 
Sensors in the region of significant diurnal variation can be corrected 
using the probe models. These corrections have been applied to both the 
probe 1 and probe 2 shallow sensor data based on the results shown in Figure 
V-C-5a. 

Figures V-C-5 b and c illustrate the excellent probe-regolith coupling 
for the long period annual component upon which the subsequent revision of 
regollth thermal conductivity and heat flow was based. The nearly constant 
offset in phase between the probe and undisturbed regollth is established 
close to the surface and does not affect the subsurface conductivity deductions 
which are constrained only by the rate of phase change with depth. The probe 
and regollth annual amplitude components can be considered virtually identical 
at all depths if allowances are made for the uncertainties present in the 
Fourier transform analysis. These results, when compared with the near 
absence of , or negative probe-regolith coupling, for diurnal variations 
shown in Figure V-C-4a are indicative of the frequency dependence of the 
probe-regolith coupling for periodic temperature variations. At frequencies 
as high as the diurnal, axial heat transfer by radiation and probe-borestem 
conduction is a primary factor in controlling the subsurface probe response. 
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The high frequency of Che lunacion cycle aicigaces Che effecc of cemperacure 
varlaciona wichin Che adjacenc regollch which are eheraally coupled co Che 
probe-borescea system only by Che highly insulacing lunar soil. Annual 
varlaciona, on che ocher hand, are of such low frequency chac ac any given 
cloe during Che year, a near equllibriua channel regime is escablished in 
which Che probe CeaperaCures are alaosc co rally conscrained by Che adjacenc 
regollch CeaperaCures . These aodelling resulcs provided che basis for a 
reevaluaelon of che regollch eheraal conduceiviey froa che Apollo 15 and 17 
annual wave daca. 
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F ormulation of che Flni te- uifferenc e S ch e me for_ the 
Probe- Bo res tem-Regolith Numerical Model 


The computational scheme utilized assumes that all thermal properties 
and node temperatures remain constant over one time step. The temperature 
change at any given node over the time step is Chen simply the total heat 
flux across all node boundary surfaces multiplied by the time step and 
divided by the heat capacity of the node. 


AT. 


At 

V o C 
i p i i 


Ni 

J-l F U 


(1A) 


where At is the time step. 

V l’ P i’ * r * tlie volume, density and heat capacity of node i, and 
Fjj is the flux into node 1 across its jch boundary, summed over each of 
its N^ boundaries. 

Zn the cylindrical coordinate system, for an Interior node with no 
free surface boundaries, Equation 1A can be expressed as: 


“l.l ■ 0?^ 


T 1-1,J ~ T l,1 T l+l,1 ~ T l,1 T 1,J-1 ~ T l.j + 
<“>1-1.1 <“>1, i (R,) 1.J-1 


T - T 

1.3 


(RV) 


i.3 J 


(2a) 


where (RR) . is the thermal resistance in the radial direction between 

J 

node centers i,j and i + l,j and (RV) . . is the thermal resistance in the 

J 

vertical direction between node centers i,j and £, j +1. 

The thermal resistances can be calculated from the node geometry under 
che assumption of constant flu* over the duration of che time step. 









A 


In the vertical direction. 
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Similarly, 


“i-i j * DZ i ' CT i.i> 

H(ro Z - ri Z ) 


DZ 1 / (2 *i.1> 1 DZ 1+1 < 2K 1,H-1> 

II (to| - ri 2 ) 


At the lunar surface, the flux boundary condition is: 


* -71 

s <Z z»o 


eoT - (1 - A) • I (t) 
s 


where 


E is the infrared surface emittance 


'J is the Stefan-Boltzmann constant 
A is the solar alt °.do : and, 

I(t) is the time dependent solar flux. 

Equation (.*^1 is expressed in finite-difference form for each of the 
lunar surface ? aes by using the approximation: 


SSSSSSS 


5T , 

~lz ] z-c 


2 (T - T ) DZ ( 2(T2 - T ^ - 2 < Tl ~ T Al ) 
_ tT i V _ u *i l P Zrl-DZ 1 DZi ’ 

dz l : dz 1 + 2 dz 2 


(10A) 


With the ubstitution of (10A) into (9A) , the non-linear surface boundary 
condition can be solved iteratively for T g using che Newtan-Raphson method, 
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For the borestem, probe and cable nodes which possess free boundary 
surfaces, radiative exchange is calculated with all other line of sight 
nods surfaces by assuming unit emittance and utilizing computed direct 
exchange areas based on the concentric cylinder formulae given in Hottel 
and Sarofim (19 67) . The net radiative flux into all such node surfaces 
is then calculated (utilizing the temperature distribution at a given time 
step) and substituted in Equation 1A. 


* 




KJ go 



d. Filtering and Smoothing Techniques: Id entification of the Annua l 

C omponents 

Revised diffusivity values resulting from analyses of the two and 
three and a half year subsurface temperature histories at the Apollo 17 
and IS sites have been published (Langseth et al., 1976). In this section 
we describe in more detail the filtering, smoothing and curve fitting 
techniques utilized in the long-term analyses. 

The primary tool usad for component separation (diurnal, annual, tran- 
sient and noise) was .. finite Fourier transform computer program 

T. e program required as input 2° equally spaced (in 
time) data points with n any integer. The temperature history data were 
converted to this form by linear interpolation between adjacent data. At 
the Apollo 15 gradient sensors, 4096 (n * 12) data points were utilized over 
approximately a three year window. For the ring bridge sensors, 512 (n ■ 9) 
data points were input into the transform program. In all cases, the output 
consisted of 2° -1 pairs of the cosine and sine components. 

At the Apollo 15 site, a two step transform process was used. First, a 
data window of exactly 38 mean synodic periods was transformed, utilizing the 
available data up to December 31, 1974 (Step 1). This transform window 
encompassed over three years of data but excluded the first few lunations 
following deployment to minimize the initial condition effects. Represen- 
tative spectra are shown in Figure V-C-6. (see Figure V-C-8 
for the time series data). At the shallowest sensors (TR11A, TG22A) , signi- 
ficant diurnal harmonics can be clearly seen at frequencies six to eight 
times the lunation fundamental (n * 38) . Significant noarandom power is 
also present at near diurnal frequencies due to the variation in the synodic 
period throughout the year. At a deeper sensor such as TG11B, diurnal com- 
ponents are still significant, although strongly att.nuated relative to the 
shallower depths. At the deepest sensors (such as TG22B) the barely dis- 
cernible diurnal component probably results from a day-night offset in the 
electronic's calibration. 

The large low frequency components are primarily the signature of the 
long-term transient rise present at all sensors. At the n • 3 harmonic, 
significant power is produced by the annual component, especially at the 
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Figure V-C-6 




shallower sensors (For a 38 lunation transform window, the precise annual 
frequency would be at n * 3,07). Since the purpose of the long-term ana- 
lysis was tho precise Identification of the annual components and pure 
transient response, the spectra such as shown in Figure V-C-6 were smoothed 
to eliminate the diurnal components and minimize the high frequency noise. 

For this purpose, the high frequency components (n > 12) , excepting the 
diurnal and surrounding frequencies, were fitted with analytic expressions 
consistent with the general trends of the sine and cosine high frequency 
components. For the sine components an exponential decay function of the 
form 

-Ottl 

S Cn) » S Q ® 

was utilized. For the cosine component, note that the high frequency noise 
has a definite positive bias. This property is a characteristic of the dis- 
crete Fourier transform utilized. The high frequency cosine noise averages 
to zero asymptotically with the number of data points used in the transform. 
For this reason, a linear fit to the n > 12, harmonics was employed. An 
example of the spectrum filtering fit, utilized for TR11A, is shown in 
Figure V-C-7. Note the absence of the diurnal and near diurnal components 
which were not used in the fits. The cosine fit (a small positive constant) 
shows that the transient power is contained primarily in the low cosine 
frequencies. As pointed out earlier, the non-zero high frequency noise 
average is an artifact of the Fourier transform computation technique. The 
experimental decay fit to the high frequency sine components is compatible 
with the theoretical spectrum of an idealized transient response (discussed 
in the next section) . Also shown in Figure V-C-7 is the inversion of the 
filtered transform showing the time series data smoothed and free of the 
diurnal components. These smoothed data were then transformed using a window 
exactly three years long (Step 2). This spectrum is shown on the right. The 
component at n • 3 contains power both from the annual variation and the long- 
term transient. Separating the two contributions was the primary goal of the 
long-term analysis. At the shallower sensors, such as TF11A (45 cm), the 
annual and transient contribution a.,e of comparable magnitude and a good 
firsc approximation of the annual component can be obtained by interpolating 
for the n • 3 transient contribution from the surrounding spectra which re- 
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Figure V-C-8 


Real time data and two year frequency spectra from Apollo 17 TG11A and DT21 data 


38 


present transient power only (see the three year transform of Figure 
V-C-7) . At the deeper sensors .however, the transient contribution at the 
annual frequency is many times larger than the annual component itself, and 
interpolative techniques cannot always be used. For example, at the Apollo 
17 site, the shallowest probe sensors were at a depth of 1.3 m. Discernible 
annual variations were present in the real time data, but did not stand out 
clearly from the transient background in the frequency spectrum. This is 
illustrated in Figure ^-0-8 which shows the time and frequency domain data 
at the Apollo 17 absolute temperature sensor TG11A and the more accurate 
temperature difference data at DT21. Note that only two years of data were 
available for frequency analyses. Thus, the fundamental annual component 
lies at the n * 2 harmonic of the transform. The primary difficulcy is 
caused by the large amplitude at the n ■ 1 frequency due to the aperiodic 
transient. The shape of the transient background spectrum is too irregular 
near the n « 2 annual frequency to allow a reliable interpolation to elimi- 
nate the transient contribution. (The modulation of the spectra seen at higher 
frequencies is a resuJ t of the short term conductivity experiments characterized 
by the splke3 within the first lunation of the time history data. The diurnal 
components seen in the absolute temperature data at TG11A result from a day- 
night bias in the electronics calibration and do not represent a real pene- 
tration of the diurnal wave to a 130 cm depth.) The most precise method for 
separating the annual components trom the transient background is an examination 
by eye of the filtered real time data. First, any diurnal variations, if 
present, and the high frequency noise were filtered out as described previously. 
The resulting real time data containing only transient and annual components 
were then replotted on as large a scale as the residual noise level permitted. 

An initial estimate of the sine and cosine annual components was then made by 
lnterpoletlon and the temperature history data replotted with the assumed 
annual components subtracted out. The resultant plot was then carefully 
examined by eye for any discernible annual periodicity remaining. This 
process was repeated for various values of assumed annual components until no 
annual fluctuation could be detected in the replotted data. The resolution of 
the trial and error method at representative sensors of varying precision is 
illustrated by the test models shown in Figures V-C-9 through V-C-12. In all 
cases, the vertical scale refers to the filtered data with no annual component 
removed. If feasible, initial estimates <- f Che sine and cosine components 








Smoothed real time data with various annual components removed 
A15 TG12B 





Smoothed real time data with various annual components removed 
A17 TC23 




are made by interpolating from the surrounding pure transient components of 
the frequency spectrum. For example, at the shallowest Apollo 15 probe sensor 
(TG11A, Figure V-C-9) the three year spectrum exhibits a distinctive cosine 
component at the n • 3 annual frequency, most of which is due to the annual 
periodicity. An interpolation from the su rounding frequencies (n - 1,2) 
yields an Initial estimate of Ac » - 0.272; As « * 0.017. Further refinement 
ia obtained by replotting the real time data with a feasible range of annual 
components subtracted. On the scale shown, annual fluctuations can bt dis- 
cerned by eye in all cases except for Ac ■ - 0.260; As • 0.005 and Ac ■ -0.275; 
As - - 0.005. (Irregularities of non-annual frequency result from the filtering 
techniques described earlier used to eliminate the diurnal components and 
minimize the high requency noi_»- ) . Further resolution of the annual com- 
ponents by plotting on expanded scale resulted in final estimates of 
Ac - - 0.265 (+ 0.010) and As - - 0.005 (£ 0.010) for TG11A. 

Limitations of the 'eyeba.’ ^ing' technique are illustrated in Figure V-C-10. 
Initial spectrum interpolation ascisuites indicate that the annual component is 
very near zero (<0.0i0) at the depth of TG12B (138 cm). (Note the difference 
in scales for the sine and' cosine spectra). Inspection of the real time data 
confirms the spectrum indications. The subtraction of any sine component 
larger than |0.0l| results in an easily discernible annual variation while no 
such variation can be seen la the time history data assuming zero annual com- 
ponent. F.v this reason, further resolutions obtained by eyeballing on ex- 
panded scale resulted in the establishment of an upper bound for the annual 
component of j Ac | < 0.003; { As | < 0.005. Note that this resolution Is an 

order of magnitude below the noise level of the absolute temperature data. 

At the Apollo 17 site, the increased emplacement depths severely limited 
the usefulness of the absolute temperature data for the determination of 
annual components. Only the uppermost probe sensors (130 cm exhibited annual 
variations detectable within the absolute temoeratt.re noise (+ 0.05”K). To 
obtain annual components at other depths, both the more accurate temperature 
difference (DT) data and the shallowest (15 cm) thermocouples were examined. 
Analysis for the Apollo 17 0T12 gradient sensor is shown in Figure V-C-ll. 

At the time of analysis, 2.7 years of Apollo 17 data were available. Thus, 
the .'ourier analysis was done over a two year window and the real time data 
were Inspected with various annual components removed ove*’ the full 2.7 years. 
From the spectra, it can be seeu that a cosine component of approximately 
-0.0035 is required. Because of the shape of the sine spectrum, the sine 
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TABLE V-C-l 


Annual Components at A-15 and A-17 Subsurface Sensors 



Depth 

A 

s 

A c 

/a 2 .a 2 

s c 

Apollo 15, Probe 1 

35 

-0.005 

-0.290 

0.290 


45 

-0.064 

-0.170 

0.182 


73 

-0.075 

-0.010 

0.076 


33 

-0.057 

-0.008 

0-. 053 


*1 

-0.037 

+0.009 

0.038 


101 

-0.030 

+0.012 

0.032 


129 

<0.01 

<0.01 

— 


139 

<0.01 

<0.01 

— 

Apollo 15, Probe 2 

49 

-0.265 • 

-0.150 

0.305 


59 

-0.140 

-0.170 

0.220 


87 

+0.015 

-0.175 

0.076 


97 

+0.025 

-0.050 

0.056 

Apollo 17, Probe 1 

15 

0.63 

0.53 

0.82 


130 

-0.020 

0.006 

0.021 


186 

-0.0012 

-0.C.36 

C.0038 

Apollo 17, Probe 2 

15 

0.64 

0.15 

0.66 


131 

-0.014 

-0.008 

0.016 


187 

-0.0009 

-0.0020 

0.0022 
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FV >1 
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component cannot be directly estimated. The real time test plots show the 
resolution available for sine components ranging from *0.003 to + 0.001. 

Even on the limited scale of Figure V-C-ll, it can be seen that a value of 
As 0.001 is required to eliminate a discernible annual variation. Further 
eyeballing anal vs is on expanded scale yields Ac * - 0.0036 (+ 0.0002) and 
As *- 0.0012 (+ 0.0002). Again, the resolution obtained is an order of 
magnitude increase over the accuracy of the data. The temperature difference 
annual components were assumed to be representative of the annual variation 
at the top gradient sensor (185 cm) . This is equivalent to assuming zero 
annual 'iriation at the lower gradient sensor, a justifiable assumption in 
view the low thermal diffuslvlty of the lunar regolith. 

The Apollo 17 thermocouple data TC23 (15 cm), shown in Figure V-C-12, are 
accurate only to + 0.5°K. However, at the 15 cm depth, the annual component 
is large enough to permit a good first esclmate from the two year spectra. 

The time series data shown have been filtered to remove a diurnal component 
and to minimize the thermocouple noise. Cosine component resolution is illus- 
trated. Further trial and error analysis on expanded scale yielded Ac ■ 0.18 
(+ 0.03) and As - 0.64 (+ 0.04)’K. 

Annual components at all analyzed Apollo 15 and 17 sensors are presented in 
Table V-C-l • Resultant bulk regolith conductivity estimates based on the 
annual amplitude depth attenuation are discussed in Langseth et al. (1976) and 
at the end of this section 

e. Annual Wave Diffuslvlty Deductions 

For a medium of uniform thermal properties the amplitude of a harmonic 
temperature variation at the surface attenuates with depth with a decay cons- 
tant 8 » A)/2jc, where u is the fundamental frequency and < the thermal dif fusi- 

vlty (Carslaw and Jaeger, 1959). For the lunar situation, departures from 
homogeneity occur primarily at shallow depths (less than 10 cm). At these 
depths, non-linearities produced by the strong temperature dependence of 
thermal properties must also be considered. However, all probe sensors over 
which temperature gradient measurements were made lie at least 35 cm below the 
surface. At these depths, diurnal variations are less than 6°K and the tem- 
perature dependence of thermal properties need not be considered. In addition, 
if one measures the attenuation from one sensor to the next, the resultant 
diffuslvlty estimate will represent, to first order, an effective bulk dlffusi- 
vity between the two sensor depths, regardless of the properties of near-sur- 
face layers. 


N» 




Figure V-C-13 






Figure V-C-13 shows the attenuation of annual temperature variations with 

depth at the Apollo 15 and 17 sites. Within the accuracy of the amplitude 

determinations, linear fits to the exponential attenuations were made, with the 

best filling attenuation curves shown in the figure. The range of acceptably 

-4 

fitting curves corresponds to a range of diffusivities of + Q.08 x 10 . 

The sparseness of the data from the Apollo 17 site is due to the relatively 
deep probe emplacement depths. Annual wave amplitudes could be resolved only 
for the absolute temperature histories of the top probe sensors. Amplitudes 
at depths of nearly 2 m were derived from the more accurate temperature dif- 
ference data with the assumption that all annual variations were due to varia- 
tions in the upper sensor temperature. The shallowest Apollo 17 data points, 
at 15 cm, were calculated from thermocouple temperature histories. The narrow 
range of linear fits again supports the assumption of regolith thermal property 
homogeneity at depths below 10 cm. 


f. Theoretical Fits to the Long-Term Transients 

Following subtraction of our best estimates of annual components from 
the time series data, the remaining transient histories, propagating down- 
ward from the disturbed surface region, were analyzed to provide an additional 
estimate of the bulk regolith thermal properties. 

In actuality, the region of perturbed mean surface temperature is made up of 
a patchwork distribution of footprint indentations, scuffed surface fines and 
some undisturbed surface elements within a few meters of each probe. The 
relative contribution of each disturbed surface element to the subsurface 
transient behavior will depend on the solid angle subtended at a given probe 
depth by that surface element: 


dT - 


B 

2rr 


J£77 

/ Vk t 



2 + A 


4 K 


+ erfc (. 


L2. , 2 

/Z + r 

2/kT t“ 



( 1 ) 


(Lachenbruch, 1957) 


In the above equation dT is the temperature change at probe depth Z and time 
t due to a temperature change B over a surface element subtending solid angle 
dfl at depth Z and located a lateral distance r from the probe axis. The 
regolith thermal dlffuslvity, k , is assumed uniform. Erfc is the complimentary 



error function. Equation (JL) indicates that at the shallower probe depths, 
transient disturbance will be weighted most heavily by the disturbed surface 
elements nearest the probe axis. Ac greater depths, the probe axis proximity 
of disturbed surface elements will be less of a factor in their relative 
effect on the subsurface transient response. 

Because the actual distribution of disturbed surface elements is unknown, 
we have fitted the subsurface transient histories to an idealized transient 
function appropriate to a uniform disk shaped surface disturbance of tem- 
perature magnitude B and radius R centered at the probe axis. For this model 
integration of Equation (1) yields 

jrr % 

AT (z , t) - B l erfc C - ) - - Z r , erfc (— LL sar- )] (2) 

2 / / Z 2 + R z 2/ <t 

where AT is the total transient effect at time t and probe depth Z due to 
the uniform disturbance of magnitude B and radial extent R. It is to be noted 
that the fitted diffuslvlty parameter, < , represents a bulk or average value 
between the surface and a given probe sensor depth. At both heat flow sites, 
however, the linearity of the mean temperature profiles and exponential annual 
amplitude attenuations Indicate nearly uniform bulk thermal properties at probe 
depths. Only the upper surface layer (2-10 cm) varies significantly with depth. 
One would therefore anticipate that the transient diffuslvlty fits should 
approach the annual wave results for the deeper probe sensors. 

In the fitting of Equation (2) to the probe transient histories, time t * 0 

(onset of 'instantaneous' change in the local near surface temperature) was 

taken to occur at the end of the first lunar night following deployment. In 

actuality, the establishment of the perturbed regions? mean surface temperatures 

takes place over the course of the first lunation cycle following probe 

emplacement. Earlier t - 0 times were tested (for example, immediately 

following the early morning deployment and at the first sunset) for the Apollo 

IS transient fits to assess the effect of the zero time choice on resultant 

diffuslvlty deductions. In no case did the results yield dlffusivltles differing 

-4 2 

by more than 0.1 x 10 cm /sec. At the deeper sensors, the effects of zero 
time choice were even less. 

In our early transient fits (Langseth et al., 1976) the transients for a 
particular probe were required to yield a constant value of the parameters 3 
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, and R as well as a depth profile of average diffusivlty, <(z) • This was 
based on the assumption that the disturbed region approximation represented 
> by a 'best' choice (in the least squares sense) of B and R in Equation (2) 
should apply equally well at all sensor depths. This, however, need not be 
the case. One would expect that the region of astronaut disruption be con- 
centrated within about one matef of the probe during the drilling of the hole. 

At distances greater than one meter, disruption occurs less frequently with 
the result that the angle averaged mean temperature change will decrease with 
distance from the borehole site. Thus, due to the solid angle effect ex- 
pressed in the weighting Equation (2) , shallower sensor transients will be 
more strongly affected by the near hole surface disturbance, while at deeper 
sensors, the more 'extensive disturbed region will bear more nearly equal 
weight for the transient response. For the refined transient fits, therefore, 
we have allowed depth variation of the best fitting disturbed region parameters 
B and R within constraints which will be described further on. 

Representative examples of the refined transient fits are shown in 
Figures V-C-14 and 15. The theoretical transient shown are meant to illustrate 
the sensitivity. of the fits to the diffusivlty parameter alone. Each transient 

was fit in the least squares sense for the parameters B and R for diffusivltles 

-4 . . 

of 0.40, 0.45... 1.5 x 10 . Values of R were restricted to 100 s R 5 300 cm 

at 50 cm intervals. The solid curves indicate the best fit obtainable by 
varying all three parameters. The dashed and dotted curves are the best fits 
obtainable using < values + 0.2 x 10 from the optimum allowing further vari- 
ation in B and R. At the Apollo 15 site (Figure V-C-14) , residuals in the range 

0.005 - 0.025*K were obtainable by optimization of all three parameters. Varying 

-4 

< by + 0.2 x 10 from its optimum value, and then optimizing 3 and R increased 
the residuals by a factor of 2-3 in most cases. At the deepest sensor (139 cm) 
comparable fits could be obtained for < • 0.50 - 0.80 x 10 (RMS 5 0.015). 
However, the lower range of these satisfactory diffusivlty values CO. 50-60) can 
be eliminated by the criteria that the average disturbance magnitude, B, cannot 
exceed that determined for the shallowest sensor. As discussed earlier, the 
shallowest sensor's transient response is heavily weighted toward the near-hole 
region of maximum disturbance. We, therefore, impose the additional constraint 
that the effective B value appropriate to the deeper sensors (which are sensitive 
to a greater extent of disturbed surface elements) must be less than, or 
equal to, that determined for the shallowest sensor at a particular probe site. 
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TABLE V-C-2 Theoretical Transient Fits 


Sensor 

Depth 

(ca) 

' • 4 *2 

(xl0 4 cn sec) 

R 

(cm) 

B 

(*R) 

RMS 

CK) 

K-Range 

Apollo 

15 

(Probe 1) 



• 


TG 11A 

35 

0.75 

200 

3.73 

0.0221 

0.65-0.90 

TR 11A 

45 

0.85 

250 

3.12 

0.0161 

0.75-0.95 

TR 11B 

73 

0.85 

300 

2.25 

0.0063 

0.75-0.90 

TG 11B 

83 

0.85 

300 

2.20 

0.0140 

0.65-0.95 

TG 12A 

91 

0.80 

300 

2.20 

0.0173 

0.70-1.00 

TR 12A 

101 

0.90 

300 

2.00 

0.0112 

0.75-1.00 

TR 12B 

129 

0.80 

300 

- 1.97 

0.0052 

0.65-0.85 

TG 123 

139 

0.75 

300 

2.20 

0.0070 

0.60-0.85 

(Probe 2) 







TG 22A 

49 

0.20 

100 

1.34 

0.0503 

0.10-0.80 

TR 22A 

59 

0.25 

100 

1.50 

0.0398 

0.15-0.70 

TR 22B 

87 

0.45 

200 

1.58 

0.0129 

0.35-0.60 

TG 22B 

97 

0.50 

250 

1.56 

0.0185 

0.30-0.65 

Apollo 

17 

(Probe 1) 




• 

TG 11A 

130 

0.70 

300 

3.90 

0.0046 

0.55-0.75 

TR 11A 

140 

0.75 

300 

4.10 

0.0122 

0.60-0.95 

TR 11B 

167 

0.70 

300 

4.46 

0.0040 

.0.60-0.80 

TG 113 

177 

0.75 

300 

4.00 

0.0022 

0.65-0.85 

TG 12A 

185 

0.75 

300 

3.49 

0.0039 

0.65-0.95 

TR 12A 

195 

0.70 

300 

4.16 

0.0014 

0.60-0,75 

TR 12B 

223 

0.70 

300 

4.15 

0.0024 

0.65-0.90 

TG 12B 

233 

0.75 

300 

4.10 

0.0025 

0.70-0.90 

(Probe 2) 







TG 21A 

131 

0.55 

300 

2.10 

0.0101 

0.50-0.70 

TR 21A 

141 

0.60 

300 

1.99 

0.0070 

0.55-0.75 

TR 2 IB 

168 

0.60 

300 

1.96 

0.0042 

0.55-0.70 

TG 2 IB 

178 

0.55 

300 

2.11 

0.0042 

0.50-0.70 

TG 22A 

186 

0.55 

300 

1.74 

0.0024 

0.50-0.70 

TR 22A 

196 

0.60 

300 

1.96 

0.0032 

0.55-0.70 

TR 22B 

224 

0.60 

300 

2.10 

0.0026 

0.55-0.70 
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Thus, for example, at the Apollo 17 probe 1 site (figure V-C-15) the best 
fitting curve at the deepest sensor (233 cm) requires an unreasonably large 
average mean temperature disturbance, which is incompatible with the transient 
fits at shallower sensors. The restriction of the fit to acceptable values of 
B ( *4.0) still results in residuals well within the precision of the filtered 
real time data ( ^O.OIK). 

Refined transient dlffuslvlty deductions for the Apollo IS and 17 sub- 
surface sensors are presented in Table V-C-2 . The parameter values presented 

in columns 3, 4, and 5 produce the best transient fits (in a least squares 

-4 2 

sense) for a given sensor. A parameter resolution of 0.05 x 10 cm /sec in 
thermal diffuslvl y and 50 cm in R was obtained. The dlffuslvlty range 
presented in column 7 shows the range of < values which can produce fits (by 
further variation of B and R) with BUS residuals within a factor of two of 
the indicated RMS best fit. 

Values of the average surface mean temperature disturbance were restricted 
to values no less than 0.5B* and no greater than 1.2B*, where B* is the B 
value deduced from the uppermost sensor data at a given probe site. This 
criterion is based on the assumption that the shallowest sensor, where the 
measured transient is largest, best reflects the surface disturbance nearest 
the borestem axis where effects are the largest. For example, the fits at 
the Apollo 15 probe 1 site require significantly smaller values of B at the 
lower sensors than at the shallow probe depths, suggesting a decrease in surface 
disturbance with radial distance from the probe axis. 

At the Apollo 15 probe 2 site, the quality and resolution of the theoretical 
fits are poor. The SMS minimums are relatively large and not sharply defined 
in the 3 parameter space. At the shallower probe 2 sensors, the transient 
histories are significantly irregular, suggesting a particularly Irregular 
distribution of mean temperature perturbation at the surface, Including regions 
of negative mean temperature anomaly near the probe axis. (Early probe 2 
transients are negative.) Such a surface disturbance results in the poor fits 
of the idealized model and the downward bias of the dlffuslvlty estimates. 

Note, however, the large range of < values which produce comparable RMS fits 
(within a factor of two) and which are compatible with the probe 1 results. 

At the Apollo 17 site, the deeper emplacement of both probes results in 
a more uniform integrated response to the surface disturbance. Consistent 
values of the disturbance parameters B and R produce the best RMS fits at all 
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depths. In all c&s^s, residuals are within the precision of the absolute 
temperature data. Probe 2 diffuslvlty estimates are slightly less than at 
probe 1, a result consistent with the annual wave < deductions. The 
slightly lower diffuslvlty estimates based on the transient fits are pro- 
bably due to the low conductivity layer near the surface. In view of the 
assumptions and approximations Inherent In the Idealized transient fits, 
the bulk diffuslvlty results are considered consistent with and supportive 
of the more reliable annual wave deductions. 

3. Regolith Thermal Conductivity from 15-230 cm: Results and Discussion 

Analysis of the annual wave and astronaut-induced transient from the 

long-term data yielded bulk thermal diffuslvlty results for each of the 

four probe sites. Conversion to effective thermal conductivity values 

required the input of bulk density and heat capacity. The specific heat of 

a number of lunar particulate samples has been measured (Hemingway et al. , 

1973) . Results indicated substantial variation with temperature but almost 

no variation from sample to sample. A value of 0.67 W sec/gm*K, measured 

at 250*K was utilized for all probe depth conductivity calculations. 

The bulk density of drill core samples was reported by Carrier et al. 

(1974) . The results from the Apollo IS drill core indicated densities that 

range from 1.7S to 1.84 g/cm^ from depths of 40-160 cm. Tue to disruption 

of the soil during drilling, these probably represent minimum estimates of the 

in situ values. Maximum density measurements on Apollo 15 soils of 1.89 
3 

+ 0.03 g/cm were obtained by Carrier et al. (1974 ). Thus, a reasonable 

range of densities to use for converting the Apollo 15 diffuslvlty estimates 

to thermal conductivities is 1.75 - 1.90 g/cm'. Similar results for Apollo 

3 

17 soils lead to an estimated range of 1.83 - 2.09 g/cm . 

Resultant thermal conductivity estimates based on the long-term data are 
shown in Figure V-C-16 along with the preliminary results based on the 
initial probe cooldown analysis and the heater activated experiments. Con- 
sideration of all the data indicates large scatter between the results of 
different methods at comparable depths. The final choice of the annual wave 
deductions as being the most reliable measurement of the undisturbed regolith 
conductivity was based on considerations of the sampling volume characteristic 
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of the different methods as well as Che inherent uncertainties of each. 

The amount of regolith material influencing the sensor response for each 
of the transients considered varies directly with the -'-a a duration charac- 
teristics of the particular transient. In order of Increasing time duration, 
the analyzed transients included; 

a. Heater-activated transients (36 hour duration) 

b. Initial cooldown equilibration (approximately 45 days) 

c. Annual wave propagation (one year period) 

d. Astronaut-induced disturbances (approximately three to eight years 
to reequilibrate) . 

The short term heater-activated experiments sampled tho adjacent regolith 
only to distances of a few centimeters from each heater location. For *his 
reason, the modeled sensor response, upon which conductivity deductions were 
based, was particularly sensitive to the designated thermal linkages of the 
probe and borestem near each gradient sensor. In addition, a significant 
fraction of the regolith material sample included the contact zone of material 
disrupted during the drilling process, casting doubt on the representativeness 
of these measurements. 

The conductivity deductions based on the cooldown data were representative 
of primarily the first 10 cm of regolith adjacent to the borestem. However, 
the uncertainties of these. measurements, in general a factor of about two, 
were directly proportional to the uncertainties of the initial excess heat 
estimates. It is to be noted, however, that the assumption of negligible 
drill heating effects (the minlmun initial temperature case), led to conduc- 
tivity results couslstent with the annual wave deductions. 

The probe-borestam-regolitb numerical model demonstrated that the probe 
sensor response to transients of long period (one year or move) is controlled 
entirely by the bulk thermal properties of the regolith. Variations in 
modelled properties characterizing the thermal linkages of the probe-bores tern 
system and disrupted contact zone had no effect on the probe response to the 
long period transients. Thus, only uncertainties related to the precision of 
component identification (for the annual wave deductions) and the magnitude and 
distribution of the surface forcing function (for the long period astronaut- 
induced transient) limited the accuracy of the long-term analysis. 
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The major uncertainty of the deductions based on the astronaut-induced 
surface disturbance resulted from the necessity to fit the data to a model 
dependent upon three parameters. The accuracy of the least square fits to the 
transient data depended upon the magnitude and distribution of the induced 
mean surface temperature changes as well as the thermal diffusivity of the 
regolith. In addition, the results were influenced by the low conductivity 
surface layer (0-10 cm), especially at the shallowest sensors. It was 
therefore anticipated that the transient deductions would be representative 
of a lower bound for thermal diffusivity, especially at shallow depths. As 
shown in Figure V-C-16, at depths greater than 100 cm the results vary little 
with depth and approach the annual »■ ve deductions. 

As discussed in a previous section, fits to the annual wave propagation 
depend only ou one parameter, the thermal diffusivity. Also, the low 
conductivity surface layer does not affect the results since only the relative 
attenuation between sensors is considered. For these reasons, we consider the 
deductions based on the annual wave attenuation tc be the most reliable of 
the in situ conductivity measurements performed to date on the moon- 
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V-D. STEADY STATE GRADIENTS AND HEAT-FLOW VALDES 
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Refinement of Subsurface Temperature P rof i les 

The Identification and removal of diurnal, annual and transient variations, 
and the corrections for shallow sensor radial fluxes allow reliable temperature 
profiles to be calculated at the Apollo IS probe 2 site and the upper section of 
probe 1. Only the probe 1 lower section results have been previously reported 
for Apollo 15. Figure V-D-l shows the steady state temperature profiles at the 
Apollo IS and 17 sites. Corrections at the Apollo IS probe 1 lower section and 
at both Apollo 17 stations due to transient and annual effects were quite small 
(less than 52) and the results are largely unchanged from those reported pre- 
viously (Langseth et al., 1972, 1974). 

The temperatures shown at 65 cm at the Apollo 17 stations are thermocouple 
measurements which are substantially less accurate than the probe sensor results 
(see Section III-B). With the exception of the Apollo 17 probe 2, 65 cm tempera- 
ture, all profiles are remarkably linear, suggesting a vertical uniformity of 
bulk regollth thermal properties at probe depths. The gradients indicated 
were calculated from a linear least squares fit to the probe sensor results. 

The anomalously low gradient at Apollo 17 probe 2 represents a lover vertical 
heat flow along this section (see Section 7-C) . The gradient Increases signi- 
ficantly at shallower depths and may indicate a local disturbance to the heat 
flow since no comparable change in thermal conductivity is indicated by our 
analysis. 

TABLE V-D-l. Gradients and Heat Flows 




Gradient 

m' 1 

Diffusivity 

10“^cm 2 /sec 

Conductivity 

Heat Flow 
10”^/ cm 2 

Apollo IS, 

probe 1 

1.85 

0.87 

1.07 + 92 

1.98 + 10* 

Apollo 15, 

probe 2 

2.52 

0.74 

0.91 + 10* 

2.29 + 152 

Apollo 17, 

probe 1 

1.35 

1.00 

1.31 + 112 

1.77 + 122 

Apollo 17, 

probe 2* 

0.79 

0.88 

1.15 + 122 

0.91 + 13* 

Apollo 17, 

probe 2 

1.30 

0.88 

1.15 + 122 

1.49 + 162 


* Based on thermometers below 130 cm only. 


Table V-D-l summarizes the gradient, diffusivity and heat-flow values 
for each probe. 
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He at-Flow Variabilit y 

At the Apollo IS site the two probes which are separated by ten meters, 

•6 2 

yield heat flows that differ by about 0.3 x 10 W cm (142 of the upper 
value). The gradients differ by 0.67°K/m. If the subsurface temperature 
gradients are extrapolated linearly to the surface, they yield temperatures 
that differ by about 2.3°K. It is very likely that the difference in gradient 
observed results from lateral variations in extrapolated temperature. The 
extrapolated temperature do as not represent the true mean surface temperature 
but rather the mean temperature at the base of the layer dominated by radiative 
heat transfer (see Section V-A) . Consequently, lateral variations in extrapo- 
lated temperatures result from variations in the detailed radiative charac- 
teristics of the very near surface layer. These affects are discussed further 
in Section VI. 

At Taurus Littrow we see a similar inverse correlation of extrapolated 
temperature and gradient if we ignore the thermocouple reading at 65 cm of 
probe 2. We have no ’direct evidence that this sensor is giving us a false 
reading. However, the thermocouples were not as carefully tested as the 
platinum resistance thermometers and the possibility of a loss of calibration 
cannot be ruled out. Extrapolating only the probe 2 temperatures (those 
deeper than 130 cm) to the surface yields a temperature difference between 

r 2 

the two locations of 1.7*K, the difference in heat flow is 0.86 x 10 W/cm , 

482 of the upper value. The difference in gradient is 0.56 *K/m. The ratio 
of subsurface gradient difference to extrapolated temperature difference is 
approximately the same as that observed at Apollo 15. This suggests that the 
variation in extrapolated temperature between the two probe locations is a 
primary contributor to the variation in heat flow. 

In past reports we have assumed the thermocouple reading at probe 2 (65 cm) 
is valid. Support for this assumption is that the gradient calculated using 
the temperature difference between the thermocouple and the lowermost probe 2 
thermometer is 1.3'K/m, the same as at probe 1. There are two alternative 
explanations for the variation In gradient with depth and for the difference 
between the two locations. The first is to assume that a large volume of high 
conductivity material lies close to probe 2 in the subsurface. In the Apollo 
17 preliminary science report we conjectured that a solid rock may lie adjacent 
to probe 2 (Langseth et al., 1973). The fact that the lcng period transients 



112 


are apparently not affected by this hypothesized rock detracts from the 
explanation. An alternative cause of the substantial decrease In gradient 
below 100 m could be an area of the surface near probe 2 with anomalously 
high temperature. This hypothesis Is supported to some extent by the smaller 
effect of the astronaut disturbance at probe 2 compared to probe 1 at Taurus 
Lit trow, see Figure V-B-2 . At all four probes the astronaut's disruption of 
fie surface resulted in a higher mean surface temperature subsequently. Thus, 
if an area of anomalously high temperature already existed near the probe 2 
Location prior to drilling the hole, the relative effect of the astronaut's 
disturbance would be lessened. 

The above discussion clearly illustrates the problems of interpreting 
only four near surface heat-flow measurements on a surface where sizable and 
inknown variations can be expected. 

The best estimate of the regional heat flow at Hadley Rille is the mean 

>f the two observed values (2.1 x 10~^W cm”*) . At Taurus Littrow we assume that 

the probe 2 thermocouple reading is valid which indicates a gradient of 1.3*K m 

and a heat flow of 1.49 +16X. When this value is averaged with the probe 1 

value a mean of 1.6 x 10 cm”* is obtained. The value at Taurus Littrow is 

probably influenced by its position in a valley. The topographic effect of 

surrounding ridges contributes about 0.2 x 10~^W/cm”* to the heat flow. The 

terrain effects will be discussed more fully in Section VI. When the correction 

is taken Into account, our best estimate of regional heat flow at Taurus 

•6 —2 

Littrow becomes 1.4 x 10 W cm . 
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VI. ANALYSIS OF HEAT-FLOW MEASUREMENTS 
AND COMPARISONS WITH EARLIER DATA. 

* 

VI-A. TOPOGRAPHIC EFFECTS 


I2B2E2E& 

Topographic effects on the moon are complicated by the fact that the mean 
temperatures of sloping surfaces are a function of the dip angle and azimuth. 
Variable temperatures over the lunar surface will Induce localized gradients in 
the subsurface and consequently must be taken into account. When correcting 
oceanic floor geothermal measurements for terrain effects, it is virtually 
correct to assume that the seafloor surface is Isothermal. This simplifies 
the analysis. In continental regions of rough topography, the surface mean 
temperature may also be affected by the dip and azimuth of sloping surfaces. 

For example, in the northern hemisphere the northern slopes of hills will be 
considerably colder than those on the south facing slopes (see Blackwell et al., 
in press) . Most heat-flow measurements on earth are made at depths great enough 
so that the effects of mean surface temperature variation are attenuated to a 
few percent of the regional gradients. However, these effects are important 
for the lunar heat-flow measurement because they are so shallow. 

Figure VI-A-1 illustrates the dependence of mean temperature ten centimeters 
below the local surface on dip and azimuth angle. The ten centimeter depth is 
chosen to include the effects of the large mean temperature gradient which 
results from the large temperature dependence of thermal conductivity within 
the first few centimeters below the surface. The mean temperature contours 
shown are based on one dimensional model calculations run over three lunations 
for dip angles of 0-30* and azimuth angles (measured east of north) of 0-360* at 
10* intervals. The contours shown represent the deviation of mean temperature 
from a flat surface. Note that north facing slopes (azimuth angle 0*) exhibit 
decreasing mean surface temperatures as the dip angle Increases. The opposite 
effect occurs for the south facing slopes (azimuth angle of 180*). 

In the above discussion the radiative exchange between surfaces and shadowing 
of surfaces by adjacent terrain were Ignored. Where slopes are large, greater than 
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10* , the exchange between sloping surfaces and shadowing can be significant 
and have a pronounced effect on the mean surface temperature. An example of 
this effect is provided by Hadley Rille. The AJLSEP experiments were set up 
approximately 1.7 km from the eastern rim of Hadley Rille. The rille is about 
1.4 km across and 400 m deep. The cross-sectional profile is roughly V-shaped 
(see Figure VX-A-2) . During the morning, the eastern slope which is facing 
west is in the shadow until the sun rises to nearly 30° above the horizon. 

However, it is receiving radiation from the western slope which, chough partially 
shadowed during the early morning, once the shadow has passed receives solar 
radiation at a relatively high incident angle compared to the eastern slope or 
the surrounding nearly level surface. 

The rather complex history of solar Insolation and radiative exchange between 
the sloping walls of the rille were modeled using a two dimensional finite difference 
model. The geometry of the model is shown in Figure VI-A-2. The computation 
was carried out in two steps: first, the mean surface temperature was determined 

by modeling only the upper several centimeters of the regolich and radiative 
balances at surface nodes during a lunation. Secondly, after the mean surface 
temperature was calculated (the values are shown in Figure VI-A-2) the effect on 
subsurface temperatures was found by relaxation techniques on a two dimensional 
grid of nodes representing the subsurface. The resulting isothermal structure is 
also shown in the figure. Notable features of the calculation are: 

1. The eastern slope actually attains a mean temperature higher than a 
level surface. As a result, heat flows into the surface here. 

2. The western wall of the rille has a temperature 13 s lower than a level 
surface and heat flow through the surface is as a result 3.1 times higher than 
the observed regional heat flow, The isothermal structure under the western slope 
is profoundly affected. 

3. The disturbance produced by the rille decays rapidly away from the eastern 
rim and the disturbance at the ALSEP site is -2% of the regional heat flow. Notice 
that the effect is much larger on the western rim, so that we benefited by the 
relative location of the experiment to this feature. 

At Taurus Llttrow, there were two features near the ALSEP site that could 
produce a prominent terrain effect. The north massif is a ridge of mountains 
forming a sloping face whose base is only three kilometers from the ALSEP site. 




Heat-flow anoaaly (above) and subsurface lsothems (°K) for a cross 
section through Hadley Rllle which passes through the Apollo IS site 
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The ridge is nearly continuous, so chat it can be modeled reasonably accura- 
tely by a bench 1.67 ka high with a plane face sloping 25*. The toe of the 
bench is 3 tan from the ALSEP site. Modeling the feature in this way allows 
the application of the oethod of Lachenbruch (1969) , to deteralne that this 
feature alone results in an 112 increase in the heat flow at the oeasureaent 
site. 

The south massif is about 6 lea froa the site and its shape is approximated 
by a truncated cone. The effect on the heat flow of this feature was modeled 
as a truncated cone, using finite difference techniques in cylindrical coordinates. 
The choice of a coordinate system transforms the problem into a two dimensional one, 
that can be more easily formulated and computed. The model showed the effect of 
the south aassif to be surprisingly small. Its presence increased the heat 
flow only 0.72 at the ALSEP site (see Figure VI-A-3) . Similar models of craters 
showed that terrain features, circular in plan, produce very small effects at 
distances equal to, or greater, chan the radius of the feature. 

Application of the Jeffreys or Birch Technique of Terrain Correction : 

A good approximation to the terrain effect in complex topography of 
relatively low relief can be obtained using an approach proposed by Jeffreys (1937) 
and later developed by Birch (1950). First, the temperature distribution on a 
fictitious horizontal plane, the reference plane passing through the surface at 
the station point, is estimated by a simple linear extrapolation of the tempera- 
ture on the surface to a point on the reference plane using the regional gradient. 
That is, the tesiperature on the reference plane , 

T (r,9,0) - T(r,9,h) + T -h , where (VI-A-1) 

r and 9 are polar coordinates . The vertical coordinate, z, at the reference 
plane is 0, h is the elevation of the topography, where positive is higher than 
Che station point and T is the regional thermal gradient. The disturbance to 
the surface gradient at the station point, r ■ 0, is given by 

ds — 

AT - / ~ • • dr (71- A- 2) 

. or dr 

where T (r) is the angle averaged temperature around a circle of radius r. 
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In practice the effect Is assessed by forcing a grid over a topographic 
map. Because the effect of a given area falls off rapidly with distance from 
r * 0, the area enclosed by the grid lines can be Increased rapidly with r 
without significant loss of accuracy. For each area In the grid, a mean 
elevation Is assigned and a mean temperature Is determined. Note that this 
method can accomodate variable surface temperatures due to the orientation of 
slopes. The map In Figure VI-A-4 for example, shows the surface temperature 
distribution at Taurus Llttrow using the graph of Figure VT-A-1. Once the 
temperature distribution on the reference plane has been approximated, using 
the grid and Equation (VI-A-1) , temperatures are averaged around each annulus 
to determine T and then the Integral in the equation can be evaluated numerically. 

The technique was applied to the topography at Hadley Rllle and Taurus 
Llttrow and the correction to be applied to the measurements, given In Table 
VI-A-1. 


TABLE VI-A-1. Terrain Corrections in Percent 
For Lunar Heat-Flow Measurements Based on the Jeffreys' Approximation 


Site 

Topography 

alone 

Surface Temperature 

Net 

Hadley Rllle 

-6 + 3 

+ 4.2 to + 8.3 

+ 0.3 + 3.7 

Taurus Llttrow 

-10 + 4 

< 0.1 

- 10 






Errors resulting from the application of the Jeffreys' method have been 
Investigated for two dimensional topographies by Lachenbruch (1969). His 
results show that the effects (determined using Jeffreys' method) of positive 
relief features of moderate slope (less chan 30*) on surface gradients at 
points farther chan one relief height from the feature, differ from exact 
solutions by only 31 or less. For example, the effect due to topography 
(calculated using Jeffreys' approximation) are essentially the same as that 
determined using the Lachenbruch bench model discussed earlier. At Hadley 
Rllle there Is a pronounced topographic effect of 61 due largely to the 
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Ms an surface teoperature (*K) in tha region of Taurus Littrow. 
The contours represent topography. 
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Figure VI-A-4 
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Apennine Front south ^ the site. However, this front is facing north and 

consequently, has a surface temperature significant'../ below that of a level 

surface. This large area of below average mean surface temperature produces 

an opposite effect to that of the topography, so that the net effect is not 

significantly different from zero (see Table VI-A-1) . The results given in 

this table are thought to be the best estimate of the terrain corrections 

and they have been applied to the raw data to give the best estimates of he.* 

flow. At Hadley Riile, there is no correction so that the value remains 
_2 

2.1 uw cm . A -102 correction is applicable to the Taurus Littrow measure- 

-2 

ment reducing its value to 1.4 yw cm . 

Using Ra n dom Walk s to_E3t imace_Togographic_Corrections 

We experimented with a novel approach to determine the effect of terrain 
on surface heat-flow values based on the probability theory. The principle of 
this method is that the temperature at a point interim to a volume on whose 
surface the temperature is defined everywhere, can be er imated by averaging 
the temperatures at points on the surface where random walks launched from 
the interior point emerge (see for example, Haji-Sheikh, 1965; Haj 1-Sheik and 
Sparrow, 1967). 

The problems associated with using such an approach aided by a computer 

are: 

1. A large number of walks are required to make a sufficiently accurate 
estimate at a single point, particularly when the point is close to the surface. 

2. Finding a method to determine when a walk had Indeed emerged, was a 
non trivial problem. Associated with this was the problem of finding a con- 
venient way of mathematically representing a complex surface, such as lunar 
topography, for rapid computer calculations of walk positions relative to 
the surface. 

3. Finding a truly random number generator is important because the 
slightest bias in azimuth of emerging points could Introduce large errors. 

Test runs using thousands of points were made to determine the distribution 
of emerging points on a plane surface. The IBM Library Subroutine RANDU was 
found to be unsuitable for our purposes. 

The first problem was largely solved through the use of ? very fast 
machine. To obtain a reasonable accuracy, runs required about 40,000 separate 
walkr from a single point. After Initial experiments in which walks were 
launched from several points along a vertical line near the surface to deter- 
mine the gradient, we found greater efficiency resulted from using the 


122 


random walks to datcrmlne the temperature distribution on a reference plane 
as in the Jeffreys’ technique then, using Equation VT-A-2 to calculate the 
gradient distribution. 

The speed of the walks can also be increased by using the so-called 
floating walk. To perform a floating walk at each step, the length is chosen 
to be the distance to the nearest point on the surface, although the direction 
cosines of the walk segment are randomly chosen. Such walks find the surface 
relatively rapidly. The path was considered to emerge from the volume when 
the direction of a step made a prescribed, small angle with the direction to 
the closest point on the surface. 

The model used to represent lunar problems was approximately a rectangu- 
lar paralle’epiped. The upper surface of the box was the topography represented 
by triangular facets. This is the only means of representing the surface that 
we tried and it may not be the most efficient. It did lend itself to rather 
easy mathematical definition. 

The floor of the box is a plane horizontal surface placed at a depth, 0, 
about five times the largest relief. . The temperature at the bottom was taken 
as T x D where F is the regional gradient. The vertical walls were also planar 
and temperature on the walls was defined using the linear relation 

T(z) - T(h) + (z + h) x (r x D - T(h)) / (D + h) 

where h is the elevation at the top of the wall and the depth, z, is zero at 
the site. 

Various tests of the method were made. The single bench was modeled. 

Random walk determined solutions could be compared with the exact solution of 
Lachenbruch (1969). Good agreement was not obtained due to the effects of 
elevated topography outside the perimeter of the model which is not well 
accounted for by the random walk method. The far topography could be easily 
and accurately accounted for t*slng the Jeffreys' technique. When the random 
walk results were added to this the comparison with the exact solution was 
perfect (+.3T estimation error) with 18,000 walks . 

When the method was applied to the topographies of Hadley Rille and 
Taurus Littrow the results were in good agreement with the results o.' 
the Jeffreys' approximation (see Table VI-A-1) . 
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VI-B. LATERAL MEAN SURFACE TEMPERATURE VARIATIONS 


As we noted In Section VI-A, topographic slopes on the moon cause variations in 
surface temperatures that induce localized gradients in near surface layers. 
Changes in radiative properties of near surface layers can also cause surface 
temperature variations on a local scale which should be considered. The most 
obvious source of an anomalous radiative material are large boulders exposed 
at the surface, which will in general maintain different mean temperature than 
the surrounding regolith because (1) their greater thermal inertia results in 
the retention of more heat during the lunar night, and (2) the conductivity of 
solid rock materials is less temperature dependent than the soil layers so that 
only a small mean temperature gradient is induced in the zone where diurnal 
variations are significant. 

At both the Apollo 15 and 17 sites, very, large mean temperature gradients 
were observed in the upper few centimeters of the soil due to the predominance 
of radiative heat transfer during the lunar day. Temperatures at a few centi- 
meters were 35-40°K higher than the mean surface temperature (Kelha and Langseth, 
1973) . Since this increase depends on the ratio of radiative transfer to con- 
ductive transfer (see, for example, Linsky, 1973) variations in the radiative 
properties of the upper part of the soil could Induce significant lateral 
variations iu mean temperature at a depth of about 10 cm, even if the mean 
surface temperature is uniform. These subsurface variations will in turn 
induce localized gradient perturbations. 

We have little data related to the possible lateral variations of radiative 
properties of the regolith and oust to a large extent depend on Indirect 
evidence and intuitive arguments. 

The effects induced by lateral surface temperature variation on gradients 
measured at depths of 1-2 m will depend on the proximity of the boundary between 
regions of contrasting temperature and the magnitude of the contrast. The 
worst case, as far as geometry is concarned, is a circular area centered over 
the measurement site, similar to the model discussed for the astronaut distur- 
bance in S action V. In Figure VT-B-1 we show the gradients induced at four 
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different depths by a -1®K temperature anomaly in a circular region as a 
function of radius of the circle. Notice cat at each depth the disturbance 
goes through a maximum. At the Apollo measurement depths (0 . 5—2 m) the dis- 
turbance is maximum for a boundary that is 1-3 m in radius. As the boundary 
moves further away, the disturbance decreases rapidly. This illustrates that 
it is principally local variations on the order of a few meters that create ' 
appreciable disturbances. 

This mechanism offers an alternative explanation for the 0.6°K/m gradient 
difference observed between the two probes at the Apollo 17 site. An anomalous 
temperature boundary a few meters from either of the probes which are separated 
by 10 m with a contrast of 2-3® is r^qui -ad based on the curves of Figure VI-B-1. 
The higher temperature would be required above probe 2. Extrapolation of sub- 
surface probe temperatures to the surface shows the surface over probe 2 to be 
1 /°K higher if we ignore the thermocouple measurements at 65 cm. 

Besides the effects of rocks mentioned earlier, recent craters a few 
meters of more in diameter may cause lateral variations in radiative properties 
by disruption and compaction of the soil. In general, the astronauts were ^.ble 
to emplace the probes far from obvious craters and surface boulders. However, 
the effects on the thermal properties of the near surface soil layer may remain 
long after the surface expression of the event has become subdued. The inability 
to exclude local disturbances emphasizes the need for further data which would 
tend to smooth effects at a scale of meters. Critical Information can be 
obtained by remote sensing of infrared and microwave radiation from the lunar 
surface, and the response of the surface to the varying insolation function. 

Such observations can be made by earth-based telescopes but far more valuable 
results' could be obtained from a spacecraft in a low altitude orbit. 
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VI-C. SUBSURFACE CONDUCTIVITY CONTRASTS 

The disturbance of heat flow by topography can be estimated with good 
accuracy because it depends on parameters that are observed. There are other 
possible disturbances due to the subsurface structure of the outer layers of 
the moon about which we have little direct evidence. One such disturbance is 
caused by the refraction of heat at sloping interfaces between bodies of dif- 
ferent thermal conductivity. Several investigators have examined effects on 
the lunar surface heat flow of conductivity structures based on geological 
models of mare basins or ejecta blankets (Conel and Morton, 1975; Mizutani 
and Osako, 1974). The results of these model computations depend on assigning 
thermal conductivities to geological units such as mare basin fill, Impact 
ejecta blankets, and the regollth material 

The existing experimental data show that the thermal conductivities of 

lunar rocks and disaggregate soils that compose the outer layers of the moon 

are mainly controlled by porosity. In Figure VI-C-1 data for lunar rock 

samples and regolith fines are plotted versus porosity. A clear inverse 

relation between porosity and conductivity is seen for rock samples (Group I) 

and disaggregate regolith fines (Group II) . The two groups are not described 

by a single relation between conductivity and porosity. For Group I samples 

the porosity can hardly increase above 251 or else the rock would fall apart 

unless it is an extremely vesicular pumice, or similar rock so chat there is 

-3 -1 -1 

probably a lower limit to the conductivity (about 1-2 x 10 W cm a K ) that 
most rock samples can attain. Mainly the porosity of these samples depends on 
the extent of fracturing of th. original glassy or crystalline matrix, and to 
a lesser extent on the nature and composition of the matrix. Sample measure- 
ments show that among the rock samples there is very little differentiation 
between basalt and breccia samples. The correlation of conductivity with 
porosity suggests that the conductivity structure of at least the upper 1 km 
of the moon depends on the history of impacts that have changed rock porosities 
by fracturing and comminution, and the extent to which the fractures are closed 
at depth as the compressive load increases. 


J — I— 4>4<d 
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,F-B ROCK samples 


GROUP I 


♦ CRUSHED ANORTHOSITE 



POROSITY (%) 


Published conductivity values of lunar materials vs. porosity. For 
laboratory measurements, values at 300®K are shown as points in the 
plot. The doubly hatched box is the range of in situ values deter- 
mined from nighttime cooldown data (Keihm and Langseth, 1973) and 
che singly hatched boxes are the results reported in this paper. 
Rock data are summarized in Horai and Winkler (1976) and soil data 

in Cremers and Hsia (1973) 


Figure VI-C-1 
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In Group IX samples ehe Inverse correlation of conductivity with porosity 
results from the smaller number of intergranular contacts in lower density material. 
The drive-tube samples show the density of regolith materials increases rapidly 
downward reaching densities of 1.8 - 2.0 g/cm^ (Carrier, 1974) below 20 cm. 
Interpretations of cooldown of the surface during lunar night also indicate a 
rapid increase in conductivity with depth in the upper 10 cm (Keihm and Langseth', 
1973). Values at 10-13 cm are similar to those we have deduced from analysis 
of long-term variations presented here. The soils at depths greater than 20 cm 
are near their maximum attainable density and, therefore, near the maximum attainable 
conductivity for Group II samples. The observed profile of conductivity in the 
regolith is no doubt determined by the impact history of particles of a few 
grams and less. The zone of rapid transition from low to higher conductivity 
corresponds to the depth where the net effects of impacts change from disruptive 
to compactlve processes. Because the porosity profile and consequently the 
conductivity profile depend on the random impacts that have struck within a 
given area many times, the conductivity profile in flat regions of the moon 
should be uniform laterally, except near recent large craters. 

The seismic velocity structure of the upper 1 km of the lunar surface 
has been defined by seismic refraction profiles made at three of the landing 
sites. These results show the acoustic velocity structure at all sites is very 
similar to a depth of about 73 m. Cooper et al. (1974) give a comprehensive 
summary of these results. Figure VI-C-2 shows the P-wave velocity structure 
reported at Taurus Littrow. Very low velocities of about 100 m/sec characterize 
the regolith layers. This velocity varied very little between the sites although 
the thicknesses of the layers at Apollo 14 and 16 sites are greater than at the 
Apollo 17 site. The velocities of underlying layers range from 230 to 350 m/sec. 
Only the Apollo 17 seismic experiment defined structure deeper than 100 m. 

The velocities to 1400 m are low compared to lunar rock samples, indicating that 
the materials are highly fractured or are rocky layers interbedded with low- 
velocity materials, even though the traverse gravimeter experiment (Taiwan! et 

3 

al., 1973) suggests these layers have high densities (3.3 - 3.5 g/cm ). The 
seismic results support the evidence that the extent of fracturing and commi- 
nution of near-surface materials determines the seismic velocity structure. 

These factors, which depend on the impact history of an area, also control the 
conductivity structure. 
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By assuming chat the velocity contrast at 30 a in Figure VI-C-2a 
corresponds to the interface between primarily disaggregate materials (Group 
II) and fractured but essentially intact rocky layers (Group I) we have esti- 
mated ranges of bulk conductivity, shown in Figure VI-C-2b andc, are shown 
bracketing temperature profiles corresponding to a uniform vertical heat flow 
of 1.6 uW/cm^. 

The laboratory and in situ measurements of the physical properties of lunar 
surface materials, though still very few, indicate that the conductivity struc- 
ture of the upper 1 km of the moon is laterally tmiform at a regional scale in 
level areas. Consequently, large-scale near-surface conductivity contrasts that 
are steeply sloping are unlikely to occur even at the margins of major geo- 
logical provinces. At depths below 1 km the lower porosity will result in con- 
ductivity contrasts between rock types such as anorthosite and basalt. How- 
ever, these contrasts are small and their refractive effects on surface heat 
flow will be muted by the overlying more uniform blanket of lower porosity 
material. Suggestions that the heat flow at Apollo 15 and 17 may be anomalously 
high by a factor of two or more because they are made at the edge of filled 
circular basins are not supported by existing data. 

Although evidence suggests that at a regional scale the physical properties’ 
structure of the upper 1 km is uniform over level regions of the moon, large 
physical properties contrast can occur at the scale of several meters. The 
large boulders scattered over the surface at all landing sites clearly attest 
to this. One cannot exclude a finite probability that measurements made to 
depths of 2 m are disturbed by refraction effects of such local heterogeneities. 
The change in heat flow with depth observed at probe 2 at Taurus Llttrow no 
doubt results from such a local effect, although it is not necessarily due to 
conductivity contrasts. 
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VI-D. COMPARISON WITH MICROWAVE OBSERVATIONS 


Ie has long baan realizad that cha opacity of tha lunar surficial material 
eo electromagnetic vavas dacreasas with incraaslng wavalangth. Thus, vharaas 
observations of tha moon's infrarad amission indicatad surfaca temperature 
variations, observations at longer wavelengths (in tha millimeter and microwave 
spectral region) yielded information related to the subsurface temperature pro- 
file. The primary difficulty in interpreting the microwave results lay within 
the problem of the interdependence of the thermal and electrical properties and 
their possible variation with depth and temperature. Additionally, most of the 
observations have been subject to large errors in absolute calibration resulting 
in large scatter among all of the data. 

In the absence of more definitive Information of the physical properties 
of the surficial material, the standard procedure for interpreting microwave 
observations has been the following: Assuming a locally plane parallel semi- 
infinite medium, the electromagnetic energy flux (expressed in terms of brightness 
temperature) , leaving the surface at angle 9 ft relative to the local normal , is 
expressed in the form 

z 

• /<, (z'T)sec 9. dz' 

T b (X, t,9 # ) - ( 1 -R(p, 0 a )]sec9 i /< x (z,T) x T(z,t)e"° A 1 dz, (VI-D-1) 

where T^(\,t,9 e ) Is the equivalent brightness temperature expressed as a function 
of the wavelength X, time of lunation t, and the emerging angle 9 ; R(p,t ) is 
the intensity reflection coefficient dependent on the state of polarization p, and 
9 # ; 9^ is the angle at which the observed radiation impinges on the surface from 
below and is related to by Snell's law; <^(z,T) is the power attenuation 
coefficient expressed as a function of depth z and temperature T. 

Equation (VI-D-1) is a generalization of the equation given by Plddington and 
Mlnnett, 1949. Tne Rayleigh-Jeans approximation is assumed valid for the lunar 
situation at microwave frequencies. Equation VI-D-1 also assumes that the 
attenuation coefficient is purely absorptive with no scattering component 
and does not take into account tha possibility of ray-path bending due to possible 
inhomogeneltles in the regolith material. Additionally, Equation VI-D-1 yields 
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eh* equivalent brightness temperature that would be observed by an antenna beam 
sufficiently narrow so that the surface area Included in the beam is small enough 
to be considered planar. In actuality, the observations are made with antenna 
beam widths which necessarily average the radiation over a large portion of the 
lunar disk Including a rang* of emergent angles. 

At longer wavelengths, the entire disk contributes to the observed radiation. 
Krotlkov (1963) has derived a formula to convert disk center brightness temperature 
to disk average temperatures assuming homogeneous regolith properties, the 
conversion formula is primarily a function of the real part of the dielectric 
constant. For wavelengths less than 3 cm, many of the observations have been 
made for the disk center; i.e., the antenna beam width are smaller than the disk. 
Because the intensity contribution away from the disk center varies according to 
a cosine law, the predominant contribution for such observations comes from the 
area near disk center. Thus, these observations can be compared directly with 
little loss of Information to theoretical models governed by Equation VI-D-1 
with ■ 9^ • 0 

Because of the t/eraglag character of the observations and the inherent 
accuracy limitations of the remote measurements, Equation VI-D-1 has been used 
mainly in the simplified form in which the lunar regolith properties are assumed 
homogeneous in depth and independent of temperature. The assumption of homogeneity 
allows the temperature profile T(z,t) to be expressed as a Fourier time series 

in terms of the harmonic amplitudes of the surface temperature variation and the 

1/2 

thermal attenuation coefficient for each of n harmonics, S Q a [nupc/2ic] , where 
u is the lunation frequency and p, c, and < are respectively the density, specific 
heat, aud thermal conductivity of the material (Carslaw and Jaeger, 1959, p. 63- 
70) . Equation VI-D-1 can then be solved for the Fourier components of the lunation 
brightness temperature variation (Muhleman, 1972). The resultant brightness 
temperature amplitude components are related to the surface temperature com- 
ponents through the empirical parameter 6 • B /<, , the ratio of the thermal 

U u A 

to electrical attenuation coefficient for the nth harmonic: 

T. (X) « (1- R(p) ] T / [1+25 + 25 2 ] 1/2 (VX-D-2) 

oti ana 

where T^ q (X) is the nth component of the observable brightness temperature at 

disk center and T is the nth component of the surface temperature diurnal 
a 

variation at disk center. 



The microwave observations have been analyzed In terms of a mean temperature 

and the amplitude and phase shift of the fundamental frequency. The standard 

1/2 

procedure has been to derive the fundamental parameter, 6^ • (wpc/2k) /<y from 
the observations of T^o^bl’ which Is the ratio of the mean to the amplitude of 
tua first harmonic: 

V T u • «,/ T i> a + 2S i * 2 S i Z,l/2 

The ratio T^/T ^ ^ can be determined directly from the observations ’with 
inaccuracies due to calibration errors minimized. The only sources of data for 
the mean to first harmonic ratio of surface temperature prior to the lunar 
missions were the infrared observations during lunar night. Daytime surface 
temperature variations are controlled mainly by the solar insolation, surface 
topography and the radiative properties of the surface. Troitsky, 1965, adopted 
a value of 6^ ■ 2. OX based on measurements between 0.1 and 3.2 cm and a value of 
T q /T^ * 1.5. Corrections for moonwide averaging effects ware applied. 

The interpretation of the empirical parameter 6^, iu terms of electrical 
properties of the surficial material c&mot be done unambiguously without an 
independent determination of the thermal properties. Prior to the more definitive 
measurements of the thermal properties from the Apollo missions, infrared obser- 
vations of eclipse and lunation cooling curves were the primary source of thermal 

property Information. The controlling parameter for the cooling curves is 
-1/? 

Y • (icpc) ” (see, for example, Jaeger ,1953). 

It was found that the eclipse data controlled by properties of the upper 

few centimeters at most, required a value of Y much higher than the lunation 

cooldown data which, because of its longer duration, reflected thermal properties 

2 

from deeper layers down to about 20 cm. Values in the range y • 250-350 (cm K 
-1 - 1/2 

V sec ) characterized the eclipse measurements and y ■ 100-200 characterized 
the lunation nighttime curve. Selection of an appropriate thermal parameter y, 
based on the Infrared observations, then permitted an Independent determination 
of the thermal absorption length once values for d and c were estimated from 
other information 


JjxrfzyH ZpcY‘ 


(VI-D-4) 
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The telescopic observations reveal that almost no variation in the micro- 
wave brightness temperature of the lunar disk can be detected at wavelengths 
beyond about 5 cm. If one, therefore, assumes that the electrical attenuation 
coefficient is independent of depth for the effective emitting layers correspon- 
ding to X > 5 cm, and that the thermal conductivity k, is uniform in depth end 
independent of temperature for these same layers, the disk center brightness 
temperature as a function of wavelength can be written from Equation VI-0-1 
as 

T. (X) - [1 - R(p) ] x K ./* (T + q • z/k)e“ K \* dz (VI-D-5) 

D A o 0 

where q is the heat flow from the interior and q/k the uniform temperature 
gradient. The integral is elementary and yields 

T b (X) - [1 - R (p)](T o + (q/k)fc x ). (VI-D-6) 

If we now assume that the reflection coefficient (or equivalently, the 
real part of the dielectric constant for the lunar situation) is independent 
of wavelength, we can writ^ the spectral gradient directly by differentiating 
(VI-D-6) 


<5T b /5X - (1 - R (p) ] (q/k) (<5Z X /6X) (VI-D-7) 

(since < x * l/£- x ) . 

For the conditions known to exist in the lunar regollth, and if we neglect 
scattering, it can be shown (Muhleman, 1972) that the electromagnetic absorption 
length may be written in the form 

l x - X/(2JIe l/2 tanA) - C] X, (VT-D-8) 

where e • the real part of the relative dielectric constant and tan A * 2 a/ tv 
is the loss tangent expressed in terms of the dielectric constant £ , and the 
electrical conductivity <7 at frequency v. 

The form of Equation (VI-0-8) , together with the assumption that tan A be 
independent of wavelength (i.e., 0 proportional to v) , has actually already been 


f 
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utilized In the assumption that <5^ is directly proportional to the wavelength. 

Incorporating Equation (VI-D-8) into (VI-D-7) yields the heat flow directly: 

q - k (ST b /<5A) c’ 1 [1 - R (p)f 1 (VI-D-9) 

where c^ - 1^2IIe^ tan A) is assumed wavelength independent. The spectral 
gradient, 5T./6X, can be computed directly from the observations at different 

D 

wavelengths. Then, by using values of k and deduced by methods described 
previously and an R(p) value based on radar studies, an expected value of the 
heat flow can be calculated. 

It is to be noted that independent deductions of the constant c^ and the 
thermal conductivity k could only be made from remote observations of the infra- 
red and short centimeter emissions which yield information characterizing only 
the upper 'v 20 cm of the lunar regollth. Thus prior to the lunar landing missions, 
extrapolation of near surface properties to the greater depths was required for 

heat flow determinations to be made. For example, Baldwln(1961) estimated an 

-6 1 

upper limit of q « 1.0 x 10 W cm* based on a homogeneous thermal model proposed 
by Jaeger (1953), a derived lower limit of the electrical absorption length ly 
and a brightness temperature measurement at 168 cm. Russian investigators, using 
a new calibration technique called the ’artificial moon’ method, made a series of 
high precision observations at wavelengths between 0.4 and 70 cm in 1961 through 
1964 at the Radiophysics Research Institute (NIRFI) in the Soviet Union. Tikhon- 
ova and Troitsky (1969), postulating two layer models with a high density, high 
conductivity rock layer overlain by a less dense, less conductive regolith, and 
using thermal and electrical material properties deduced from the Infrared and 

short centimeter (0.1 < X < 3.2 cm) observations, estimated a lunar heat flow 

-6 2 

in the range 2.7 -4.6 x 10 W/cm based on NIRFI observations. Additionally, 
they suggested that a regollth layer thickness between 3 and 13 m best fit 
observations of a decreasing spectral gradient for X >30 cm. 

Observations made with the Areclbo radio telescope in 1969 suggest the 
possibility that at wavelengths between 69.8 and 406.5 cm the spectral gradient 
may be negative, i.e., the observed brightness temperature decreases with 
wavelength. However, the calibration difficulties involved in the detection of 
lunar thermal emissions at these long wavelengths make any strong reliance on 
these data extremely dubious. This point will be examined in more detail. 
Nevertheless, Salibury and Fernald (1971) assuming a power absorption length of 
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i^ m 40 \ , and accepting the observed negative spectral gradient as real, esti- 
mated a thermal temperature gradient of - 1/3 Km for depths corresponding to the 
observation wavelength emissions. They interpreted this result in terms of heat 
escaping the moon not by thermal conduction, but by the transpiration of volatiles 
from the lunar interior through the porous surface material. We postulate that, 
if the negative spectral gradient is real, it is not produced by a negative 
thermal gradient, but by the processes of scattering from debris layers tens of 
meters below the surface. Only che near linear portion of the brightness 
temperature spectrum (5 < X <25 cm) ca.. be utilized in che determination of che 
lunar heat flow. Observations at longer wavelengths are most valuable in the 
interpretation of che structure of che deeper surface layers. 

E lectrical Abso rp tion L ength 

In Figure VI-D-1, data are shown for the parameter l^/\ for a number of 
different wavelengths between 0.1 and 63 cm. The data of Ade et al. (1971) 
correspond to a density of p« 1.25. The data of Clegg et al. (1972) are t'or 
a density of p • 1.13. The values at 0.5, 1.25, and 3.2 cm were calculated 
from Equation VT-D-3 using the data of Bassett and Shackelford (1972) for the 
dielectric constant at P- 1.3 and a reported value of tan A- 0.015. The span 
of absorption lengths at 63 cm were taken from che data of Gold et al. (1973) at 
a density of 2.0 g/cm^. The particular densities chosen correspond to the esti- 
mated densities at che depths of tne primary emitting layers for a given wave- 
length. It can be seen from Figure VI-D-1 that the validity of Equation Vl-D-8 
assuming a nearly constant value of £ over the relevant wavelength range, 
requires a wavelength dependence of che loss tangent with a variation of can A* 
0.040 at \ • 0.1 cm to can A « 0.002 at X- 63 cm. Such a decrease in tan A 
with wavelength is not expected based on measurements of dry dielectric materials 
which are characterized by a nearly frequency independent loss tangent for fre- 
quencies above 10^ Hz (Strangway, 1969). Gold et al. (1973) have pointed out 
that even very slight atmosphere contamination may significantly increase the 
losses for a given sample measurement. Such contamination could be responsible 
for the low absorption lengths measured at the millimeter and short centimeter 
wavelengths. Clegg et al. (1972) postulate that scattering processes may be 
responsible for che absorption behavior at the millimeter wavelengchs. Husky 
(1973) points out that, at millimeter wavelengchs, one may be seeing the low 
frequency tail of far infrared resonances and chus a frequency dependence of the 
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Ratio of the amplitude of variation to the mean brightness temperature over a 

I* ?? Ve J? n *5 h ‘ P 16 curves X ’ 7 are resu ^s of theoretical calculations 
based on Apollo 15 thermal properties and 7 different dielectric property pro- 
files. 1) data of Bassett and Shackelford (1972). 2) l\ - X , 3) l\ • 3 \ 

^ 4 X " f ** 5) Zx ■ 12 X . Curves la used the results of Ade et al., 1971* 

and lb the results of Clegg et al. (1972) for X - 0.1 - 0.2. 


Figure VI-D-2 
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absorption length. It will be shown that the absorption behavior shown In 
Figure VI-D-1 fx.3 the remote observations quite well. Relatively low absorp- 
tion lengths are required to fit the 0.1 - 3.2 cm observations and significantly 
higher absorption lengths, near those measured by Gold and coworkers, are required 
to fit the 5-25 cm wavelength emission observations. 

Compatibility of the Measured Physical Properties and the 0.1 to 3. 2 c m R emote 

Observations 

Millimeter and centimeter microwave emissions originate in the upper few 
centimeters of the lunar surface material where the thermal and electrical 
properties vary significantly with depth and there are large temperature fluctu- 
ations during the lunation cycle. Therefore constant property models cannot be 
used and Equation (VI-D-1) must be solved numerically, using physical property 
relationships derived from the available data. To obtain the temperature profile 
T(z) , as a function of time of lunation, the heat conduction equation was solve a, 
using the density and thermal conductivity profiles deduced for the Apollo 15 
heat flow site. The functional relationship between 2^ and density was taken 
from a linear fit to the calculated values of 2.^/X vs density based on the 
data of Bassett and Shackelford (1972) . From these data we derive 

3 

l^/\ • 8.0 - 1.5p (Ping/cm ). 

The curves shown in Figure VI-D-2 of the ratio of the fundamental amplitude 
to the mean T^/T^ , as 4 function of wavelength result from the numerical 
solution of Equation (VI-D-1) with 0^ * 9^ * 0. The observational data points 
are taken from the review paper of Hagfors (1970) . The comparison with disk 
center observations is quite good, considering the scatter in the data. The 
fit for the short millimeter wavelengths can be improved by incorporating the 
data of Ade et al. (1971) or Clegg et al. (1972), for X - 0.1 - 0.2 cm into 
the numerical solution of Equation (VI-D-1) ( see curves 2,3 in Figure VI-D-2). 
Curves 4-7 result from electrical property models in which 2^/X does not vary 
with depth (density) and are shown to illustrate sensitivity to the T. ,/T 

Ox oo 

ratio of the electrical absorption length. 

Shown in Figure Vl-D-3 are the curves using the same 2^/X functions of 

density but incorporating the thermal conductivity and density profile derived 

for the Apollo 17 heat flow site. It can be seen that the agreement is still 

good, although the higher T, ,/T. values obtained at the longer wavelengths 
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imply chac slightly longer absorption lengths are required to best fit the 
observational data. In fact, model 4 fits the observations best using the 
* 8.0 for 0.5£l£3.2, a value only slightly higher than the values 
of Z^/\ obtained from Bassett and Shackelford's data. The fact that the 
thermal properties derived from in situ data and the electrical absorption 
lengths measured for returned lunar samples produce good agreement with the 
0.1 to 3.2 cm microwave data lends strong support to the representativeness 
of these physical property determinations for a large portion of the lunar 
surface to depths of 10-20 cm. 


T he 5-500 cm Emission Spectrum. Heat Fl ow D e du ct ions a nd _ t he Mean Depth o f 


Th e Lunar Regolith 

Diurnal variations in microwave emission for X >5 cm are very difficult 
to detect. This indicates that the primary emitting layers for wavelengths 
greater than 5 cm lie below the depths at which substantial diurnal variations 
and mean temperature changes due to the temperature dependence of the thermal 
conductivity are felt. The subsequent increase in brightness temperature with 
wavelength up to X a 30-50 cm is thus due to the steady state thermal gradient 
in the lunar regolith. 

Table VI-D-1 shows tne observations upon which our deductions will be 
based. The data up to 70 cm wavelength is based on the artifical moon techni- 
que. The observations at 69.8 cm and three longer wavelengths were made at 
Arecibo by Salisbury and Fernald (197i). The stated accuracy of the Arecibo 
data is highly questionable. Problems of absolute calibration and substantial 
sky background have discouraged most radio astronomers from attempting measu- 
rements of the lunar thermal emission at wavelengths beyond 70 cm. The Arecibo 
data are included only as tangible evidence for the possibility of a negative 
spectral gradient existing at meter wavelengths. 

The primary Inferences drawn here depend only on the reliability of the 
Soviet data. The Soviet measurements represent the only observations made at a 
number of different wavelengths in the microwave band using a consistent calibration 
method. The Soviet data agree quite well, in absolute temperature levels, with 
models based on in situ thermal property determinations and measured subsurface 
te.v°.rature levels, expecially when the large scatter among ocher observations 
is considered. It must be pointed out, however chat the importance of observed 
brightness temperature measurements to inferences of global heat flow and physical 
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TABLE VI-D-1. Lunar Disk Average Brightness Temperature Observations at Wave- 
lengths Between 3 and 500 cm 


Wavelength 

(cm) 

Brightness/Temperature 

(K) 

Error estimate 
(K) 

3.2 

213 

+ 4 

9.6 

218 

+ 5 

11.0 

219 

+ 6 

14.2 

221 

+ 5 

20.8 

225 

+ 5 

25.0 

226 

+ 6 

30.2 

227 

+ 7 

35.0 

236 

+ 8 

36.0 

237 

+ 6 

40.0 

224 

+10 

50.0 

229 

+ 6 

58.0 

224 

— 

70.16 

225 

+ 7 

69.8 

256.9 

+ 7 

152.7 

241.9 

+ 5 

'269.1 

228.2 

+10 

406.5 

214.8 

+14 


structure of the upper crust lies in the change of brightness temperature with 
wavelength. The spectral gradient between 5 and 20 cm wavelengths, upon which 
any Inferences regarding global heat flow must be based, cannot be considered 
definitively accurate from a statistical viewpoint because of the short range 
of wavelengths relative to the absolute temperature uncertainties (see Table 
VX-D-1) . In this regard, conclusions drawn of lunar global heat flow must be 
considered tentative until additional, more refined measurements can be made. 

Inferences drawn regarding the average depth of the lunar regolith layer 
depend on the turnover in the spectral gradient between 30 and 70 cm wavelength 
suggested by the Russian data. The sharpness or precise location of the ipectral 
peak is not as important to the implication of regolith thickness as the accuracy 
of the 5-20 cm spectral gradient is to Che estimates of global heat flow. 

Because of the scatter in the four longest Russian observations, a more precise 
positioning of the spectral peak cannot be made, and only a range of feasible 
regolith thickness can be estimated. 
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Two features of the Russian data stand out. The near linear Increase in 
brightness temperature with wavelength between 5 and 20 cm (see expanded scale 
of Figure VI-D-5) is consistent with a homogeneous regolith over the depth range 
from which the 5-20 cm waves originate. If regolith homogeneity were maintained 
to depths on the order of 100 m, a continuation of a linear emission spectrum 
to meter wavelengths would be expected. The observations, however, indicate a 
flattening of the emission curve beyond 20 cm followed by a possible decrease 
In brightness temperature at longer wavelengths. One explanation for a flattening 
of the emission curve is the possibility of increasing limb darkening effects at 
longer wavelengths. If the regolith density increased significantly over depths 
relevant to the 5-30 cm regolith material reaches a state of near maximum com- 
paction within the upper 10-30 cm of the surface, lending support to the assump- 
tion of uniformity with depth for all but the very near surface material of the 
regolith (cf. Carrier et al. (1973)). Also, the thermal conductivity measure- 
ments at depths in the range 91-233 cm at the Apollo 15 and 17 sites showed only 
a slight trend toward higher conductivities at greater depths (Langseth et al. , 
1973) . Most of the variation was more likely due to the presence of centimeter 
size rock fragments at different locations along the heat flow probes. Such 
small fragments would only negligibly affect the propagation of the longer wave- 
length emissions. 

It should be pointed out that the seismic data (Kovach and Watkins, 1973) 
indicate that no widespread solid rock layers exist within 20-200 m of the lunar 
surface. However, the multilayered models are utilized only to demonstrate the 
feasibility of reflective scattering processes contributing to the observed 
negative spectral gradient for \ > 30 - 70 cm. The actual scattering is most 
likely due to a random distribution of large rocks and boulders within the 
rubble layers. The simplified models can, however, provide information related 
to the extent of the scattering and the average depth of layers with significant 
amounts of boulder size rocks. 

The solution for the surface intensity observed for a two-layer model of 
regolith material overlying bedrock is presented in the paper of Tikhonova and 
Troltsky (1969) . The general problem of N emitting layers is more complicated 
but can be solved using basically the same ray-tracing approach as was illus- 
trated by Tikhonova and Troltsky. 
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Figure VI-D-5 


Wavelength dependence of the steady-state disk average lunar brightness temperature. Data points are taken) 
from earth-based observations* The theoretical curves shown are based on multilayer models of the moon’s * 
upper crustal structure. The regollth thermal gradient utilized, 1.3 K/m Is based on a mean heat flow of 
3 xlO“ 2 w/m 2 and a mean thermal conductivity of 2.3 x 10" 2 v/att m-lfT* measured at depths below 1 m at the 
heat flow sites. 
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The general method of solution is to replace the continuous distribution of 
sources in a given layer with two discrete sources, one located at the top of 
the layer radiating upward and one at the bottom of the layer radiating downward. 
Multiple reflections are taken into account. The next step is to calculate Che 
power attenuation of each source as its energy propagates to adjacent layers 
above and below. These layer- to-layer transmission factors may then be used to 
calculate the contribution to the observable intensity emergent at the surface 
from each of the discrete sources. The total observable intensity is simply 
the sum of the intensities due to each of the sources. The multilayer solution 
for a given emergence angle 9^ will depend on the absorption coefficient 
real dielectric constant £ (through the reflection coefficient appropriate to 
each boundary) and the temperature distribution within each layer. 

As discussed previously, the earth based observations for wavelengths 
greater than 5 cm necessarily average the emission over the whole disk. Thus, 
for a valid comparison of theory and observation, theoretical disk average 
brightness temperatures must be computed. 

* e 

\ (X) " » iVfi ///z /2 T b x cos '^ , cos 2 <f> 4* 44, ^ VI - D - 10 > 

where f, (X) is Che disk average theoretical brightness temperature observable 

D 

from the earth and T, ($,4*,X) is the theoretical brightness temperature observable 

0 

at lunar latitude $ and longitude t|> (relative to the subearth point) at an angle 

3 ■ cos” X (cos $ cos t|> ) relative to local normal, 

e 

For each particular multilayer model, solutions for T^(4>,tJ>,\) were carried 

out for 81 different sets of lunar coordinates at 10* intervals in and ( p . 

Regolith and subsurface physical properties and layering configurations were 

considered invariant with $ and \p . Temperature profiles at different lunar 

coordinates were adjusted for latitude. 

In all of the theoretical models considered, the physical properties and 

spacing of the rock layers below the regolith were the same. The general 

geometry and physical properties which remained unchanged from model to model 

are shown in Figure VI-D-4. Values of £ ■ 6.5, • 10X and q/k • 0.027K/m were 

chosen for the rock layers based on the measurements of Gold at al. (1973) for £ 

-6 2 

and JZ.^; Horal and Simmons (1972) for k; and a mean heat flow of 3 x 10 V/cm 
based on the Apollo 15 and 17 measurements (Langseth et al., 1973) A value of 
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e ■ 3,Q wta chosen for all che soil layers including che regolich based on the 

3 

data of Gold et al. 0-9731 at a density of p * 2.0 g/cm and che blstacic radar 
measurements of Tyler OL9.681 . The particular spacing used for the rock soil 
alternating layers vss chosen to produce the negative spectral gradient of *\» 

■"10 K/a observed at Arecibo for wavelengths beyond 70 cm. The particular 

«* * 

physical parameters chosen fob che subregolich reflecting layers are noc critical 
In chat a different set of parameters can be made to produce che - 10 K/m 
spectral gradient simply by varying the thicknesses of the reflecting layers. 

The spaclngs and physical properties of the subregollth layers have negligible 
effect on the 5-20 cm spectral gradient and location of the emission peak. 

Only the regolich properties and thickness substantially affect the emissions at 
wavelengths where the spectral gradient is positive. 

In all models, che local temperature profile was considered to cake the form 


T (*,♦) - T (♦) + / * J dr (Vl-O-ll) 

O OK 

If each layer Is homogeneous, the temperature in any layer n, can be written 


T a <z,$) 


n-i 

T (4) + I 
0 1-1 


a a + i 

k i k 
K i 1 n 


(z - 


d n-l> 


(VI-D-12) 


where T ft (4) " tr.e adjusted mean surface temperature at latitude p , 
a^ - che thickness of Isyer i, 

d^ > the depth to the lower boundary of layer 1 from the surface, 
k^ - the thermal conductivity of layer i, and 
q - the steady state heat flow. 

Mean surface temperature variation with latitude 4 was calculated using the 

thermal model established for the Apollo IS heat flow site(Keithm et al., 1973) 

and the solar insolation function appropriate to a given latitude, 'Adjusted' 

mean surface temperatures T (4) were calculated by adding tr.e mean temperature 

0 

difference observed in the upper few centimeters at Hadley R. 4 lle (*v 4GK) to the 
surface mean calculated from the Hadley Rllla thermal model applied at various 
latitudes. The large mean temperature gradient is confined to che upper few 
centimeters of surface material and has negligible effects on calculations of 
lunar emission for X >5 cm. Thus, the approximation of a constant regolich 
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Figure VI-D-6 
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The same data and models as ahc<m In Figure VI-D-5, with the exception that the models Incorporate a re 
gollth thermal gradient of 2.0 K/m, *\* 50% higher than the mean gradient measured at meter dephs at the 
Apollo 15 and 17 heat flow sites. 
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thermal gradient is still valid. However, the total mean temperature difference 
must be added to the actual mean surface temperatures to avoid a constant error 
^ 40 K in the subsurface temperature profiles 

Theoretical results are shown in Figures VI-D-5 to VI-D-8 with the Russian 
and Arecibo data plotted for comparison. The inset of each figure repeats the 
5-2$ cm spectra on expanded scale to show more clearly the data relevant to global 
heat-flow determinations. The curves of Figure VI-D-5 are produced using the 
following regolith parameters: 

q/k « 1.3 K/m ; 50X [models of Figure VI-D-5] 

These values are based on the average gradient measured at the Apollo 15 and 
17 heat flow sites (Langseth et al., 1973) and the absorption length measurements 
on Apollo 15 and 16 samples (Gold et al., 1973). The three models shown are dis- 
tinguished by varying regolith thicknesses. For a fixed absorption length and 
thermal gradient, the regolith thickness controls the level and wavelength at 
which the emission maximum is reached. At wavelengths beyond the wavelength 
of maximum brightness temperature, the effective emitting layers lie substantially 
within the rubble or reflecting layers. Thus, energy loss through multiple 
reflections limits the observable longer wavelength Intensities despite the fact 
that the effective emitting material is at higher temperatures. It can be seen 
from Figure VI-D-5 that a regolith thickness of ^ 20 m produces a disk average 
brightness temperature maximum between 40 and 50 cm wavelength, compatible with 
Che observations. 

In the same vein, it can be seen from the inset of Figure VI-D-5 chat the 
regolith thickness also controls the range over which the spectral gradient 
remains linear. At wavelengths for which the effective emitting layers lie 
substantially above the reflecting layers, the spectral gradient remains nearly 
constant in accordance with Equations (VI-D-7) and (VI-D-3). In this respect, 
the models of Figure VI-D-5 indicate that a regolith thickness of 20-30 cm is 
required to maintain a near linear emission spectrum out to 20 cm wavelength. 

The linear portion of the spectrum must be used to estimate global heat 
flows and regolith properties. The magnitude of the near linear spectral gradient 
is controlled by the thermal gradient q/k, and the electrical absorption length 
2^. All of the models shewn in Figure VI-D-5 produce spectral gradients that are 
less than that observed by the Russians over the 5-20 cm wavelength band. 
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Figure VI-D-7 






The same data and models as shown In Figure VI-D--5, with the exception that the models utilize a power fit 
to the regollth absorption lengtl from the data of Figure VI-D-1 




Equation VT-D-7, used as an approximation to the disk average spectral gradient, 
indicates that a larger spectral gradient requires either a larger thermal gra- 
dient or longer absorption lengths. 

The emission curves of Figure VI-D-6 are produced using a thermal gradient 
of 2.0 K/m, larger by a factor of one half than the mean value measured at the 
heat flow sites. The 5-20 cm spectral gradient now matches the Russian data 
quite well if a 20-30 m regollth thickness is used, although the absolute bright- 
ness temperatures lie near the upper error estimate of the observations. The 
absolute brightness temperatures are dependent on the global distribution of mean 
surface temperature assumed in the models. Errors of 5-6 K could result from 
extrapolation of the Apollo 15 and 17 observations. The most serious discre- 
pancies of models 7 and VI lie in the spectral range 20-70 cm where the high 
brightness temperatures reached exceed the Russian observations by 15-25 K. 

A method for increasing the 5-20 cm spectral gradient without significantly 
Increasing longer wavelength brightness temperatures is illustrated by the models 
of Figure 7I-D-7. In these models, an electrical absorption length of ■ 80A 
is assumed. This value is near the waft mum measured by Gold et al. (1973) at a 
68 cm wavelength on lunar samples. Model IX produces good agreement with the 
Russian spectral gradient for 5 < A <20 cm although a mean regollth thi.-.kness of 
at least 30 m is now required to maintain the linearity of the spectrum out to 
20 cm wavelength. This is due to the fact that, because of the longer absorption 
lengths, the effective emitting layers begin to lie within the reflecting rubble 
layers at shorter wavelengths. Such a regollth thickness is two to three times 
greater than that inferred from active seismic data (Kovach and Watkins, 1973). 

It is to be noted that all of the models shorn in Figures VI-D-5, 6 and 7 
apply electrical absorption length data at 68 cm wavelength to the entire spectral 
band from 5 to 500 cm. It is quite possible that the absorption constant i^/A , 
varies substantially over this large wavelength band as suggested by the data of 
Figure VI-D-1. To examine the feasibility of a systematic increase in V A 
with wavelength, we have fit a power curve to the data of Figure VI-D-1 to 
obtain 2,^/A * 5.8A 1 ' 1 * 8 . The resultant theoretical curves are shown in Figure 
VI-D-8. Mote that the 5-20 cm spectral gradient is much less than the obser- 
vations. A therual gradient of ^ 3.0 K/m would be required to match the data, 
indicating either a global heat flow more than twice that measured at the Hadley 
Rllle and Taurus Llttrow sites or a mean regollth conductivity less than half the 
mean value measured at the heat-flow sites. If the lnterpolatlve formula for 
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2.^/X does apply to the average regolith material, however, a mean regolith 
thickness of 10-15 m would be sufficient to match the wavelength at which the 
observational maximum occurs (30-70 cm). Such a regolith thickness is in good 
agreement with the seismic data. 

One point is quite clear from all of the models. If the Russian data are 
accepted at face value, the regolith thermal gradient measured at the two lunar 
heat flow sites lies at the lower end of the range of likely moonwide average 
values. Higher gradients could be explained without increasing the heat flow 
if over much of the nearside of the moon, regolith conductivities were 30-50Z 
lower than those measured at the heat-flow sites. Slightly less compact rego- 
lith material than that found at the mare heat flow sites would be characterized 
by lower thermal conductivities and longer absorption lengths, effectively 
increasing the spectral gradient between 5 and 20 cm wavelength. 
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VII. SIGNIFICANCE OF THE HEAT-FLOW RESULTS 
A. HEAT-FLOW AND SURFACE THORIUM ABUNDANCES 

The heat flows reported In Section V-D Indicate that the measurement at 
Hadley Rllle Is significantly higher than chat observed at Tautus Littrow. 

There is no reason to prefer one measurement relative to the other as repre- 
sentative of the global heat-flow based on an assessment of local disturbances 
(see Section 71) . The orbiting gamma-ray experiment flown on Apollo 15 and 16 
(Metzger et al., 1974; Bielefeld ec al., 1976) showed the western mare plains 
of the moon's nearside to have relatively high abundances of thorium at the 
surface. Over the eastern margin of Imbrium, the surface abundance of thorium 
is over two times higher than over Taurus Littrow which suggests chat the 
higher heat flow at Hadley Rllle may be related to the higher thorium abundances. 
These high relative abundances are probably confined to a near surface layer. 

If the abundance detected by the gamma-ray experiment represented abundances 
uniformly distributed in a 60 km chick crust, then the very high surface 
concentrations such as those observed over the Archimedes and Fra Mauro regions 
cannot be representative of average lunar crustal abundances as they would lead 
to unreasonably high heat flows and near melting temperatures at the base of 
the crust. In Figure VII-A-1, the dashed line shows the heat flow at the sur- 
face of a 60 km crust as a function of thorium abundance assuming a K/U ratio 
of 3200 and a Th/U ratio of 3.7 higher than the earth's average heat flux! 
Furthermore, the temperature at the base of a 60 km crust would be nearly 1200°C 
which would place material there very close to the solidus. Consequently, 
the high thorium concentrations observed by the gamma-ray experiment over 
Oceanus Procellanum must be a surficial enrichment in thorium rich materials. 

Interpretations of geothermal measurements on earth have revealed a 
persistent correlation of heat flow with the surface abundance of heat producing 
isotopes. The correlations are most striking in regions that are tectonically 
older (Mesozoic or older; see, for example, Birch et al., 1968 and Lachenbruch 
1970). The correlation within a large geographic region is linear and described 
by the relation q » q* + D A(o) ; where q* is the intercept and A(o) the heat 
production. When a line is best fitted to the heat flow and heat production 
data, it extrapolates the zero heat production axis to yield a positive heat- 
flow value, q*, interpreted as the heat flow into the base of the crust. 
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Two models have been proposed to explain the observed linear relationship. 

The first proposed by Birch et al. (1968) has the heat producing Isotopes 
concentrated uniformly in a surface layer of finite thickness D. This 
thickness is given by the slope of the linear fit to the data. For the 
eastern United States the average thickness determined is 8 to 10 km. The 
second alternative model was proposed by Lachenbruch (1970). It assumes chat 
the abundance of heat producing isotopes decreases exponentially with depth 
in the crust with the logarithimlc decrement D. Lachenbruch (1970) has also 
pointed out that a linear decrease of isotopes would satisfy the relation. 
Lachenbruch further shows that, of the three models, only the exponentially 
decreasing source with depth remains valid under differencial erosion, which 
is a highly likely circumstance on earth. 

It may be overly optimistic to apply the relation to two points for the 

moon; nonetheless, this is done in Figure VII-A-1 and is shown by the solid 

-2 

line. The value of q* is seen to be 0.9uW cm and slope D is 19 km. The 
-2 

value of 0.9 uW cm for the heat flow into the base of the crust is within 
the range of values calculated by another analysis given in Section VII-B. 

There is already evidence that the value 0 cannot be applied to the lunar 
crust In a manner chat fits well for the earth's crust. Firstly, Metzger 
et al. (1971) have noted the very strong inverse correlation between surface 
thorium abundance and elevation of the lunar surface. The very lowest abun- 
dances are observed over the highest farside regions. Unless there are 
enormous lateral variations in the crustal abundance of radioisotopes, this 
correlation suggests that the heat production increases with depth, at least 
to depths of a few kilometers. This is exactly opposite to what is observed 
on earth. Secondly, the surface abundance of thorium seems to correlate most 
strongly with the extent of mare forming volcanism on the surface and inversely 
with the age of the mare volcanism. 

It seems a reasonable conjecture that the surface thorium abundance cor- 
related with intrusive type volcanism in the lunar crust as well. It has been 
suggested by several investigators that the upward extent of magma flooding 
or injection was controlled by internal stress (essentially hydrostatic) and 
equipotentlal surfaces at the time of injection. This model would fit the ob- 
served inverse correlation of elevation and surface abundance of thorium. Lower 
elevations simply have more volcanic intruslves than the higher areas which 
extended above the reach of magmas at the time of their emplacement. 




Heat flow vs thorium abundance; Che dashed line shows the expected 
heat-flow vs. thorium concentration if thorium is uniformly distri- 
buted through a 60 km crust. The surface thorium abundance in the 
regions Archimedes and Fra Mauro defined by Metzger et al. (1974) 
are indicated by arrows for reference. The stippled zone and line 
connect the observed heat-flow values (squares). 


Figure VTI-A-1 



Th« foregoing discussion is highly speculative; yet it is set down here 
to show the enormous and profound results that could be gained from the global 
mapping of thorium abundance, elevation and heat flow on the moon. They would 
provide data relevant to the chemical stratigraphy of the crust on a broad 
scale, the lateral heterogeneity of the crust, the extent both laterally and 
vertically of volcanisa subsequent to lunar crust formation. Such measurements 
have been proposed as part of a polar orbiter mission which was until recently 
an active candidate mission in the foreseable future. I 'fortunately, due to the 
costs of other activities ~y NASA, among them the Space Shuttle, such a scien- 
tifically valuable mission has been posponed indefinitely. 
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VII-B. TEMPERATURES IN THE UPPER 300 KM OF THE MOON 

The contents of this section are largely extracted froa a paper given 
at the Lunar Science Conference no. VIII. 

An abundance of the geophysical and geocheaical data froa Apollo is 
relevant to the present theraal regime of the lunar interior with varying 
degrees of significance. The selsalc data strongly imply the absence of 
any widespread melting in the outer 800-1000 km and a high Q upper aantle that 
extends to 300-S00 km (Nakamura at al. , 1976; Dainty et al., 1976). The 
observed anlsostasy of the circular mare oascons leads to rigidity requirements 
for the crust and upper aantle that have been interpreted in terns of viscosity 
and lithospheric thickness (e.g. Arkani-Hamed, 1973). The magnetics data, on 
the other hand, suggest that temperatures within a few hundred degrees of 
solidus may occur at depths as shallow as 200-300 Von (Kuckes, 1974; Sonett 
and Duba, 1975; Dyal et al., 1976). However, quantitative interpretation of 
the seismic and electrical properties, as well as the gravity data, in terms of 
temperature, are subject to the uncertainties of modelling assumptions, compo- 
sition, and laboratory measurements relating elastic, viscous and electrical 
properties to temperature. 

Data relevant to the total abundance and radial distribution of heat 
sources are provided by the Apollo 15 and 17 heat-flow measurements and the 
gamma-ray mapping of surface thorium abundances over the Apollo 15 and 16 
groundtracks . These >bservations are not subject to prior modeling or compo- 
sitional assumptions but do suffer the uncertainties of extrapolating measure- 
ments of limited areal coverage to global averages. We thus have available 
observations relevant on essentially a global scale (seismic, gravity and 
magnetics aata) which are subject to Interpretational uncertainties, and 
measurements of limited areal coverage (crustal thorium abundance and heat 
flow) which may or may not be globallv representative. It can be quite ar* 
bitrary which observations are considered 'constraints' on the models and which 
observations are to be tested by the model results. However, since the thermal 
models themselves require bulk and radial distribution of heat sources as 
Inputs, the Inferred crustaJ thorium enrichment and measured heat flew will be 
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referred to as constraints in the present work. It should be kept in nind, 
however, that this choice is only one of convenience. None of the observations 
can be considered definitive. The essential question is: Can a model of the 

present day thermal regime of the lunar interior be constructed which is consistent 
with all of the relevant observations? 

The one constraint which we do assume Is that the moon is presently thermally 
at steady state, in the sense that the heat generated Internally 1? balanced by 
the surface loss to space. Evidence in this regard is bared on theoretical as 
well as observational considerations. The absence of surface features related 
to thermal expansion or contraction implies that the mean lunar temperature has 
not changed considerably since the tine of mare emplacement. Present day high 
temperatures suggested by the magnetics data at 200-300 km and the subsolidus 
requirement imposed by the seismic data to i 00 km imply that radial temperature 
gradients in the interior are small. This, in turn, implies either that signi- 
ficant heat sources are absent from the deep interior or that heat now generated 
is efficiently transferred to the outer layers by a process of *olid state con- 
vection. In either case, a near steady-state thermal regime v„uid result. From 
a theoretical standpoint, the measured surface heat flow implies bulk heat 
source abundances one to two orders of magnitude larger than that required for 
the occurrence of convective transport, assuming a realistic temperature depen- 
dence of the material rheology. With the initiation of t solid state convective 
process, the time required to reach steady state would be greatly reduced, to a 
small fraction of the lunar age. 

In this section, we examine the coupling of global heat flow, crustal heat 
source enrichment, thermal conductivity and crust and upper mantle temperatures 
for the outer layers of a steady-state moon where conductive heat transfer 
dominates. By utilizing a range of heat source distributions chat would result 
from the formation of the crust by differentiation of the upper mantle and a 
thermal conductivity profile based on experimental data, comparisons can be 
made with model dependent constraints Imposed by other geophysical and geo- 
chemical data. 

O bservation s Relevant to Total Heat Source Abu nda nc e_and Radial_Distribution_ 

Heat flow -- For a moon near steady-state, the global heat flux can be 
interpreted directly in terms of the total internal heat generation. In situ 

' measurements carried out at the Hadley Rillc (Apollo 15) and Taurus Littrow 

2 

(Apollo 17) sites yielded regional heat flows of 2.1 and 1.4 uW/cm respectively. 


II 
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Earth based microwave observations of the lunar thermal emission (Troitsky, 1965), 
interpreted using laboratory measurements of regolith samples' electrical absorp- 
tion length (Gold et al., 1973), lead to subsurface thermal gradients near the 
center of the lunar disk consist-ut with the in situ gradient measurements of 
1.3 - 1.75 ®K/m (Keihm and Langseth, 1975) Combined with evidence for large 
scale homogeneity of the upper few meters of regolith, the microwave data 
support the global heat flow estimate. Based on these results, we have examined 

the implications constraining the range of average surface heat flow to valuer. 

2 

between 1.4 and 2.1 uW/cm . This constraint restricts the total heat source 

*8 3 

abundance of the lunar interior to values between 2.4 and 3.6 x 10 uW/cm 

With the total internal heat production constrained by the range of mea- 
sured heat flow, the present day radial distribution of heat sources can be 
estimated from geochemical and geologic evidence relevant to the total amount 
of differentiation which has taken place. Thermal history calculations have 
long pointed to a need for very high, possibly liquidus , temperatures in the 
outer layers of the moon early in its history to explain the episodes of mare 
volcanism and differentiation of the lunar crust (e.g.. ToksBz et al., 1972). 

Recent geochemical modeling studies (Walker et al.,1975; Hubbard and Minear, 

1975) suggest chat the early high temperature layer was on the order of 200-500 
km thick. A thermal analysis, based on the lack of extenslonal or compressional 
features on the surface of the moon, suggests that the early high temperature 
layer was limited to a thickness between 100 and 300 km (Solomon and Chaiken, 

1976) . These resul ? i_piy that the lunar crust was differentiated from a sub- 
stantial outer layer (the magma ocean), which would have a profound effect on 
the radial distribution of heat sources. A direct constraint on the amount of 
heat source differentiation which has taken place is provided by estimates of the 
mean abundance of heat generating isotopes concentrated in the lunar crust. 

The orbital gamma-ray determinations of the surface thorium abundances 
along the Apollo 15 and 16 groundtracks provide the most relevant data. These 
data have been shown to be consistent with the measured abundances of returned 
lunar samples from the Apollo sites (Metzger al. ,1974) . Bielefeld et al. (1976) 
recently updated the surface thorium data. Abundances in regions of high albedo 
(nonmare areas) are shown in Figure VII-B-1 as a histogram of values within 
5* x 5° areas. 
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Figure VII-B -1 j 


Histogram of the surface thorium abundance distribution over 5° x 5° areas of the Apollo l 1 ) and 16 
ground tracks from the analysis of Bielefeld et al. (1976) of the gamma-ray data. Only high albedo 
(nomuare) areas are included in the compilation 
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Measurements from mare have been excluded since observations from these 
areas are probably not representative of crustal abundances because the basalts 
that fill the mare basins are products of partial melting in the mantle or deep 
crust of the moon. Uranium and thorium would be preferentially concentrated into 
the liquidus fraction of a partially molten mantle. The most commonly found value, 
0.6 to 0.7 ppm, is observed over the farside highlands; the second peak at 1.0 
corresponds to limb highland areas. These values are highly biased toward 
Apollo IS observations because of the larger area covered. How to relate the 
surface distribution to the average crustal abundance is far from clear. A 
strong inverse correlation has been observed between surface thorium abundance 
and surface elevation (Trombka et al., 1974). This could be Interpreted as 
indicating that the crust is stratified relative to thorium in which case the 
values from the elevated farside regions would not be representative. On the 
other hand, it would be surprising if stratification in the highland crust could 
have survived the large scale mixing chat would have accompanied the intense 
bombardment by large bodies, unless the increase of thorium with depth was em- 
placed subsequent to 4 b.y.b.p. by widescale Intrusions into the lover part of 
the crust. Subsequent intrusion would not explain the relatively high values 
observed in the nearside terrae; therefore, the possibility of significant 
lateral heterogeneity of the crust cannot be dismissed. Until we have a more 
credible model of the crustal composition and stratification, an area weighted 
average of the surface abundances in nonmare regions offers the most unbiased 
estimate of the thorium abundance. As shown in Figure VII-B-1, this is 0.9 ppm 
for all regions (see also Taylor, 1977). Because of the uncertainties discuss, id 
above, we examine the consequences of constraining the average crustal thorium 
abundance to values between Che peaks of the frequency distribution CM). 63 - 
1.0 ppm). The lower limit admits the possibility chat the enrichment in Incom- 
patible elements may decrease with depth as it does in continental crust on earth 
and the upper limit reflects the uncertainties associated with the nearside data. 

It will be shown in a later section that the inferred range of crustal thorium 
abundance, combined with the heat-flow constraint, has important implications for 
the amount of differentiation which has taken place and present day temperatures 
in the lunar lithosphere. 
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The orbital data and sample data have established with a fair degree of 
certainty the isotopic ratios of potassium to thorium and thorium to uranium. 

We have used a K/U of 2500 and a Th/U of 3.7 to calculate heat production in 
terms of the thorium abundance . Because potassium is so volatile, the K/U 
ratio inferred for surface rocks may not be representative of the lunar interior. 
However, since present day radioactive heat generation is dominated by the 
uranium and thorium isotopes, interior K/U ratios even as high as terrestrial 
values would not significantly affect the model results. 

Description of Steady-State Models 

Thermal conductivity profile: — For a given heat flow and radial heat 
source distribution, temperatures in the conducting lunar lithosphere depend only 
on the thermal conductivity and its variation with depth and temperature. Rea- 
listic models of lithospheric temperatures must take into account the insulating 
properties of the lunar crust (relative to an olivine-pyroxene upper mantle) . 
Simplified models which assume a mean conductivity for the lithosphere appro- 
priate to an olivine-pyroxene mantle (0.04 W/cm *K is a commonly used value) 
will lead to temperature estimates for the upper mantle hundreds of degrees 
lower than when the insulating properties of the crust are taken into account. 

Laboratory measurements of lunar crustal rock samples show that porosity 
is the most important factor controlling the thermal conductivity (Figure VII- 
B-2, top). Mainly, the porosity depends on the extent of fracture of the original 
glassy or crystalline matrix, and to a lesser extent on the nature and composi- 
tion of the matrix. This suggests that a measure of the extent of fracture and 
crack closure with depth, such as the crustal p-wave velocity profile, can 
provide a good estimate of the thermal conductivity variation in the crust. A 
semiquantitative correlation is shown at the bottom of Figure VII-B-2. The large 
increase in V^ from just below the unconsolidated regollth to a depth of 20 km 
is believed to result from the closing of fractures as the compressive load 
increases. A porosity of about 5%, corresponding to a conductivity of 0.01 
W/cm°K, is believed to be an upper limit for a compressive load greater than 
1 km. Conductivities in the range of 0.018 - 0.020 W/cm®K at depths of 20-60 
km are based on the assumption that cracks and fractures are almost totally 
closed in the lower crust, as suggested by the p-wave data. Comparable crustal 
conductivities were obtained by Mizutani and Newbigging (1973) using a similar 
analysis. Crustal conductivity models + 30% removed from our nominal profile 
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were also examined co decarmine Che resultant uncertainties in the derived 
temperatures . 

Below the crust, the most important variations in thermal conductivity 
are probably due co temperature. Seismic and geochemical models suggest an 
olivine-pyroxene composition for the upper mantle with possibly a slight in- 
crease in Fe/Fe + Mg ratio with depth. The data shown in Figure 7II-B-3 are 
from the laboratory measurements of Schatz and Simmons (1973) on lunar mantle 
simulant materials. The samples show a slight decrease in conductivity to 
temperatures of about 1200°K, followed by a conductivity increase as radiative 
effects become important at higher temperatures. For cur nominal mantle con- 
ductivity model, we have adopted a theoretical conductivity versus temperature 
relation proposed by Schatz and Simmons for a mantle of FO^q (solid curve of 
Figure VII-V-3) . We have also examined mantle conductivity models + 0.01 W/cm 
®K removed from the nominal curve, thus bracketing nearly all of the experimental 
data in the temperature range of interest. 

Heat sources distribution: — A wide range of radial heat source distributions 
were examined by employing two parameter models of the early lunar differentiation. 
A given model was characterized by the depth, d, to which heat source differen- 
tiation took place and the fraction, a, of the initially uniform distribution 
of heat sources which were redistributed into the crust. For example, a model 
with d * 300 km and a * 1.0 has all heat sources originally distributed uniformly 
in the upper 500 km concentrated in the present day crust. For a given model, the 
total heat source abundance is determined by the assumed heat flow. 

For all models presented, an average crustal thickness of 60 km, based on 
the seismic velocity structure (Latham et ai., 1972; ToksBz et al., 1974), is 
assumed. The effects of different average crustal thicknesses are essentially 
equivalent to variations in the constraint of present day crustal thorium abun- 
dance. The abundance range we have assumed, 0.63 - 1.0 ppm, probably encom- 
passes the range of uncertainty of average crustal thickness. 

Each steady-state model was characterized by a choice of thermal conductivity, 
heat flow and differentiation depth and efficiency. Lithospheric temperature 
solutions for the radially symmetric heat conduction equation in spherical coor- 
dinates were obtained iteratively, assuming constant thermal conductivity over 
10 km Intervals. For all calculations, a fixed outer boundary temperature of 
310 # K was assumed. This is based on the porosity-conductivity relationship 
shown in Figure V1I-B-2, and the p-wave velocity structure inferred for the 
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1 ka at Che Taurus Littrow site (Cooper ec al., 1974). A flux boundary condition 

was calculated at the base of the lithosphere using the assumed global heat flow 

and radial heat source distribution. Models designated compatible were chose 

2 

with heat flow in the range of 1.4 >2.1 UW/cm and differentiation processes 
which resulted in present day crustal thorium abundances between 0.63 and 1.0 
ppm. 

Formulation of the Model: — Using the surface heat flow, F gur ^ the average 
volumetric heat production of the moon is calculated by the simple steady state 
relation for a sphere 


3 F J 
surf 


where H is the average volumetric heat production and is the radius of the 
moon taken as 1740 km. Petrological and chemical differentiation in the outer 
layers of the moon are assumed to produce three layers; the crust, depleted mantle 
and the primitive interior. The radial distance from the center of the moon to 
the upper boundary of each layer r^ , i ■ 1,3 could be varied in the model. The 
volumetric heat production in each layer was calculated using the following 
simple relations: 

For lunar crust: » H (1 + a) (r^ - r^ )/ (rj* - r^) 

For depleted mantle Hj ■ § (1 - 3) 

For the lunar interior * H 


Alpha is the fraction of heat producing isotopes depleted from the mantle. 

The thermal conductivity was defined within 10 km intervals based on 
the model shown in Figure VII-B-2. The initial conductivity in the mantle is 
given an assumed value which is uniform throughout. With the assumptions of 
thermal steady state, and spherical symmetry, the temperature at the base of 
any 10 km interval T^ can be computed from the temperature T^ at the top 
of the i - i layer by: 


W + V 38 !-, - 28 i-,- R i ,/6 K i-, R i- 
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where _ » Che heat flow at the top of layer l - i . 

equals the radial distance to the top of layer i, and 

h • Hi «« h > r 2 , 

H 2 for r 2 2 R i > r 3 » and 
H 3 for < r 3 , 

The flux into the base of layer i can be computed from the flux at the surface 
of the layer by 

f i - f i - , ■ r i - , /R i - « 3 i . , - < 

Using these formulae, the calculation is made layer by layer starting from the 
surface. 

The temperature dependence of conductivity in the mantle is assumed to 
have the form proposed by Schatz and Simmons (1973). 

K (T) « A - BT + DT 3 

This dependence is included in the model by recalculating for each mantle 
layer letting T » (T^ _ t + T^) /2 and repeating the computation of tempe- 
ratures until the sum of differences in temperatures of all of the layers was 
less than a degree on subsequent runs. 

R esults and Discussion 

For the nominal crustal and mantle thermal conductivity functions, the 
relationships between heat flow, differentiation depth, crustal thorium abun- 
dance and 300 km temperatures are shown in Figures VII-B-4 a, b, and c. The 
300 km temperatures are shown as dashed contour lines. The solid contours 
illustrate the relationship between crustal thorium enrichment and heat flow 
and differentiation thickness. The stippled region shows the range of tem- 
peratures at the base of a 300 km lithosphere which are compatible with the 
heat flow and crustal thorium observations as well as the subsolidus require- 
ment (less than 1600 a K) . Figures VII-B-4 a and c differ only in the assumed 
heat source differentiation efficiency, a. Note that for highly efficient 
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heat source differentiation ( a greater than 0.7S; Figures VII-B-4a and b) , 'magma 
ocean' depths of 300-800 km are implied, compatible with the geochemical require- 
ments. If early outer layer melting was restricted to depths of 100-300 km, as 
suggested by the thermal expansion constraint of Solomon and Chalken (1976), the 
present day average crustal thorium abundance must be near the lower limit 
(0.63 ppm) inferred from cite gamma- ray data. Implying that the heat source en- 
richment may decrease significantly in the lower crust. Much deeper differen- 
tiation layers can be compatible with the heat flow and crustal thorium cons- 
traints (Figure VII-B-4c) ; however, much less efficient upward migration of 
heat sources must be invoked. 

For a wide range of compatible heat flow and differentiation models, 

temperatures in the range of 1100 - 1600®K result at the 300 km depth. The 

requirement of a subsolidus upper mantle restricts the acceptable heat flows 

2 

to values less than 1.9 yW/cm . Observations of the magnetic field in orbit 
and on the lunar surface yield determinations of the electromagnetic response 
of the whole moon to transient variations in the ambient magnetic field. The 
response has been examined in the time frequency domain by several authors (e.g. 
Sonett and Duba, 1975; Dyal et al., 1976) by comparison with electrical conducti- 
vity models. All of these analyses require a rapid increase in electrical con- 
ductivity with depth in the outer 100 - ZOO km, and values on the order or 
-4 -3 

3 : 10 - 5 x 10 mhos/m are reached at 300 km. Uncertainties remain in the 

analysis of the magnetic data which stem from asymmetries in the enveloping 
electric field produced by the solar wind and the effects of higher modes of 
Induced response. A recent analysis by Dyal et al. (1976) of the transient 
response when the moon was shielded from the solar wind in the earth's magnetic 
tail, yielded results in general agreement with those of other workers. 

Correlation of inferred electrical conductivity and temperature is based on 
the strong temperature dependence demonstrated in the laboratory by common 
silicate rocks. Consideration of all data shows a relatively large scatter and 
led to a large uncertainty in early interpretations of mantle conductivities 
in terms of temperature. More recently, however, the role of oxygen iugacity 
has been more fully recognized and the measurements of Duba et al. (1974, 1976) 
on olivine and orthopyroxene samples under simulated lunar conditions have 
considerably narrowed the uncertainties in the temperature deductions. Mantle 
temperatures based on the electrical conductivity models of three different 
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Three hundred kilometer temperature contours plotted versus average 
surface flux (F s ) and differentiation depth. The stippled area des- 
cribes the range of 300 km temperatures and differentiation thick- 
nesses which satisfy the heat flow and crustal thorium abundance 
data and the subsolidus requirement. The three sets of contours 
are distinguished by the assumed differentiation efficiency, at: 

(a) a • 1.0; (b) a » 0.75; (c) a • 0.50. 


Figure VII-B-4 c 
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investigators are shown In Figure VII-B-5 along with a range of salenotherma 
produced by models compatible with the heat flow and crustal thorium abundance 
constraints. Note that all three electrical models require temperatures within 
200*K of the solidus (Hodges and Kushiro, 1974) at depths greater than 200 km. 
The representative selenotherms based on the present work are consistent with 
these results, although temperatures as low as 1100*K at 300 km can also be 
compatible with the heat flow and crustal thorium abundance constraints. The 
approach towards solidus temperatures at about 300 km is also supported by the 
seismic model of Nakamura et al. (1976). A large negative shear wave velocity 
gradient above 300 km suggests large temperature gradients In this region. The 
decrease In shear wave Q from 4000 to 1500 at about 300 km may indicate a sof- 
tening of the material as the solidus Is approached, causing the lower mantle to 
bshave more plastically to stress. It has been argued previously that high 
seismic Q values precluded temperatures within a few hundred degrees of solidus 
in the upper mantle. Especially relevant to this point are the high p-wave Q 
values of "v 5000 to depths of approximately 500 km inferred from the analysis 
of seismic surface events (Dainty et al., 1976). These results conflict some- 
what with the seismic models of Nakamura et al. (1976), which require a lower 
s-wave Q below 300 km. Based partly on the higher seismic Q values of Dainty 
et al., lower temperatures and significant temperature gradients were expected 
to depths of at least 500 km. However, the laboratory measurements of Tlttman 
et al. (1976, 1977) indicate that waterless lunar materials can exhibit very 
high Q values even at elevated temperatures. Recent measurements up to 2.5 
kbar and 800*K (Tlttman, personal communication) suggest that these effects may 
remain significant at pressures and temperatures characteristic of the lunar 
upper mantle. Seismic property measurements up to 15 kbar and 1200®K could be 
crucial in resolving the high temperature, high seismic Q requirements of the 
lunar crust and upper mantle. 

A second objection to high upper mantle temperatures has been the rigidity 
requirements imposed by the anisostasy of the circular mare mascons. Kuckes 
(1977) has suggested that the load of the mascons co'*ld be supported by flexure 
of a relatively thin elastic lithosphere over a viscous mantle. He finds that 
the moon's elastic lithosphere must flex as a plate about 60 km thick to support 
the mascons. On earth, the oceanic lithosphere flexural thickness has been 
shown to be about 24 km (Watts and Cochran, 1974). Temperatures of about 750*K 
are reached at this depth in the oceanic lithosphere greater than 50 m.y. old 



; i 177 

i i 



Figure VII-B-5 j‘ 


« f 




178 


(see, for example, Langseth and Von Herzen, 1970). Thus, comparable temperatures 
might be expected at 60 km in the sic , that is at the base of the lunar crust. 

It Is seen from Figure VII-B-5 that the compatible selenotherms require crustal 
base temperatures between 700 and 850°K. 

Because of the temperature dependence of the thermal conductivity below 

the crust, thermal gradients are limited to values between about 1.8 and 3.2 a K/ka 

for compatible models. These gradients ara not large enough to account for the 

measured shear wave gradient of - 0.0013 km/sec/km (Nakamura et al., 1976), 

without invoking an additional effect, such as an increase in the Fe/Fe + Mg 

ratio with depth. Experimental data by Chung (1971) gives 5 V g /5T ■ -3.2 x 

10 Sm/sec/*?., <$V /6p * 3.7 x 10~^km/sec/kbar; and Sv /<5% (Fe/Fe + Mg) * 

s s 

• 0.015 km/ sec/ X (Fe/Fe + Mg) for forsterite samples. Combining these values 
with the observed shear wave gradient of - 0.0013 km/ sec/km yields a thermal 
gradient of 4.6 °K/km if no compositional changes are significant. Thus, if 
the actual gradient is aoout 2.5 °K/km, an increase in the Fe/ (Fe + Mg) ratio 
of 0.032 km may be implied for the upper mantle. 

The discussion to this point has been based on results obtained using 
the nominal thermal conductivity profile described earlier under the 'Des- 
cription of Steady -State Models'. The effects of the conductivity uncertainty 
are best illustrated by assuming extreme cases of uniformly higher and lower 
crustal and mantle conductivity profiles and examining the effects on the thermal 
regime still subject to the heat flow and crustal f jrium abundance constraints- 
The selenotbr ms shown in Figure VII-B-6 are for variations in crustal conduc- 
tivity of + 30Z from nominal and variations in mantle conductivity of + 0.01- W/cm 
°K, which nearly brackets the laboratory data for olivine and pyroxene samples. 

The model parameters are d * 500, a * 1.0, F *1.5 and Th _ * 0.91 ppm for 

s crust 

all five cases. Variations in crustal base and 300 km mantle temperatures of 
+ 150*K result. The effects of the crustal and mantle conductivity uncertainties 
are approximately additive so that errors of up to + 300°K at 300 km could result 
if both the mantle and crustal conductivities were uniformly in error by the 
amount of their stated uncertainties. 

Conclusions 

Pv utilizing the steady-state criterion and a thermal conductivity model 
for the crust and upper mantle based on experimental data, the relationship 
between global heat flow, crustal heat source enrichment and lithospheric 
temperatures can be quantitatively evaluated. If the measured ranges of crustal 
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thorium enrichment and heat flow are globally representative, the features of 
a lunar thermal lithosphere that extends to at least 300 km can be summarized 
as follows: 

1) At 300 km, temperatures are In the range 1100 - 1600°K. Magnetic models 

of the upper mantle electrical conductivity and the requirement of a subsolidus 

outer regime favor temperatures in the range of 1300 - 1500 *K and a global heat 

2 

flow of 1.4 - 1.8 yW/cm . Three hundred kilometer temperatures significantly 

2 

lower than 1200*K require a global heat flow less than 1.4 yW/cm or a crustal 
heat source enrichment greater by a factor of two than that Inferred from the 
gamma- ray data. 

2) Models of h- it source differentiation require an outer layer extending 
to depths of 300 - 600 km if 100Z of the heat sources in the differentiated zone 
have migrated upward into the crust. If 751 have migrated upward, the layer 
could be as much as 300 km thick. The 100Z differentiation model compares 
favorably with the 'magma ocean' model thlckne -s of 200-500 km derived from 
broad geochemical constraints and 100 - 300 km based on Solomon and Chaiken's 
(1976) analyses of constraints on the allowable thermal expansion and contraction 
of the moon compatible with the absence of related surface features. Alternate 
models of deeper heat source differentiation can be compatible with the heat 
flow and crustal thorium constraints only if a significant fraction of the 

heat sources remain in the deep interior. 

3) Upper mantle gradients in the range of 1.8 - 3.2 “K/km are required at 

200 - 300 km depths. These values are not high enough to account for the observed 

seismic shear wave gradient and suggest that an additional effect, such as a 

small increase in iron abundance with depth. Is required. The mantle gradients 

imply significant heat generation below 300 km. Heat fluxes in the range of 
2 

0.7 - 1.3 yW/cm at 300 km indicate that 1/3 to 1/2 of the measured surface heat 
flow is derived from deep interior heat sources. 4n efficient solid state 
convection process is required to maintain subsoil 'us temperatures to depths of 
800 - 1000 km. 

The primary uncertainty for the steady state models is the thermal conduc- 
tivity and its variation with depth and temperature in the outer, nonconvecting 
layers of the moon. A realistic range of possible conductivity error leads to 
uncertainties of 100 - 200®K in the derived upper mantle temperatures for a given 
model. The range of compatible heat source differentiation models is not strongly 
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affected by the conductivity uncertainty. Neither is the general conclusion from 
this analysis that the measured heat flow and inferred crustal thorium abundances 
are consistent with the high upper mantle temperatures Interpreted from the 
magnetics and seismic data, a differentiated outer moon supported by broad 
geochemical constraints and a convecting subsolidus lunar Interior. 
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VII-C. IMPLICATIONS FOR THE MOON'S BULK CHEMISTRY 


As we have stated in other sections of this report, there is 
strong evidence that at present the moon is close to thermal steady state. 

Solomon and Chaiken (1976) have drawn attention to the lack of compressional 
or extensional surface tectonic features on the moon, indicating radially 
averaged internal temperatures have been nearly constant since the time of the 
last mare flooding, about 3.2 b.y.b.p. The extremely low level of mean surface 
seismic activity (Latham et al., 1973) indicates that thermally induced stress 
levels are very small at present. Thermal history models of the moon containing 
radioisotopes at concentrations sufficient to explain the surface heat flow and 
which fit other well established constraints (see, for example, ToksBz and 
Johnston, 1974; ToksBz et al., 1976; and, ToksBz et al., 1977) indicate that 
most of the radioactive heat sources have been concentrated in the outer 100 km 

9 

where their generated heat would equilibrate with the surface flux in 10 years 
or less, or the deep interior is convecting so that the heat transfer is enhanced 
and steady state conditions are more rapidly approached. 

For a moon at steady state, the heat generation In the interior is 
related to the average surface heat flow by 

h -ai 

v r 

where F is the mean heat flow and r is the lunar radius. If the thorium to 
uranium ratio is assumed to be 3.7 and the potassium to uranium ratio for the 

3 

whole moon is 2000, then the volumetric heat production in W/cm is related to the 


uranium abundance [U] in gm/gm by 


H - 6.73 x 10 [U] 

v 


From these two simple relations we can estimate the global heat flux and the 

uranium concentration in ppb. In Section VII-B we showeo chat for a moon at 

2 

steady state hear, flow values in the range of 1.4 to 1.8 uW cm are consistent 
with other data relevant to interior temperatures in the moon. 


i 


-s 6o 
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If an analysis similar to that described In Langseth et al. (1976) Is 

carried out, a similar range of values is found. The analysis used the two 

measured values of heat flow, a mean heat production in the crust, H and 

c 

determinations of crustal thickness along the Apollo 15 and 16 tracks to cal* 
culate a mean heat flow for these areas. Isostatlc adjustment of the crust 
is assumed so that lunar elevations can be converted into crustal thicknesses 
if a density contrast Ap, between the crust and mantle is assumed. If we start 
from the seismically oj.ermined crustal thickness of 60 km over eastern Procel- 
lanum and use the mean densities assumed by Solomon (1975) , a crustal density 
of 2.8 and a mantle density of 3.4, Ap - 0.06 and the isostatic assumption, then 
the thickness T over highland areas and nonmascon mare plains is given by: 

T - 73 + 5.7 h 

where h is the elevation relative to the mean lunar radius of 1738 km. In areas 
of mascons the surface mass concentration of 10^ gm/cm is assumed to be uncom- 
pensated and, 

T - 56 + 5.7 h. 

In nonmascon areas that are not enriched in surface thorium, the heat flow 
within an area was calculated by the relation 

F « F u + T • H 
M c 

where F u is the flux from the mantle (the value at Taurus Littrow) . Values 
determined in Section VII-B can be used here. In areas where the surface 
thorium is enriched 


F » F„ + T 
M 


H +F .. 
c add 


"»«• * D ' • (H o - 5 => 

H q Is the surface heat production calculated from the surface thorium abundance 
and D' is the slope of the line in Figure VII-A-1. 
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Once the values In each area are determined by this model, the average 
heat flow Is calculated by 


N 

■ Z 

u-o 




where F_ and A are Che heat flow and area of each area of available daca. 
a n 

Using che areas defined by Meczger ec al., (1974) or Bielefeld ec al. (1976) 

Che values ac Hadley Rille of 2.1 yW/cm 2 and Taurus Littrow of 1.4 yW/cm” J can be 
used co calculace a mean value. For values of H £ becween 0.9 x 10 " 13 W cm” 3 
Co 1.55 x 10‘ “ 13 W cm -3 mean heac flows of 1.4 co 1.8 yW/cm ” 2 are obcalned. 

This range of global heac flows would correspond co average maximum 
concencs of 36 co 46 ppb for che case of a moon ac sceady stace. More detailed 
calculaclons of che Chennai evoluclon of che moon from lnlclal condlclons and 
Including posslbllicles of conveccion lead co heac flows of abouc 1.7 yW cm” 2 
If a uranium concent of 35 ppb uranium la assumed (Toksflz ec al. t 1977). 

However, che compucaclons by ToksSz ec al. lead co substancially lower tem- 
peracures In Che lunar incerlor ac presenc chan our sceady scace models wi,th 
Internal conveccion implies. The model of ToksHz et al. Indicates chat tem- 
peratures of 1500 a K are reached at abouc 800 km compared with 300 km deduced 
from che sceady state model. The Tokstiz ec al. model also concludes Chat solid 
scace conveccion is confined to a relatively small region of che lunar interior 
deeper chan 1000 km. These models assume crustal conductivities chac are equal 
to Chose of che mantle. This assumption would greatly decrease gradients in 
Che crust and enhance che heat flux from che mantle to che surface during lunar 
history. The thermal conductivity is not known with sufficient certainty to 
settle on a 'best' model ac this time. 

An important implication of uranium abundances of 36 ppb is thac they are 
very similar co abundances now implied for che earth's mantle. 

Recently, estimates of che earth's average heac flow have increased Co 
account for che contribution from che seafloor spreading process. Williams and 
Von Herzen (1974) have given a new estimate of the global heat flow from the 
earth of 8.4 yW/cm 2 , which is 302 higher than the previous estimate of 6.3 
yW/cm 2 , based on an average of all heac flow observations. For an earth ac 
sceady scace, a uranium abundance of 32 ppb is required to produce a flux of 
8.4 yW/cm 2 if we assume che usually accepted K/U and Th/U ratios. If all of 
che heac sources are concentrated in che mantle then che abundance of uranium 
in che mancle is 42 ppb. 
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VIII. RECOMMENDATIONS & rdTURE MEASUREMENTS OF LUNAR HEAT FLOW 


The present status of our knowledge about the heat flow from the interior 
of the lunar surface consists of direct measurements of heat flow through 
ehe lunar surface layers at two widely separated locations in the northern 
part of the nearside, numerous measurements of the microwave brightness 
temperature versus wavelength using earth-based radio telescopes, and a 
large library of new data on the thermal and electrical properties of lunar 
surface materials that resulted from the lunar sample program. 

Indirect observations that set limits on heat flow at depths of a few 
hundred kilometers come from the study of seismic wave propagation in the 
coon and the whole moon's interaction with the time variable magnetic field 
of the solar wind as measured by surface and orbiting magnetometers. These 
observations indicate gradients of seismic velocity, attenuation coefficients, 
and electrical conductivity in the lunar interior which can be related through 
experimental data to temperatures and radial temperature gradients in the lunar 
mantle (Nakamura et al.,1976; Sonett and Duba, 1975). 

The two direct surface measurements of heat flow made during the Apollo 
program were made in similar physiographic settings; i.e., on the edge of 
two of the moon's major mascon basins (Langseth et al., 1974). How repre- 
sentative these values are of the true global heat flow from the moon is not 
known because the variability of flux over the surface is not known. 

The gamma-ray results taken from orbit indicate larger variations of 
thorium over the surface (Bielefeld et al., 1976) as well as variations in 
critical isotopic ratios. Far infrared emission spectra also give evidence 
for compositional variations over the lunar nearside. The orbiting x-ray 
and gamma-ray results show similar variations over the belts covered by the 
Apollo 15 and 16 Command Modules. If these compositional variations continue 
to considerable depth in the moon, they would result in variations in surface 
heat flow. 

Near-surface measurements of heat flow are also subject to near-scrfaca 
disturbances. One of the principal disturbances for which we have little 
data to estimate the effect is the refraction of heat flow at sloping inter- 
faces between bodies of different thermal properties. If, for example, dense 
rock rises close to the surface bereath rhe mountainous rings of the circular 
mare basins, they could have a considerably higher conductivity than the 
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overlying regolith, and consequently, a pronounced effect on the distribution 
of heat flow near such structures (see also, Conel and Morten, 1975; Mitzutani 
and Osako, 1974). Consequently, it Is extremely important to obtain further 
data that may determine the variability of heat flow over the lunar surface. 

Microwave Spectra of the Moon ; 

One class of observation that bears directly on the variability of heat 
flow over the lunar surface is that of the microwave spectra between the wave- 
lengths of J and 100 cm. The dielectric properties of the moon's surface are 
such that it is semitransparent tc se wavelengths and, thus, the intensity 
of any given wavelength is determined by the integrated temperature profile 
to depths roughly equal to absorption length of the given wavelength. Con- 
sequently, microwave spectra can be interpreted in terms of the vertical 
temperature gradient if the dielectric properties, in particular the loss 
tangent is known. The first effort to determine the heat flow from the moon 
using microwave spectra was made by B^' in (1961) based on a vary long 
wavelength observations (178 cm). At a>out the same time, similar obser- 
vations were being made by Soviet astronomers. One of their publications 

-6 -2 

(Troitsky and Tikhonova, 1970) inferred heat flows 2.9 to 4.0 x 10 cm . 

Recent improvements in radlc astronomy have greatly impreveu the resolution 
cf modem telescopes. Conceivably, earth-based telescopes could map the 
lunar nearside in the wavelength rar ge from 3-20 cm at sufficiertly gcod 
resolutions to detect regional variations in the brightness temperature 
gradient. Measurements of this type are currently being carried out by 
scientists of the Jet Propulsion Laboratory. 

Observations of the variation of the microwave brightness temperature 
spectra with solar phase angle provide a direct of the ratio 

between the thermal and electrical absorption len^i.-i vs, the microwave 
data contain information on the variability of the 9 iperties over 

the surface. However, an indep indent estimate of either iacmal diffu- 

sivity or the electrical absorption length is require • .^it microwave 

brightness temperature gradients to temperature gradi? <• ;r heat flow. 

These estimates must be based on laboratory measurements of Apcllo samples 
and in situ measurements of thermal parameters. The. uniformity of conducti- 
vity values deduced from the heat-flow experiment are particularly significant 
In this regard. 
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Microwave Measurements from Orbit 

It has long been recognized char passive microwave observations from a 
lunar orblcer would map heat flow over the entire lunar surface if, as in 
the case of earth-bound observations, the thermal or dielectric properties 
can be Independently estimated. Such aeasureme :s were considered by the 
1965 Falmouth meeting described earlier. More recently, NASA has been con- 
sidering a Lunar Polar Orbiter mission as a follow-up to the Apollo Program. 

If this mission carried passive microwave radiometers, it could map the 
brightness temperature gradient and the parameter, b, which is the ratio of 
microwave to thermal extinction length over the lunar surface. Such an 
instrument could map the variability of heat flow at the scale of 100 km 
or less. 

An instrument to be flown on a polar orbiter has been proposed by a 
group at the California Institute of Tacimolog/. Work is now in progress 
on the problem of interpreting these types of data. Maps of global heat 
flow over the entire lunar surface, whe**-. combined with simultaneous mea- 
surements of the surface composition b> gamma-ray spectrometry and other 
techniques, premises a synoptic view of the moon’s surface and bulk chemistry 
which would rival our knowledge of the earth. 

Future In Situ Measurements of Lunar Heat F low 

Remote observations relevant to lunar heat flow will always suffer from 
the problem of o\';iag a unique interpretation from the data because of uncer- 
tainty in estimates of certain crucial in situ physical properties. Thus, 
additional in situ measurements of the thermal properties of the lunar surface 
layer and the heat flow through the lunar surface would be extremely valuable 
because it would allow a refined interpretation of the microwave results*. 
Especially important is an improvement in the absolute value of the heat flow. 

Based on our present knowledge of the lunar surface 1 'yer and its thermal 
regime, accurate surface heat-flow measurements could be made from a soft 
landing station of moderate size and cost. 

R ecommendations for the Design of Futur e | Inst ruments to Mea sure t he Lunar 
S urface Heat Flow 


The merhanieal and thermal properties of the lunar surface layer were 
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reviewed la Section III-A. Those properties that most strorgly Influence 
future heat-flow Instrument designs are reviewed hern. 

1) The regolith is covered nearly everywhere by a thin 2-15 cm layer of 
relatively low density, fine grained ro r* and glassy powder that ha? extre- 
mely low thermal conductivities on the order of 2 - 8 x 10 * watt/cm *K. This 
layer has two important effects on the time-varying thermal structure of the 
regolith. Firstly, radiative transfer is important and a steady state tempe- 
rature difference of several tens of degrees will be established across this 
layer. Secondly, the very low thermal conductivity will greatly reduce the 
periodic variations of temperature due to diurnal fluctuations of the solar 
flux at the surface; see Section V-A where the depth dependency of the am- 
plitude and phase of diurnal variations is shown. The effect on the mean 
temperature of the -rad' tive transfer in the uppermost layer based on data 

at the Apollo 15 and 17 sites are also shown. 

2) Below th low density surface layer, the reg-^ith fines are densely 
packed, approaching th maximum possible density (Carrier, 1974). The con- 
ductivity at these depths is 4-7 times greater than the near surface material. 

-4 -1 

Conductivities of 1 to 1.5 x 10 watts cm*K were deduced from the attenua- 
tion of the annual wave with depth (Langseth et al., 1976). Ine regolith exhi- 
bits some stratification at the millimeter and centimeter scale in the lunar 
regellth, baseu on the grain size distribution -’nd compcsitlon. Consequently, 
seme variability of the thermal conductivity can be expected. Despite the 
density of subsurface material radiative heat transfer may still be impert?’ “ 
and any thermal axperlments such as those needed to measure the thermal con- 
ductivity should produce temperature differences less than three degrees. 

3) Almost any disruption of the lunar surface will result in a change 
in the surface radiative end thermal properties and, therefore, the mean 
surface temperature. This change will produce a transient disturbance that 
propagates slowly downward from the surface. See Section V-C which shows 
examples of the disturbance induced at various depths by a 1 # C .hange in 
surface Temperature over a one meter radius circle. To avoid disturbance of 
any measurements by : v .-se effects, the measurement period should either be 
constrained to a short period after installation; or a very long period cf 
one year or more so that transients can be observed and eliminated from the 
data. 
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4) The measured heat-flow values at Hadley Rille and Taurus Llttrow 

2 

of 1.4 and 2.1 watt/cm induce vertical temperature gradients on the order of 
a degree per meter in the subsurface to depths of about 2-3 meters. 

5) Heat flow may vary laterally over the lunar surface at the scale of 
a few meters due principally to variations in the radiative properties of 
the very near surface layer and the lateral heterogerei ty cf subsurface 
thermal properties. Thus, within any given region numerous measurements 
should be made to measure and average out this variability. 

6) The regolith layer is relatively easy to penetrate with a tool of 
small ross sect .<"<al area and a sharp hard tip as was amply demonstrated by 
numerous core tubes driven intc the soil and penetrator tests carried out 

on the .unar surface by the astronauts. 

The characteristics descr-'bed above apply over a very large proportion of 
the lunar surface. There are some isolated regions, however, where hard rock 
and higher thermal conductivities lie only a few centimeters below the sur- 
face. The IR hot spots (see e.g. Shorthill, 1972) identified from earth- 
based telescopes and lunar orbit (Hendell and Low, (1975) are one such area . It is 
also possible that on slopes greater than 10° the regolith may be considerably 
thinner than over more level areas. 

Concep tual Design fo r a Lunar Station Heat Flow. E xper iment 

The thermal and mechanical properties of the lunar soil layer described 
above are such that an accurate measurement of the mean temperature gradient 
(related to the steady outward flow of heat) could be obtained by a probe 
inserted only one meter into the lunar surface. The gradient measurement 
would be made in the lower 70 cm. 

Figure VIII-1 shows a conceptual design for such an experiment. The 
main component is the long needle probe that makes the temperature and con- 
ductivity measurements. This probe could be made cf titanium alloy tubing 
with a sharpened carbide steel tip. If a relatively high conductivity material 
such as titanium is used, the needle probe must have a large length to dia- 
meter ratio. In the design, we recommend a ratio of 200, which is more chan 
adequate. 

The temperature profile would be measured by a single small cemperature 
element near the tip of the needle (the prime cemperature element). Use of 
a 3ingle element would reduce the stringent requir -nents on the accuracy of 
an array of elements such as those used in the Apollo instrument. Thus, a 
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sensor such as a thermistor could be used. 

The temperature profile would be measured by driving the probe, which is 
initially retracted inside the housing, slowly into the soil. A motor driven 
lead screw and vibrating head at the top of the needle probe accomplish this. 
The probe makes ten or more predetermined stops on the way down. A wait at 
each stop is made long enough to allow calculation of the equilibration tem- 
perature of the adjacent regolith material. For a probe 0.5 cm in diameter, 
made of titanium alloy, the equilibration time is about a half hour. 

Additional sensors along the probe (K-sensors) would be used in con- 
junction with a single filament wire running the length of the probe to mea- 
sure the thermal conductivity. The technique for conductivity measurement is 
the well known needle probe technique used in seafloor geothermal studies. The 
absolute accuracy of these sensors is not important. A resolution of about 
0.003 K is required for the measurement and stability over a six hour period. 
The conductivity measurements require about six hours. Measurements at all 
locations would be made simultaneously. 

Measuring the gradient with the 'moving sensor' technique can only be 
used below depths of about’ 30 cm, because of the sizable variations intro- 
duced by the diurnal wave in the near surface regolith and heat transfer 
along the probe from the housing on the surface. 

The instrument should also include surface temperature thermometers. 

Lunar surface thermometers are hard to design. The primary temperature 
element sees thermal radiation, and during the lunar day reflected solar 
radiation, from the lunar surface over a field of view of 2 * steradians. 

Over the other 2iT steradians, the sensor sees the interior of blackened 
copper cavity, which is kept at as low a temperature as possible. The 
average interior temperature of the cavity is measured by additional 
thermometers . The cavity is kept relatively cold by shielding it from direct 
sunlight and giving it the largest possible view of outer space. 

The greatest uncertainty arises from the solar albedo which will be a 
factor only when the sun is up. Thermometers with an accuracy of a few 
tenths of a degree would be adequate to give temperatures with far greater 
accuracy than heretofore. The surface temperature observation is required to 
determine the forcing temperature function at the surface. 

Surface temperature observations at a single point should be taken 
during one lunar day. This would allow determination of the variation of 


1 


| 194 

i 

surface temperature over a complete lunation cycle. From these data, . • 

./, i 

the mean surface temperature and an estimate of the parameter (K pc) j 

of the upper 10-15 cm of material can be determined. ; 

Penetration of the regolith would be accomplished during either the ‘ • 

lunar day or xunar night, depending on when the most power was available. 

The subsurface measurements could be accomplished in approximately ten earth 
days. Temperature profiles measured by the heat-flow probe would have to be 
corrected for transient disturbance caused by the diurnal and annual varia- 
tion propagating downward from the surface. The instrument would obtain 
sufficient data to make these corrections quite accurately. 

This approach could be used on any planet with an extremely tenuous 
atmosphere and penetrable surface layer. 

Other Possible Techniques to Measure In Situ Heat Flow: 

The method of heat-flow measurement described above could only be carried 
on a soft landing spacecraft. Another type of mission that has been considered, 

! is the so-called penetrator mission. By this concept, small probes would be 

I 

| launched from an orbiting spacecraft and bury themselves in the lunar surface 

i - , 

S layer. The body of the probes and internal instruments would have to be of 

| rugged design to withstand impact and penetration. There is a major tech- 

nical problem to using the penetrator technique on an atmosphereless planet. 

A rather accurate orientation and guidance technique is required to assure 
vertical penetration, or the vehicle will be destroyed by shearing forces 
on penetration. 

i Such a penetrating probe would provide good penetration at most areas 

! on the moon; deep enough so that undisturbed gradients, due to heat flux 

! 

from the interior could be measured, after the initial heating during pene- 
tration has dissipated. The dissipation would require at least two to three 
lunations. A serious problem is the heat generated by any dissipating cir- 
cuits or power sources in the penetrometer. Because of the extremely low 
conductivity of the regolith material, very high temperatures would result 
from very low power levels. At the present time, the penetrator mission 
does not appear to offer an attractive possibility for accurate lunar heat- 
flow measurements. 

I 

! ' i 
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APPENDIX A. EARLY EXPERIMENTAL CONCEPTS 

* 

The test of this section Is largely extracted from a 1966 report by 
Arthur D. Little of a study of various concepts for measuring the temperature 
gradient la the regolith and the thermal conductivity. 

To guide the design of a probe to determine lunar heat flow, we (ADL) 
established a list of "desirable characteristics." The relative ranking of 
many alternative concepts and selection of the probe design was guided by 
these characteristics. The desirable characteristics are first outlined in 
general terms and then a number of specific designs and their advantages and 
disadvantages are discussed. 

DESIRABLE PROBE CHARACTERISTICS 
Fall Safe 

An ideal probe is "fail safe." That is, it is highly desirable that the 
probe should be capable of yielding useful data even if one component fails. 
For example, there is an increased possibility of failure of a probe with 
moving parts, however, by appropriate design such a probe would not become 
completely useless. 

In situ Calibration 

Both absolute and differential temperatures are to be measured. The 
heat flow experiment requires very high accuracy in the differential measure- 
ments. It is advantageous to determine, at least intermittently, whether 
the differential temperature signals are real or spurious. A simple check 
is to force the differential sensors to read the same temperature; a zero 


*Design Definition of Heat Flow Probe - Summary Report by Arthur D. Little, 
Inc. 68399 
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signal indicates that any transmitted signal is, in fact, due to a temperature 
difference. In the absence of an in situ calibration the sensor has to main- 
tain the original calibration during launch, traverse to the moon and obser- 
vation period. 

Thermal Conductivity Measurements 

Thermal conduct! /ity of the lunar material has to be measured in the 
region where the temperature gradient is to be determined. Depending upon 
the local conductivity two techniques are to be used. For low-conductivity 
materials, an integral heater-temperature sensor is necessary. For high- 
conductivity materials, the heater and temperature sensor must be separated 
axially. To determine axial inhomogeneities in the lunar material, thermal 
conductivity is to be measured at different axial locations. Thus, in any 
probe design at least one heater and two temperature sensors have to be pro- 
vided. One of the temperature sensors can be the same as that used to measure 
the temperature gradient. 

Maximum Interval for Temperature Difference Measurement 

It is desirable to measure the axial temperature difference over the 
total probe length. Additional measurements over a portion of this length 
or multiple probes are desirable if weight limitations are not exceeded. 
Simplicity 

The design of the probe should simplify the experiment. Certain sensors 
require simple electronics; others do not require motion of any part relative 
to another. Simplicity will increase reliability, require less development 
work, and thus shorten development time and cost. 

Emplacement 

The probe has to be designed to be compatible with astronaut human 
factors, emplacement tools, and packaging requirements. Both emplacement 
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at the required dopth and radial positioning within the hole have to be 
accomplished by an astronaut. 

Low Axial Conductivity 

To measure small temperature differences accurately in the lunar material, 
the probe must not behave as a thermal short in the axial direction. The axial 
thermal resistance of the probe vail, insulation, lead wires, etc., must be 
large. 

Weight 

The probe weight should be as small as possible. 

Resistance to Vibration and Shock 

Each probe must be designed to meet ALSEP environment specifications 
without decreasing the operational effectiveness of the probe. 

Conceptual Probe Designs: 

THE STATIC PROBE 

The static probe is the simplest type. Temperature sensors are placed 
at or near the probe ends and held apart by a lightweight, strong, thermally 
resistive structural member. A schematic diagram of this type of probe is 
shown in Figure 1. Two additional temperature sensors would be needed about 
8-12 cm from the end heater-sensor units to measure temperature when operat- 
ing the thermal-conductivity experiment for high-conductivity lunar materials. 

Additional temperature sensors can be incorporated within the probe to 
measure the gradients between different locations, or two probes can be 
connected in tandem to determine temperature differences over an additional 
distance in the drill hole. 

The principal disadvantage of the static probe is the dependence on 
the absolute sensor stability. There is no way to calibrate the probe in 
situ to test that measured temperature differences are real and not spurious. 


50 cm 
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Lead Wires 


LEGEND: 

□l Isothermal Section 
E3 Insulation Section 
S Sensor 

H Heater 

Sj£ Thermal Conductivity Sensor 


SCHEMATIC DIAGRAM OF STATIC PROBE 


Appendix A- Figure 1 
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This type of probe can only be used when Che long-term stability of the 
two end sensors can be ensured. However the thermal conductivity sensors 
on the probe wall do not need long-term stability since they are only used 
to measure a temperature rise during a short time Interval. 

DYNAMIC PROBES 

In a dynamic probe at least one moving part is incorporated within the 
probe. The dyanmlc probe concepts can be subdivided according to the required 
movement: long, short, or very short. Long-movement devices require relative 
movement of one part over the entire probe length, e.g., to be accomplished 
by a motor at the lunar surface. Short-movement devices require movement 
over a few centimeters, e.g., by a solenoid mounted on the probe. Very-short- 
movement devices limit motion to a fraction of a centimeter, e.g., a magnetic- 
latched reed switch. Lead wire bending should be kept to a minimum to prevent 
generation of spurious signals. It is convenient to simplify dynamic probe 
designs we have considered into two general groups. In one type, the temp- 
erature sensors are thermally short circuited to obtain a zero reading. In 
the other type, a single sensor is moved vertically up and down inside the 
probe to sense the temperature at various stations, and to indicate any off- 
sets of temperature sensors. 

Thermal Switch Probes 

The thermal switch probe concept is Illustrated in Figure 2. The two 
end sensors have been extended by "Isothermal" bars to points near the center 
of the probe. The open circle area Indicates a thermal switch zone. To avoid 
long recovery times, the two ends of the isothermal extensions should not be 
thermally shorted; therefore, a third "sampling" sensor must be incorporated. 
Such a sensor oust be light-weight (low thermal capacity), insulated from the 
probe walls , and mounted so as to sample first one end and then the other end 
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of the extension bars. This sensor could be differentially compared to a 
second sensor located in either the top or bottom isothermal section. If 
the top were chosen, then when the thermal switch is closed the sampling 
sensor is thermally shorted to the top section and no signal should be 
received from the sampling sensor. An automatic zero-offset calibration 
is thereby achieved. 

When the sampling sensor is shorted to the lower section, the probe 
temperature difference should be obtained. This difference may be compared 
to the value read by the fixed sensors and thus the offset of the latter may 
be checked. The thermal mass of the sampling sensor is sufficiently small 
so that no significant variation In the temperature of the upper and lower 
sensors occurs during switching. 

The fixed sensors are redundant if the sampling sensor operates satis- 
factorily. ’lowever, if the sampling sensor fails to operate, the fixed 
sensors would provide a back-up so the experiment need not be aborted. 

The thermal switch can be of several types; we considered two types: 

The sampling sensor can be centered so that an enclosed reed switch is 
magnetically latched to contacts on either the upper or lower extension bars 
(see Figure 2s) . 

The extension bar of the upper sensor can be made into an annulus with 
an insulated bar extending through its entire length (see Figure 2b) . Hear 
the end of this insulated bar is a high- thermal-conductivity section. The 
lower bar has a pocket for receiving the rod and the sampling sensor is 
provided with an annulus to allow passage of the rod. The insulated rod 
can be driven by a solenoid on the fop of the probe. When the sampling 
sensor is bridged by the high-conductivity section of the rod to the upper 
Isothermal arm, it is thermally shorted to measure the upper sensor temperature. 







(a) 

Reed Switch 


(b) 

Solenoid Switch 
(Shown Shorting 
Sampling Sensor 
to Top Section) 


LEGEN D: 

0 Isothermal Section 

0 Insulated Section 

S Sensor 

H Heater 

Sk Thermal Conductivity Sensor 


SCHEMATIC DIAGRAM OF PROBE WITH THERMAL 
SWITCH AND SAMPLING-SENSOR 


Appendix A-Figure 
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Downward movement of cha rod causa* eh* high-conductivity section to nova 
through cha sampling sensor and thermally short to - . lower extension bar. 

One disadvantage of the solenoid switch Is the necessity of holding the 
sampling sensor In position rigidly. This can be done by axial springs. 

Alternatively, the sensor could be eliminated and Incorporated into 
the high-conductivity section of Che solenoid-actuated rod. Shorting can 
be accomplished as before. The short travel of the solenoid will reduce 
bending effects on the leads. 

These concepts satisfy the requirement of providing a zero-offset cali- 
bration for the temperature gradient sensors. Probes Incorporating these 
concepts require considerable development work, would be heavier and have 
a lower axial thermal resistance. 

Moving Sensor 

This concept would use a single movable senspr that Is raised and lowered 
Inside the probe sheath. To utilize this method, techniques must be available 
to move a sensor over distances of SO cm without introducing effects caused 
by the continual bending of lead wires. Figure 3 shows the moving sensor 
concept. Two fixed sensors are placed in the top and bottom isothermal 
sections. The center of the probe is enclosed by a thin-walled tub* inside 
of which another thin-walled tube moves vertically. A small sensor is attached 
to the lower end of this tube. This sensor is an absolute temperature-measuring 
device with a low heat capacity so that it rapidly reaches the upper or lower 
sensor temperatures without significantly affecting the sensor temperature 
at those locations. To calibrate, the sensor Is lowered to the bottom station 
of the probe. It Is then raised to the top station and another reading taken. 
The difference in temperature should agree with that read out on the fixed 


sensors . 
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As an alternate feature, other heat stations c.*n be added at Intervals 
along the probe. Betveen calibrations, the moving sensor can be located 
at any of these stations to yield additional temperature data. In addi- 
tion, the tubes can extend to much higher levels and the moving probe can 
be used to monitor amplitudes and phases of the "AC" thermal wave propagation 
downward from the lunar surface. The thermal conductivity can be determined 
at the top and bottom locations with two heaters. If the moving sensor falls, 
the moving probe degenerates to a static probe. If the fixed sensors fall, 
the device can still operate to produce data. 

The moving sensor probe may be used without fixed sensors at the probe 
ends but the reliability of providing the motion over long periods of time 
is such as to make It desirable to have the probe degenerate to a stacic 
probe if failure of the moving parts occurs . 

TEMPERATURE SENSORS 

Platinum Resistance Bridge Sensors 

For the static probe only the platinum resistance bridge Is a suitable 
sensor capable of meeting the stringent stability requirements. Figure 4 
shows the sensor design selected for the measurements of temperature gradient. 

The sensor design has been developed specifically for use in the probe. 
The physical configuration of the sensors, in the form of a hollow cylinder, 
allows them to be mounted In the probe with minimum difficulty and to sense 
temperature with a time constant not limited by the sensor time response; 
it also allows Interconnecting leads for other probes or sensors to pass 
unhindered . 

Each cylinder contains two platinum resistance sensors of nominal 
500-ohm resistance. Each resistance element is bifilar wound of 0. 0015- 
inch platinum wire on an inner platinum cylinder. Each element consists 
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of a long coll, 0.0 10- inch in diameter, laid on a ceramic subs Crate in 
such a way Chat only the bottom portion of each turn of the coll is em- 
bedded In the ceramic. This particular mounting technique ensures that 
the vire is pure after annealing and will remain pure. The design of the 
coll provides resistance to dimensional changes of the substrate and Induced 
strain changes of resistance, and assures a rugged mounting of the sensing 
element. 

It is this mounting technique to which we may attribute the outstanding 
sensor performance. The stability of the sensors in a bridge configuration 
was 2 millidegrees or better after meeting the LEM environment shock, and 
random and sine vibration requirements as well as thermal cycling between 
195°F and 355°F (See Figure 5) . 

The element is protected by an outer gas-tight platinum case, and the 
sensor cavity is filled with helium. The sensor pinch-offs are designed to 
minimise both leakage resistance and gas loss. The epoxy sealing technique 
allows a gas-tight seal with low electrical and gas leakage. 

The smaller temperature sensor is preferable for the thermal-conductivity 
measurements. It must be installed on the surface of the probe. T’ present 
design consists of platinum wire wrapped on a ring-shaped platinum sleeve and 
embedded in a ceramic cement. Because this sensor has to measure the temperature 
rise during a short time interval, no long-term stability is required. 

Later testing of these sensors indicated that their stability was nearly 
as good as the more elaborate sensor described above. The temperature his- 
tories of these "ring" sensors on the lunar surface indicate that they have 
held their calibration within a few thousandths of a °K. 

Alternative Temperature Sensors 


We also evaluated several different movable probe-type temperature sensors. 



PLATINUM SENSORS 2A, 2B, 3A, & 3B PER FORM ANCE -ROSEMOUNT DATA 


Temperature Change - Millidegrees 


Preliminary Test 
Preliminary Test 
Reference 

Random Vibration 
Mechanical Shock 
Acceleration 
Thermal Shock 
Random Vibration 

20 g Sine 
30 g Sine 

40 g Sine 



Appendix A-Figure 5 
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Because this concept permits zero-offset calibration, long-term stability 
of the temperature sensors can be relaxed. The sensors considered include 
the following: 

- resistance thermometers (platinum, nickel, carbon) 

• thermocouples and thermopiles (iron-constantan, chromel-constantan, 
copper-cons tantan, gold-cobalt-copper 

- solid-state devices (thermistors, doped silicon, germanium, galli- 
umarsenide) 

- mechanical resonators (quartz crystals) 

' property variation sensors (gas thermometers, nuclear magnetic resonance) 

/hen we applied selection criteria such as documented performance data 
In the specified mechanical and lunar environment, compatibility with elect- 
ronics, development status, and availability from reputable suppliers to 
meet delivery schedules, none of the temperature sensors could approach the 
ovir-all suitability of the platinum resistance temperature sensors. The 
silicon carbide semi-conducting sensors and quartz crystal oscillators were 
the most promising of those considered although they did not meet several of 
the selection cri f . <ria. When installed in the probe cable, thermocouples 
can be used to measure the temperature changes within the borehole. 

The Preferred Probe Design : 

Because of its simplicity and the demonstrated stability of the platinum 
resistance; thermometer the static probe design seems the best approach at this 
time. The design of a static probe has been further advanced and is described 
in more detail below. 

Tvj separate static probes are connected by a flexible cable so that one 
is irectly above the other when installed in the borehole. This configuration 
can provide the desired data with greater reliability than a single static probe 
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The two probes, of Identical construction, are connected with a flexible 
cable which permits the whole assembly to be bent into a "0" shape. In this 
manner the probes and cables can be packaged to meet the requirements of the 
ALSEP scientific equipment bay. 

The probes consist of epoxy-fiberglass tubular structures which support 
and house temperature sensors and their associated electrical wiring. The 
probe structure is designed to hold the components in a fixed relationship 
and to protect the sensors, heaters, and internal leads from the Imposed 
mechanical and lunar environments. The structural 'materials have been chosen 
to provide a minimum axial thermal conductivity consistent with minimum weight 
requirements. The probe structure consists of two concentric thin-walled 
filament fiberglass-epoxy tubes, spaced by five annular braces which are 
bonded in place. 

Thermal Conductivity Measurement 
Design Criteria 

A method for determining lunar material thermal conductivity in the 
borehole is a required part of the heat flow experiment. The criteria which 
we have used in evaluating various design possibilities are: 

(1) K Measurement 

moon .3 .5 * 

Range: 4 x 10 to 2 x 10 watt/cm K 

Accuracy: greater than 20Z 

(2) Power Available 

100 mw. (for 1 month) 

I watt (for less than 14 days) 

(3) Time Available 

1 month or less if possible 

* This range recognizes the possibility that an experiment might be surrounded 
by a solid rock in the regolith. 
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(4) Temperature Measurement Cabability 

Range: Q-2°K differential (possibly two fourfold-amplification 

stages) 

o 

Accuracy: +0.001 K 

o 

Maximum Environment Temperature Allowed: ^ 400 R 

(5) Temperature Read-Out Intervals 

* 

10 min (for single point read-out: 1 min) 

(6) Design 

The device must be a compatible part of the heat flow probe, 
and should meet all operational requirements for the static probe. 
Because the probe will not be in physical contact with the lunar 
borehole, heat transfer between the probe and the lunar material 
will be by radiation only. If necessary for obtaining a conductiv- 
ity measurement, solid contactors could be added to portions of the 
probe . 

(7) Multiple Determinations 

If possible, to indicate either homogeneity or stratification. 
Experiment Concepts : 

fc 

Thermal Conductivity (k) and/or thermal diffusivity (~ c ) may be determined 
by placing a controlled heat source within a medium and measuring temperature 
at a point in the system. A heater-temperature sensor arrangement which appears 
most suitable for the present application is shown in Figure 6. Since the 
possible range of lunar material thermal conductivity values extends over two 
orders of magnitude, tiro modes of operation are required. 

In the "Low K" mode the heater is turned on and its temperature is measured 
as a function of time. A typical power level of 1-5 milliwatts produces about 
a 2°K dynamic range requirement for the temperature sensor. The steady-state 
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temperature level Is related to both the lunar material thermal conductivity and 
the resistance across the radiation gap between the heater and the moon. If 
the radiation resistance can be estimated with sufficient accuracy, then the 
thermal conductivity may be computed. 

Should the lunar material be highly conductive, then the radiation gap 
resistance becomes large — relative to the effective lunar resistance. A 
small error in predicting the radiative flux can be amplified by as much as 
twenty times In the thermal conductivity estimate. Although the gap resistance 
could be reduced by the use of solid contactors, the uncertainty in estimating 
contact conductance becomes very great due to the strong dependence of the 
conductance on the unknown structure of the borehole surface. 

Therefore, if the Low K mode temperature level indicates a high lunar 
conductivity (which, however, cannot be resolved -within a + 10 per cent accuracy) , 
a second type of thermal conductivity experiment Is carried out. In Hi-K 
mode, the heater power Is suddenly raised to a continuous level of about one 
watt. Temperature sensing is initiated at a location on the probe about 8-10 
cm away from the center of the heater (again for a 2®K dynamic range)- For 
the Ideal case of a perfectly Insulating probe, the temperature monitored in 
the moon at a specified distance away from the heat source will be determined by 
only the rate of propagation of the thermal wave through the lunar material. 

The time derivative of this temperature-time curve will have a maximum value 
which is inversely proportional to the density-heat capacity product (pc) of 
the moon and a peak width or position on the time axis which is related to 
thermal dlffusivlty. If these two quantities can be measured with sufficient 
accuracy, thermal conductivity can also be determined. 



In a praccical system, several real effects decrease the sensitivity of 

Che measurement. The probe section between the heater and the Hi-K mode sensor 

must be designed to prevent the effects of thermal shorting down the probe fzom 

masking Che signal in the moon. In addition, the thermal mass of the Hi-K mode 

sensor must be small, so chat thermal inertia in the sensing device does not 

alter the transient response of the data signal. Both of these real effects have 

been considered and it appears feasible to use the Hi-K mode in the conductivity 
-3 -4 

range of 4 x 10 to about 4 x 10 watts/cm°K. For lower thermal conductivities, 

the shorting effect of even a well designed probe becomes significant. The 

time required for the measurement in this case becomes appreciably longer than 

several hours and the large amount of energy transferred to the moon requires 

an excessively long time to dissipate (several weeks). For these reasons, use 

of the Hi-K mode is suitable only for measurements at the high end of the 

possible conductivity range. However, this is the range where Low-K mode 

measuring capability is poorest. Used together, Che two methods complement each 

ocher and can meet the experimental requirements. 

The time required for making the thermal conductivity measurement is less 

than a month for both methods. For low conductivity, measuring times can 

become long, since an approach to a steady-state temperature level is required. 

The time scale is longest for the lowest thermal conductivity to be determined 

— about 10-20 days. In the Hi-K mode, the time scale of the experiment is 

greatly reduced, both because of the higher thermal conductivity and because 

transient, rather than steady-state, data are obtained. However, for a lunar 

-4 

material thermal conductivity of about 4 x 10 watt/cm a K a heating period of 
about 6 hours at a 1 watt source level is estimated. The time required for 
the gradient resulting from the heating to decay to five percent of the undis- 
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Curbed lunar gradient ranges from about 20 days for the case just mentioned down 
to about two days for a material of high thermal conductivity. 

Design of Thermal Conductivity Apparatus 

The components required for measurement of lunar material thermal conductivity 
by the methods described are: 

Lo-K mode : A heater which is integral with a stable, sensitive temperature 

sensor and is radiatlvely coupled to the moon. 

Bi-K Mode: A heater located about 8-10 cm away from a sensitive temperature 

sensor. The sensor should be lightweight (M. gm) , but long-term stability is not 
required. The heat flow down the probe section should produce a temperature rise 
at the corresponding point in the moon — the probe wall should not act as a 
thermal short. 

The design- is indicated*' in Figure 6. One such conductivity instrument is 
located at each probe end. The heater is mounted directly on one of the main 
platinum temperature sensors. The sensor has been assumed to have a dynamic 
range of 2 9 K and a sensitivity of 1:2000. The heater winding forms a cylin- 
drical heat source of about 2 cm diameter and 2.5 cm in length. Initial es- 
timates, which will require further refinement in terms of either a detailed digital 
or analog simulation of the actual system, indicate that the power dissipation in 
the heater will be about 1-5 milliwatts for the Lo-K mode operation (low conduc- 
tivity range) and about 1 watt for the Bl-K mode (high conductivity range) , On 
Che basis of a 29 volt (maximum) power supply, the heater winding should have 
an electrical resistance of about 840 ft. 

The Bi-K mode sensor, located about 10 cm away from the center of the 
heater, can have the same range and sensitivity as the main bridge sensor. 
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The smaller sensor will also be a platinum bridge designed to have a total weight 
of about two grams, with stability required only for periods of several hours 
(as opposed to the larger main sensors which must be stable for a year) . The 
sensor will be mounted In a conducting band (width about 0.5 cm) which encircles 
the body of the probe. 

The two small sensors on each probe (one coward each end) will each form half 
of a platinum bridge. Only one heater will be run at a time and the farther 
sensor will act as a reference. The wires connecting the bridge will be run 
along the inside of the probe body wall to prevent equilibration with the main 
cable bundle. Power and signal leads to the bridge will be introduced at the 
midlength point where thermal shorting from the heater down the main cable 
bundle will be at a minimum. 

The inner tube which holds Che cables may also act as a thermal short. To 
reduce '■his effect, the surfaces of the cable sheath and the inside of the 
probe body are silvered in the region between the heater and the remote sensor to 
prevent radiative transfer to the probe body from the cables. (These surfaces 
are black in regions outside of the Hi-K mode experiment zone) . The outside of the 
probe body is black to permit heat loss from the probe wall to the hole wall, 
thereby greatly reducing the thermal shorting down the probe. 

Pinal Design. Calibration and Data Interpretation 

Although initial probe design has been based on valid, but approximate, 
models which can be treated analytically, a detailed model including all the 
interrelated effects of heat transfer between and within probe components as well 
as between the probe and various types of possible lunar materials is reouired 
for orobe calibration and accurate Interpretation of data. 


1 . 
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COMPUTED HEATER TEMPERATURES FOR MODE A CONDUCTIVITY EXPERIMENT 
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Such « model can b« el char t digital computer program for the actual system 
or an electrical analog simulation. The latter method has much to offer In 
terms of ease of formulation, low cost and Inherent flexibility, but results 
are limited In accuracy by the use of components of + 1 percent precision. 

Since the temperature sensors are designed for a 1:2000 resolution, the electrical 
analog la not sufficiently accurate for measurements requiring full resolution 
over the entire dynamic range of the sensor. In the Hi-K mode conductivity 
experiment, full temperature resolution Is required so chat accurate values 
of the temperature derivative can be obtained. However, errors in evaluating 
Che temperature derivative can be avoided in the electrical analog by direct 
measurement of current (proportional to the derivative of voltage) . The 
electrical analog therefore appears to be suitable at least for final de&ign 
purposes and for calibration. A digital model may also be used, and although 
model formulation probably would require a major programming effort, the 
computer simulation would be more accurate and more suitable for final inter- 
pretation of results. 

A digital computer program, written and run by Professor Sydney Clark 
(coinvestigator. Tale University) has yielded temperature-time curves for 
various system designs over the anticipated range of lunar parameters. 

Although the results have been utilized to demonstrate the effect of various 
design parameters for different lunar materials the model is approximate jince 
the probe has been created as a solid cylindrical rod. 

Figure 7 presents curves computed by the Clark program for an effective probe 
conductivity of 10 cal/cm sec*K. This is equivalent to a probe thermal resis- 
tance of about 800*K/watt cm which is quite close to the anticipated thermal 
resistance for the actual probe. (The probe body is the main conduction path 
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since the cables and center rod are weakly coupled to the wall in this region.) 

In the Lo-K mode operation Clow thermal conductivity range) , steady-state heater 

temperature level is used to find thermal conductivity. Although the program was 

only run for 14 hours (real time) , it is evident that steady state has been reached 

for the ' rock f material and is being approached quite closely for the f vesicular 1 

material. At steady-state , lunar density should not affect the temperature level. 

The difference due to density for the 9 powder’ material appears large, but is a 

manifestation of strong transient effects which would disappear after about 

twenty days (estimate of steady-state time based on computation for similar spherical 

source). The stea< y -state temperature level for this case would be about 3*K. 

A large separation in steady-state level due to variation in thermal conductivity 

-5 -4 

occurs in the range of 10 <k < 10 cal/cm sec°K. However, resolution by the 
Lo-K mode becomes poor in the high conductivity range. The initial temperature, 
indicated at ten minutes, is due primarily to the temperature drop across the 
radiation gap between the heater and the moon. Small errors in predicting the 
radiation resistance can vary the steady-state level in this f rock f material 
range so much that a + 10 percent estimate of lunar conductivity cannot be made. 
Although increasing heater power would increase the temperature amplitude of the 
results, the relative error would remain the same. 

Consequently, Hi-K mode operations should be used if data from the Lo-K 
mode indicate a high thermal conductivity. The Clark model was used to esti- 
mate temperature at the probe surface for various distances between heater and 
sensor locations. Figure 8 shows the time derivative of computed temperatures 
at a point 8.075 cm away from the certer of the heater* The original temperature 
curves at t lis point are similar in shap'_ to those in Figure 7. For a one watt 
heater output, the peak amplitudes of the derivative curves range from about 0,010 
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to 0.035‘K/oin ^related to density values for lunar material from 0.5 to 2.0 cm/cc. 
Since a one minute time Interval Is small compared to the time scale of the peak, 
temperature differences should be taken on the basis of time Increments of several 
minutes to reduce error due to the + 0.001*K uncertainty In temperature deter- 
mination (l.e., the temperature- time curves should be smoothed before their 
derivatives are taken). 

The time at which the maximum occurs (or some other characteristic time, such 
as that elapsed between the half-amplitude and full-amplitude points (can be used 
to compute thermal diffusivlty. Therefore, thermal conductivity Is related to a 
ratio bet wetu the time scale of the peak and the peak amplitude. 

-3 2 

In Figure 8, the thermal diffusivlty for the lower peak is 3 x 10 cm /sec 
-3 2 

and for the higher, 10 cm /sec. The separation between p&aks is In about the 
same ratio. The. Clark model Indicates that thermal shorting through the probe 
Is not serious. Less resolution would be obtained for a real probe because of 
two effects: (1) the concentrated thermal mass of the sensor, and (2) the 
temperature averaging due to the finite width of the sensor ring. With a sensor 
weight of 1-2 grams and a width of 0.5 - 0.8 cm, adequately resolvable data 
should be obtained. However, data interpretation using a fairly exact computer 
Cor analog) model would be necessary to compensate for these effects. 



NOTE: Prof. Clark's Model, Runs No. 56005 and 57005. 

Sensing Location Is 8.075 cm from Heater. 

Probe (Assumed as Solid Cylindrical Rod): k = 10"* cal/ cm sec°K 
Equivalent Probe Resistance = 800 ®K/ Watt - cm Length 



COMPUTED TEMPERATURE DERIVATIVES FOR MODE B CONDUCTIVITY EXPERIMENT 
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APPENDIX B. THE DEVELOPMENT OF THE LUNAR SURFACE DRILL 

What follows Is a brief history of the development of the Apollo Lunar 
Surface Drill (ALSO) . This tool went through several stages of development 
and was the main preoccupation of the Lunar Heat-Flow Experiment. The basic 
design for the heat-flow probes was decided by the end of 1966 and, compared 
to the drill development and testing, vent smoothly. With the static probe 
design which was preferred it was clear that a hole in which it would be 
placed had to be predrilled. Based on the range of conductivities that were 
thought likely to characterize the lunar regolith, it was determined chat 
the meter long probe sections should be burled between two and three meters 
below the surface. At these depths the surface variations of temperature 
would be barely detectable. 

Two problems loomed large and were difficult to define: 

1. Would solid rock be encountered; and if so, how much should the ALSO 
be designed to penetrate to guarantee a three meter hole? 

2. Would the hole collapse after the drill was extracted? 

Results from Surveyor Missions Indicated that rocks larger than a few 
centimeters were prevalent. What the distribution of solid rock fragments 
in the subsurface was, was a matter of conjecture. If they were as numerous 
as they are on the surface, then there was a finite chance that a rock 20 cm 
or less would be encountered when drilling a three meter hole. Thus, we 
specified that the ALSD should be capable of penetrating at least a 20 cm 
section of solid rock with a hardness similar to terrestrial basalts. 

The soil mechanics tests carried out on the Surveyor Mission indicated that 
the lunar soil was not very cohesive; consequently, there was some likelyhoed 
that the hole would collapse after the drill was extracted. The probe design 
selected, with two tandem probe sections flexibly connected would be extremely 
hard to insert into an open bare-walled hole i.. any case. Consequently, we 
concluded chat some means to sleeve the hole would be required. 

H ard Rock Drilling 

From the start of the drill development, the working design put forward 
by Martin Marietta ( the prime contractor for the ALSD) was an open faced bit 
of rather narrow kerf. The outside diameter of the hole would be about one inch 
and the ID of the drill was 0.8 inches. Such a design would, of course, obtain 
a core sample of the material through which it drilled. 
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Volume-Energy Relation for 60° and 90° Wedge Bit* In 
Drop Teat* on NASA Standard Basalt Sample. 
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For the bit to bo able to cut through solid crystalline rock, using 
relatively low power and a small downward force Can astronaut on the lunar sur- 
face could apply a downward force of only approximately 30 lbs) , a rotary- 
percussive action was used. If properly designed, this allowed the bit to 
actually chisel the rock face and not grind it away as rotary diamond drills 
do. The optimum design for a bit face was achieved after considerable testing. 

An excerpt from the Martin Marietta final report follows, describing the 
developments and testing of a core bit and conclusions. 

" Cora Bit Optimisation 

The objective of these tests was to develop a core bit with a cutting tip 
geometry and material which would drill reliably (minimum degradation) and 
efficiently (maximum penetration rate with minimum power consumption) when 
powered by a light-weight, low energy power head. These characteristics were 
required for all lunar material simulants such as dense basalt, vesicular 
basalt, scoria, pumice, and conglomerates. 

Optimization of the core bic was primarily accomplished using the dense 
basalt as a standard because of its tougher drillability characteristics and 
its homogeneity, which was • conducive to the performance of repeatable tests. 
Subsequent tests indicated that the vesicular basalts, although easier to drill, 
can often result in greater degradation of the core bit chan the dense basalt. 
This factor had to be taken under consideration during design finalization. 

Dense Basalt 'Standard* Characteristics - The NASA furnished dense basalt 
standard, which represented the maximum drilling energy requirement of the 
lunar material simulants, was subjected to a series of tests (15) to determine 
its characteristics. The uniaxial, ultimate, compressive scrangch of the sample 
was measured at 21,732 psi (average) and the average shore hardness was 41.6 
with a standard deviation of +23. Drop tests were conducted in accordance 
with standard procedures using 60° and 90 s wedge bits in the energy range of 2 
to 20 foot pounds as shown in Figure 1 . . Using the data for the 90° bic, it 

was calculated chat Che ALSD operating at a percussive blow rate of 2250 bpm . 
would require a bit delivered energy of 27.6 inch-pounds per blow tu maintain 
a 2 inch /minute penetration rate. Assuming a power head output energy of 30 
inch-pounds, a cuttings scavenging loss of 202 and a 102 energy loss per extension 
tube, the following empirical ALSD penetration rates were calculated: 
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Number of Extension with Bit Penetration Rate (la. /Min.) 

1 1.26 

2 1.09 

3 0.93 

4 0.81 

These predictions did not take into account the benefit derived from 
optimization of blow indexing, which, at the time of initial power head 
design, was believed by most researchers to be of little or no benefit. 

Core Bit Test Summary - Throughout the performance of this program, 
approximately 100 development bits were tested. In general, the bits were 
designed and fabricated in groups of 6 to 15 each, tested, and the results 
of the tests and redesign recommendations forwarded to the core bit contractor 
for Incorporation in subsequent test groups. Major design variations included 
tip geometry and tip material (tungsten carbide) chemistry. Minor variations 
in tip brazing techniques and core bit body geometry were also studied prior 
to selection of the final design. 

Most of the optimization tests (16-21) were conducted using Che equipment 
illustrated in Figure 2 . With this arrangement it was possible to accura- 
tely control the axial bit force and power head input voltage while monitoring 
rpm, penetration rate and power consumption of the system. 

Typical bit configurations tested included those illustrated in Figure 
3, where 3, 4, and 5 tips per bit of a single geometry and tungsten carbide 
chemistry were tested to investigate the effects of index variations. The 
advantage of indexing was clearly demonstrated during these particular tests. 

Tests clearly demonstrated that the 3 and 5 carbide tip bits outperformed 
the 4-tip version due to improved indexing. At this point in the ALSO develop- 
ment program, a decision was made to replace the original 4-tlp design bits 
with the 5-tip configuration. Although the 3-tip design indicated slightly 
superior performance, the bit tended to 'chatter' excessively during operation 
and it was anticipated that rapid degradation of the dps would occur. 

Another factor demonstrated during the performance of the core bit test 
program was the importance of optimum axi.-J. bit force to obtain maximum 
penetration rate at a minimum power consumption. Data clearly illustrate 
the significance of optimum axial bit force in dense basalt. Bit force is 
not as critical in softer rocks such as scoria or vesicular basalt, because 






of che lover energy threshold required for predictable fracturing. 

F inal Core Bit Selection - The core bit test program indicated that a 
compromise design vas mandatory in order to meet the ALSD drilling require- 
ments in all lunar simulant materials. The harder tungsten carbides (91.3 RA) 
possess good wear characteristics In dense basalt, but tend to fracture if 
operated under light loads. Tip geometries which perform at high pene- 
tration races and low power consumption tend to degrade rapidly. The final 
basic configuration selected is detailed in MMC drawing PS600100023-11. A 
minor modification was incorporated, as shown on MMC Drawing 467A8050000-0H, 
to ensure that a uniformly round material (89.3 RA) consists of 872 WC and 
132 cobalt. ". 

Drill Stem 

To drill a three meter hole additional tubing would be needed to produce 
a drill stem somewhat longer than three meters. The wall thickness of the 
tubing had to be less than the fit's kerf width of 0.1 inches (2.54 mm). To 
remove cuttings from the hole while drilling, spiral fluting was required on 
the exterior wall of the tubing. To fit in the stowage bay on the LEM it 
would be necessary to break che drlllstem down into lengths no greater than 
22 inches (55.88 cm). These requirements called for a high strength. A 
titanium alloy was then found to be the best suited to all of these require- 
ments and it had the additional advantage of light weight. 

A critical part of the design was the joint to connect and disconnect the 
drill stem sections. The Martin Marrietta design was a threaded type joint that 
could be made in a single turn. When mated, the flute continuity was preserved 
and che stem walls remained flush. A special lubricating coating was applied 
to che joint surfaces so chat they would not seize up during drilling. 

Raving chis type of joint produced some serious human engineering problems 
for the astronaut. To uncouple che drill string as it was retracted, a 
special treadle was designed. This treadle vas restrained by the astronaut's 
foot as he drilled. It incorporated a hole to guide the stem and a built-in 
clamping arrangement that restrained the drill stem at the treadle while the 
astronaut disconnected the corehead and che stem sections filled with core 
material. 
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Powerhead 

To provide che rotary and percussive motions to the drlllstem, a battery 
operated powerhead was developed by Black and Decker Manufacturing Company 
for the ALSO. The powerhead consisted of an electric motor of about a half 
horsepower which provided torque to the drlllstem coupling and recocked a 
spring driven hammer, which fell on an anvil coupled to the top of the drill- 
stem during each revolution. The unit was housed In a pressure-tight case 
so that It could maintain nitrogen gas at atmospheric pressure to help main- 
tain the Internal parts at low temperature. 

S aving the Hole 

Originally, the drill system was designed without any provision to Insure 
that the hole stayed Intact so that the heat-flow probes could be inserted. 
Consultations with scientists working with the surveyor results Indicated that 
cohesion of the regolith was too low to support a bare-walled hole and we were 
advised to provide some type of sleeving. 

The first attempts to provide such sleeving used chin walled fiberglass 
tubes chat joined together by force fitting matching tapers cut into the ends 
of the tubes. The bottom tube was fitted with a conical head. The sleeving was 
designed to be powerdrllled Into place after che hole had been completed. 

Some typical designs that were tried are shown In Figures 4 and 5. 

R egolith Model for Testing the ALS O 

Early in che development program It was clear that a full scale test bed 
of che lunar regolith would be needed to evaluate che drill performance. 

The first such test bed was constructed at che Martin Marietta plant in 
Baltimore, Md. Lamont-Doherty Geological Observatory undertook che design 
of the three meter deep regolith monel. The basic Ingredient of the model 
was powdered basalt. Using particle-sire distributions determined from 
Surveyor phocographs, the appropriate size fractions were calculated in terms 
of standard sieve and screen sizes. The components were Chen produced by a 
local (Haverstraw, New York) crap rock company and stored in 55 gallon(208. 18 
liter) drums. Large blocks of basalt, to be placed at known spots in the 
model, were ordered from Oregon since che flood basalts there had been used 
as a standard in earlier tests. 
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A large rectangular steel tank has made at Martin Marietta to contain the 
model. A diagram of the model is shown in Figure 6. The unsorted part 

of the model was made by remixing the components according to the Surveyor 
recipe in a standard cement mixture. Layers of roughly a foot were poured 
Into the tank, large chunks CIO cm or greater) were thrown in by hand. Each 
layer was tamped down which made the fabrication of the model dusty and un- 
pleasant work. The completed model was used extensively to test the ALSO. 

Tests in this regolith model showed that the drill was capable of pene- 
trating three meters Including 50 cm of vesicular basalt within the prescribed 
time and power available from the silver cell battery pack. However, the 
design for installing the casing was far less successful. The casing would 
frequently break, or large torques (enough to stall the motor) would build up 
before the casing sleeves had penetrated the full depth of the drilled hole. 

Field Tests 

We were also desirous of carrying out tests under field conditions where 
we wbuld be drilling into an unknown subsurface. Finding an area that would 
be a reasonable analog to the lunar surface was challenging. It would Le 
important that the soil contain virtually no moisture to depths on the order 
of three meters. We looked at seme areas around Flagstaff, Arizona which had 
been selected by the Lunar Geology group :t Flagstaff as lunar analogs. There 
were ash flows and lava flows but the water content of the soil was just too high. 

We decided instead to try Death Valley, hoping that the extremely dry 
conditions would result in bone dry soil layers. 

We found three very interesting targets. 

1. Small pedestals of calcareous lacustrine deposits were abundant In 
the Shcshone Valley. Since the diameter of these pedestals was no more than 

10 to 15 meters and the heights were about five meters, there was a good chance 
that these very fine grained sediments were quite dry. Of course, the composition 
of these sediments was not lunar-like. 

2. Large fans of loosely compacted conglomerates, which had formed at the 
mouth of small ravines cut into the hills bordering the Shoshone Valley and 
Death Valley . 

3. Alluvial and fluvial deposits of silts and sand in arroyos. Vertical 
banks had been cut into these deposits by the most recent streams, so that a 
vertical section of what was drlllea could be seer.. 
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3- Meier Conglomerate and Pumice Model Construction 
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Drilling was attempted in the lacustrine deposits. At about 2.S meters 
the drill experienced stalling torques. This was because of significant 
moisture at those depths. We found it extremely difficult to remove the drill, 
it could not be moved by hand, so we jury-rigged a large lever to force the 
stem out of the ground. Once the drillstem was removed, the hole stayed 
intact so that no casing was needed. 

In the fanglomerates drilling was more difficult because the drill would 
occasionally encounter a cobble 10 cm across, but the holes could be drilled 
and the stem removed without problems. The holes made were not preserved. 

We then tried reentering with the fiberglass sleeving and encountered the 
same problems we had earlier in the model at Baltimore. 

The main result of these field tests was that the drill was capable of pene- 
trating the full three meters in a variety of soil types that in most cases 
presented more difficult drilling conditions than on the moon. 

The methods for casing the hole were not yet reliable and design Impro- 
vements were needed. 

The Boron Reinforced Borestem and the Solid Faced Bit 

The principal investigator's group saw a possible means of overcoming the 
problem of the hole collapsing. 

One means of accomplishing this would be by using a drillstem that does 
not have to be retracted from the lunar subsurface. Such a drillstem must 
have a low thermal conductance, so as not to degrade the heat-flow experiment 
performance, yet have an axial modulus high enough to effectively transmit the 
percussive energy of the powerhead to the bit. To utilize a low conductance 
drillstem that can be left in the lunar subsurface it is, of course, necessary 
to develop either a solid-face bit that does not allow material to enter the 
drillstem, or, if a coring bit is used, to develop a core retrieving system 
that can remove material from inside the stem. 

The approach that seemed most feasible was to replace the titanium drill- 
stem with a stem made of composite material. An epoxy resin was chosen for 
the matrix of the stem tube because of its low thermal conductivity. The 
epoxy is strengthened by glass fibers wound circumferentially in the tube walls. 
The axial stiffness is increased by an order of magnitude by including fila- 
ments of crystalline boron (YoungLs modulus 60 x 10^ psi) , aligned with the 
tube axis. 
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We hoped eo develop a solid-face bit by a modification of the existing 
AL5D coring bit designed by Martin Marietta and Chicago-Latrobe , thus bene- 
fiting from the experience and testing that vent into that bit development. 

Our approach, which proved successful, was to insert a central cutter inside 
the annulus of cutters of the core bit. When drilling solid rock, the core 
bit leaves a small cylinder of material standing above the cutting face, this 
cylinder is broken up by the carbide central cutter. 

An alternative to the solid-face bit is to remove the core material of 
the subsurface core samples, a system that could retrieve the sample would 
be worth the extra astronaut effort required. To this end, we began develop- 
ment of a 'core retrieval system. ' We decided to utilize a system similar to 
that used in standard drilling techniques. The basic feature of the system 
is a thin-walled core liner that is locked in the bottom section of the drill- 
stem. These core liners would be Inserted and locked in place. After drilling 
into the lunar subsurface the length of the core liner, it is unlocked 
and removed by the astronaut. For emplacement ’nd retrieval of the liners 
a special tool would have to be designed. This liner serves as a core sample 
container on the return flight. 

The AVCO Space Systems Division, under a subcontract from Lamont-Doherty 
Geological Observatory, was initially chosen to perform the design and develop- 
ment of the composite drillstem. The work was done under the guidance of 
Lamont-Doherty. The subcontract with AVCO consisted of two tasks: 

Task I: To conduct a design and manufacturing feasibility study. 

Task II: Fabrication of prototype sets of drillstem. 

Phase II was completed under a separate subcontract with Arthur D. Little, 
Inc. This work was also done under the supervision and direction of Lamont- 
Doherty Geological Observatory. 

The development solid-face bits were purchased by Lamont-Doherty from 
Chicago-Latrobe, Inc., based on Lamont-Doherty Geological Observatory drawings. 

Program Developments and Test Results 

The boron filament reinforced fiberglass drillstem: 

The boron drillstem developed under this program is shown in Figure 7. 

The tubular body of the drillstem has a sandwich construction (the details of 
which are. shown in Figure 8. Innermost are two layers of epoxy fiberglass 
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with the glass filaments helically wound at + 45* to the stem axis. Three 
layers of boron filaments are wrapped over these layers that are aligned with 

i 

the long axis of the tube. Two additional layers of fiberglass are wound over 
the boron at angles + 45® to the stem axis. This configuration gives a 
nominal outer diameter of 0.990" with an inner diameter of 0.875". 

The flutes to transport drill cuttings to the surface are helically wound 
around the stem body as shown in Figure 7. The flutes are also made of 

j fiberglass and form a double helix with a one Inch pitch. The outer surface 

of the flutes are Impregnated with silica powder to increase abrasion resis- 
tance. Each flute Is 0.10" to 0.15” across and stands about 0.45" above the 
stem body. As can be seen in Figure 7 , the diameter of the stem body is 

increased at the female joint, whereas the flute diameter is constant, and as 
a result the flute depth over the joint area is very small; i.e. , about 0.01"., 

Conductivity measurements on the boron stem were made at Arthur D. Little, 
Inc. by request of the Lamont-Doherty Geological Observatory. These results 
give a value of conductivity equal to 0.0096 watts/cm°C at 27°K and 0.018 
watts/cm’C at -73°K. 

Designing a workable taper joint for the boron filament drlllstem proved 
to be the most difficult task In the development of the boron drlllstem. The 
initial taper joint design. was a simple taper of 0.012 inches/inch for a 
length of 1" machined Into the outer glass layers of the male and female pieces. 
The Inner part of the tube containing the boron filaments was ground to form 
a flat shoulder that butted when the joint was made up. This type of joint 
failed during testing due to lack of adequate hoop strength; also the boron 
filaments sheared off at the shoulders. 

These tests led to a basic revision of taper joint design. The steel 
reinforcing rings were removed, and the female taper was reinforced by 
overlaying it with short lengths of boron filaments. In addition, circum- 
ferentially wound fiberglass was applied around the female taper joint to 
increase hoop strength. The taper angle was increased from 0.012 inches/inch 
to 0.018 inches/inch to decrease radial stresses in both the male and female 
sections of the taper. This joint configuration has proven successful, and 
is the basis for the final design. 

Many tests of the taper joints were made In the course of development. 

For the final taper joint design we ran an endurance test that lasted for 22 
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Appendix B-ri^ure 9 

bevlldlh °i S ° f lid faCed dri11 bit - Note ths cutter preluding 
beyond the kerf cutter at left. ^ * 

Appendix 3- Figure 10 Solid faced drin bit, final design. Three of the six elliptical 

exit ports for core cuttings are visible. 
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olnueas of drilling In danse baaalc with no apparent joint degradation. 

The Solid-Faced Bit 

1. • Experience gained In developing the ALSO coring bit provided Infor- 
mation on bit design parameters. Point pressure exerted on the rock by the 
carbide kerf cutters oust be kept above a finite threshold value to obtain 
chipping at the rock face. If point pressures are below the threshold, the 
rock Is pulverized rather than chipped and the drilling rates will be much 
lower than those obtained with the ALSO bit. 

2. The design basis of the solid-face bit was the ALSD coring bit. In 
the first solid bit design, a single tungsten carbide blade was placed inside 
the bit shell to remove the core. On the first prototype this blade protruded 
beyond the kerf cutters (see Figure 9) . 

Drilling rates obtained with this bit were extremely low because of low 
point pressures. A second test was made with this bit to determine the effi- 
ciency of the core cutter. An ALSD coring bit was used to drill a hole In 
dense basalt, leaving the central core undisturbed. This hole was redrilled 
with the solid-face bit as a test of core cutter efficiency. The central 
core was removed at a high rate. This test showed that it is relatively easy 
to break up the central core once it is standing about half an inch above the 
cutting face. The reason for this is that the free cylindrical surfaces lead 
to fracturing of the column both in shear and spallation. This discovery lead 
to a final, highly efficient solid-face bit design. 

The second prototype solid-face bit also used the ALSD coring bit. How- 
ever, the blade was recessed about 3/8 of an inch behind the kerf cutters of 
the core bit. Six elliptical holes were machined into the body of the bit at 
the base of the core cutter to serve as exit points for the core cuttings. 

(see Figure 10) . 

Tests of this bit were successful, giving rates of 4 inches/minute in vesi- 
cular basalt and 1 inch/minute in dense basalt. This bit configuration has 
been chosen for the flight hardware with only minor modification. 

A Core Retrieval System 

The core retrieval system was designed to collect a subsurface core sample 
without the necessity of removing the drlllstem from the borehole. This 
system was never used as part of the ALSD. This is accomplished by locking a 
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cor« liner cube Into che lower end of che drill stem, drilling co a depth 
approximately equal to the length of the core liner cube and then removing 
the tube. 

The core liner was designed to function as both a core liner and a sample 
return container. Strict design requirements are placed on it as a result of 
this dual function. The core tube must interface with the drlllstem, core lock 
and sample return container and meet the aseptic requirements for returning 
lunar samples. 

The prototype core liner is fabricated from thin-well, high strength 
aluminum alloy tubing. Testing of early designs showed that the tubes were 
deficient in two areas: 

1. The presence of rock particles which became lodged between the core 
tube outer surface and the inside of the drlllstem made removal of the core 
assembly difficult; and, 

2. core material was frequently lost during removal of the liner from 
the stem. 

Core liner jamming was eliminated by Installing a felt-flocked paper seal 
at .the bottom of the core liner (see Figure 11). This seal effectively 

limits the entry of rock and rock dust particles Into the core liner-drlllstem 
interface .id cleans that interface during core liner Insertion. Loss of 
cored material was minimized by using a core liner fitted with a core catcher. 

Core Lock 

A core lock was required to hold the coring assembly in place during the 
drilling operation. 

Our initial core lock design utilized phosphor bronze spring fingers , which 
were forced against the drlllstem wall at an acute angle. A screw driven cam 
actuated the spring fingers. Upward force on the coring assembly caused the 
spring fingers to slightly penetrate the drillstring wall thereby preventing 
movement of the coring assembly. Tests showed that these fingers were not 
nearly strong enough to hold the core liner in place. 

A second development core lock used a ball detent mechanism. Two 0.250" 
diameter steel balls were contained in the core lock body. Detent holes were 
made in the drlllstem wall to accomodate the steel balls. A screw driven cam 
forced the balls outward into the detent holes. Reversal of the cam motion 
allowed the balls to retract during core retrieval. Tests in dry rock powders 
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Appendix 3-rigure 12 p arts 0 f a prototype core lock of a design similar to the final 

design. 


, „ The final design core lock with half of knuckle joint attached. 
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Appendix B-Figure liThe knuckle joint^core liner tube on left and core lock on right 
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revealed chac the core liner exerted such a strong upward force during 
drilling that the steel balla actually tore the fiberglass stem upward from 
the decent holes. 

A third and final core lock design utilized an expandable cylinder with 
a knurled surface to engage the interior surface of the stem. In the final 
design, the core lock utilizes only one taper plug and is articulated at only 
one end. The final core lock design is shown in Figures 12 and 13. 

The core lock and core liner are joined with a knuckle type joint mated 
by rotating one of the pieces through a 90° arc, picture in Figure 14. 

Emplacement-Retrieval Tool 

Emplacement, retrieval and core lock actuation are performed with a 
special tool. This tool must be designed to transmit tension, compression 
and bidirectional torque to the core lock capstan. The prototype unit 
consists of a torque limiting handle, five extension sections, and one adapter 
section. The handle is connected to the adapter section directly or with one 
to five extension sections interposed. The interstem joints are not separable 
once they have been connected except with a special tool. Ordinarily, there 
is no need to separate the sections on the lunar surface. 

The handle limits the torque applied to the core lock to a preset amount. 

Summary of Development Tests of the Core Ret ri eval System 

Development tests were carried out at the Lamont-Doherty Geological Observatory 
test facility using a Black and Decker model 723 rotary hammer run at a reduced 
voltage to simulate the ALSD powerhead. A fifty-five gallon drum, filled with 
material from the model Martin Marietta Corporation was used to simulate the 
lunar surface. Hard rock drilling was done in either a 17 inch thick block of 
vesicular basalt or an 8 inch thick block of dense basalt. 

The results of these tests shoved the amount of fine material obtained by 
the corer was dependent on drilling rate, with optimum results giving about 
75Z recovery if a rate of 1 cm/ sec was maintained. 

In the end, the core retrieval system was not incorporated into the ALSD 
system. It was deemed to be too complicated for reliable astronaut operation. 

On hindsight, this judgement was probably correct. 

The boron filament reinforced fiberglass drillstem and the solid-faced bit 
were used successfully as part of the ALSD system. This type of stem was used 
to emplace the heat-flow probes. 
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TK« experience ac Apollo 15 showed that the larger diameter and loss 
of continuity of fluting at the taper joints between the sections did not 
allow cuttings to pass up the hole and jammed the drill's operation soon after 
the first joint had entered the lunar surface. For subsequent flights the 
taper joint was replaced by a titanium joint section chat was bonded by 
Epoxy into ends of the stems. These redesigned stems worked very well at 
Apollo 16 and 17 (see also Section IV). 
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APPENDIX C. DESCRIPTION OF THE LUNAR HEAT FLOW EXPERIMENT INSTRUMENTATION 

The material that appears in this section is only slightly modified 
from a paper prepared by Bruce D. Smith of the Bendix Corporation. 

EXPERIMENT THEORY 

The heat-flow Instrument performs measurements to determine the mean 
vertical temperature gradient dT/dz and the effective thermal conductivity 
k of the material across which the measured gradient is developed. Conducted 
heat flow F diffuses down a temperature gradient in accord with the relation- 
ship between these two quantities in one dimension, 

F—k (dT/dz) (1) 

Table 1 lists ranges of predicted density of heat-flow rates and 9oll 
thermal conductivities (Langseth et al., 1968; Langseth et al. , 1970; Langseth, 
1968) with the corresponding limits of average temperature gradients. 

To meaningfully extract the average gradient from the composite temperatures 
that may exist at practical lunar measuring depths, the temperature-measuring 
instrument must have a wider range than the gradients in Table 1 indicate, read- 
ings must be recorded frequently for more than a year, and absolute measurement 
accuracy must equal the relative accuracy requirement for the lowest mean temp- 
erature gradient expected. 

Two different approaches are used in measuring lunar subsurface thermal 
conductivity by means of the heat-flow instrument. In the first approach, the 
thermal response of in situ lunar material to known heat sources is tested. 

In the second, vertical strings of temperature sensors record the character- 
istics of the periodic propagations into the surface to determine diffusivlty; 
with this diffusivlty information and good estimates of soil mass density and 
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TABLE 1* Predicted Experiment Ranges 


Temperature 
Gradient dT/dz, 
* K/ meter 

Hear Flow J, 

watts/cm"2 

Thermal Conductivity K 
watts/cm ®K 

0.24 

31 

1.0 x 10“* 
6.3 x 10‘° 

4.2 x 10"* 
2.1 x 1<T° 
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SYMBOLS 


A 1...6 platinum bridge constants determined by calibration 
Ei... 6 platinum bridge constants determined by calibration 
C constant 

dTfdz mean vertical temperature gradient 

G system gain 
Is bridge excitation current 

J heat flow 
k thermal conductivity 
m mean value 
n point of constant input ratio 

digital numbers 
N\ ># *4 JV- 4096 
Rq resistance at 0*C 
R t ... 4 Wheatstone bridge arm resistances 
R l0 o resistance at 100*0 
Rg total bridge resistance 
Rf idealized output ratio 
Rt resistance at T * C 
T temperature in degrees Centigrade 
7> reference-bridge temperature 

TC thermocouple 
u digit resolution (2.4412 millivolts) 

v\ V* positive voltage 
y\ negative voltage 


** 

sensed reference-bridge excitation voltage 

»0 

reference-bridge output voltage 

"ocr 

error offset voltage referred to amplifier input 

v s 

bridge excitation voltage 

vi 

sensed bridge excitation voltage 

Vi 

potential difference across 20-ohm resistance for cur- 
rent Is. 

Vo 

bridge output voltage 

v r 

potential difference across 2-ohm resistance for cur- 
rent Is 

v x 

voltage at connection point X 

v Y 

voltage at connection point Y 

*!...♦ 

calibration constants for reference bridge 

z 

vertical distance 

Z 

attenuation ratio Ysl^S 

<*» 

solar absorpunce 

IS 

Van Ousen calibration constant 

s 

Cailendar calibration constant 

nr 

temperature difference in degrees Centigrade 

e 

system error offset 

«o 

output -sense-line offset 

«i 

excitation-sense-line offset 

•tR 

infrared emittance 

a 

standard deviation 
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specific heat, thermal conductivity can be calculated. 
EXPERIMENT DESCRIPTION 


The optimum configuration in the lunar surface for the Lunar Heat Flow 
Experiment is depicted in Figure 1. Two hollow fiberglass borestems, 2.5 
centimeters in diameter, are drilled into the lunar surface to a depth of 
3 meters at a distance of 10 meters from one another. Four sets of temperature 
sensors, spaced along a probe consisting of two flexibly joined rigid sections, 
occupy the bottom meter of each hole. The flexible joint permits the probe to 
be folded for transportation to the moon. The sensors, which are primarily 
radiatively coupled to the borestem and lunar soil, are connected electrically 
by 8-meter woven cables to a package of electronics on the surface . Each cable 
carries fo*ir precisely located thermocouple junctions in the borestems above 
the probes. The electronics unit is connected by a flat ribbon cable, 9 meters 
long, to the Apollo Lunar Surface Experiments Package (ALSEP) Central Station, 
which formats all the data from and controls the operation of all the scientific- 
station instruments on command from earth. 

The heat-flow instrument returns data giving average-temperature, differen- 
tial-temperature, and low- and high-thermal-conductivity information from four 
locations on each probe, with the thermocouples supplying readings for tempera- 
ture determinations in the upper part of the boreholes. Instrument performance 
requirements for these measurements are summarized in Table II. In the normal 
operating mode, the heat-flow instrument gathers ambient and high and low-sensitiv- 
ity differential temperature data from the "gradient" sensors situated at the ends 
of each half-probe section and samples the thermocouple outputs during the 7.25 
minute measurement sequence. Various subsequences can be selected (e.g. measure- 
ments on one probe only) , but most of them would not normally be used. Low- 
conductivity experiments are performed on command, with each heater activated 
in turn to 0.002 watt for about 40 hours. The normal measurement sequence is 
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unchanged. The high conductivity node of operation requires the selection of 
measurements on the remote sensors in any half-probe section, the type of data 
returned alternating between high-sensitivity differential and absolute tem- 
perature measurements. Either of the adjacent heaters at the ends of the 
probe half may be activated by command. Each heater should be on for about 
six hours, but this depends on the conductivity experienced. 

PROBE AM> SENSOR CONSTRUCTION 

The platinum resistance sensors used in the heat-flow Instrument probe 
contribute significantly to the quality of the measurements obtained. The 
approach to sensor construction is not new: Siemens (1871) described a pla- 

tinum resistance thermometer, with wire wound on a clay cylinder enclosed in 
a wrought-iron tube, at a meeting of the Royal Society of London in 1371. 
Whereas present commercial thermometers have stabilities of 10 to 100 milli- 
degrees Kelvin per year with normal use (Klaven et al. , 1970) the lunar gra- 
dient measurements require stabilities of the order of 1 mlllidegree Kelvin. 
These have been achieved in standards laboratories for some time, but with 
sensors susceptible to quite low levels of shock and vibration, and with sub- 
stantial precision equipment required tc obtain accurate readings. 

Two types of platinum resistance thermometers are used in the heat-flow 
instrument, the so-called 'gradient' sensor and the 'ring' or 'remote' sensor. 
The gradient sensor, the construction of which is detailed in Figure 2, in- 
corporates a unique method of supporting the resistance wire to reduce ins- 
tability normally induced by mechanical or thermal stress. Pure annealed 0.04 
millimeter diameter platinum wire is mounted, coiled in a 0.3 millimeter dia- 
meter helix extended to 0.07 millimeter pitch, on a glass-insulated platinum 
mandrel. The base of each loop is arranged so that only 10 percent of the turn 
is embedded in the substrate. The mandrel and glass have the same expansion 
coefficients as the coll, and the assembly is annealed at 673°K for 15 hours 
before sealing. The platinum coll is isolated from contamination by an 









atmosphere of pure helium contained within e gold-seeled platinum outer case. 
Platinum-wire coaxial leads extend through ceramic-insulated tubes for silver 
brazing to Svanohm and Manganin connecting wires. Each sensor assembly houses 
two separate elements, effectively bifilar wound, with nominal resistances of 
500 ohms at 273.16°K (0°C). 

The remote sensor, the construction of which is illustrated in Figure 3, 
consists of two 500- ohm-nominal- re sis tance platinum wires set in a ceramic 
glaze around a thin platinum ring. Because its intended use Is in the short-term 
high-conductivity experiments , it has a less stringent stability requirement 

o 

(0.002 K/6 hours) than the gradient sensor. Unexpectedly, however, it has 
demonstrated a long-term stability comparable to that of the gradient sensor. 

The sensors are mounted as shown in Figure 4, which illustrates the con- 
figuration common to the ends of all half-probe sections. The gradient assembly 
is epoxied at the mounting bushing at the sensor-lead exit end to the Inside of 
the probe end-sheath. The small tube on the opposite end (also shown in Figure 
2) is supported by a snugly fitting fiberglass bushing, which permits strain- 
free differential expansion with good mounting support. Associated with each 
gradient sensor is a 1000-ohm Karma wire heater, wound concentrically with the 
gradient sensor on the thin section of the end sheath. The ring-sensor platinum 
band is partially cemented internally to a filler sheath, which, in turn, is 
attached to the outer sheath of the probe. A tube joins the two end pieces to 
form a probe half-section. 

The structural components providing the span between the sensors are manu- 
factured from low-conductivity thin-walled filament-wound epoxy fiberglass. 
all cabling is carried through the probe body inside the split inner sheath, 
which is shielded with multilayer insulation to reduce radiative coupling between 
the wiring and probe wall, particularly during conductivity experiments with a 
heater on. The areas around the sensors are partially enclosed with guards for 
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pro taction during handling. 

A proba asaembly la 1.09 meters long whan unfoldad at tha closely collad 
axtanaion spring chat joins tha two half-sections. The coop late unit la coated 
with a matt -black thermal control paint. An assembled probe, with Its 8-mater 
35 conductor connecting cable, weighs less than 0.5 kilogram. 

Gradient and ring sensors are each Interconnected within a probe half with 
AWG 23 Evanobm wire to form bridges, the opposite arms of which are physically 
situated in the same sensor assembly at a common temperature. A schematic of 
this resistance-bridge arrangement is shown in Figure 5 (left), where T ^ la the 
temperature of one sensor assembly and the temperature of the other. The 
gradient-sensor assemblies, which form one bridge, are separated by 47 centi- 
meters; the remote-sensor assemblies, which form the other, are located 29 
centimeters apart and 9 centimeters from the heater windings. 

Six wires connect each bridge to the electronics unit. Those from a lower 
half-probe are conveyed in the upper section through the hollow sensor assemblies. 
Evanohm wire is used because it has a remarkably small temperature coefficient of 
resistance (0.00002/°C) , with closely controllable resistivity and a thermal 
conductivity that is low for an electrical conductor. 

SENSOR CHARACTER! S TI CS 

In Callender's empirical parabolic equation (Cal lender, 1887 and 1871) 

T - 100 [ (R t -R 0 ) / (* 1Q0 -* Q ) ] * 3 1 (T/100)-l J vs/100) (2) 

the constant 6 defines the characteristics of an individual platlx urn resistance 

0 o 

thermometer over the temperature range 0 C to 630 C, where T La temperature 
in degrees Centigrade, R 7 is resistance at temperature T, R g is resistance 

O (J 

at 0 C, and resistance at 100 C. The Callender constant is determined 

o 

by calibration at three fixed points: the triple point of water (0.0100 C) , 

O 0 

th<r* steam point (100 C) 9 and the boiling point of sulfur (444.600 C). For use 
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between 0°C and -183°C, a correction tarm devised by Van Duaan la added to 
the Callander aquation to give (Robertson and Welch, 1962) 

T - 100 [ (* r R 0 )/Ol 10Q -R Q ) ] * 5 [(T/100)-l] (T/100) 

+ 8 ( (T/100)-!] (T/100, ’ (3) 

where the constant 8 for an individual thermometer is found by calibrating 

o 

at the boiling point of oxygen (-182.97 C). The resistance of one platinum 

o 

element at temperature T relative to its resistance at 0 C la therefore: 

" 1 + (i/ioo) t <* IOjI /*w ] 

(r- 6 [ (T/100) -l ] (T/100) -8 t (T/100) -1 ] (T/100) 3 } (4) 

Referring to Figure 5, the bridge voltage ratio is defined 

V V (R i YW ' 1 ( YV Y***' 1 (5) 

where R^, Ry and R^ are Wheatstone bridge arm resistances, is the 
excitation voltage across the bridge, and v Q is the output voltage. By 
combining equations (4) and (5) and eliminating small terms, the following 
simplified expression is obtained, 

\ + A 3 T 2 + A 4 dT A j AT (AT + 2T) + 

( 6 ) 

A 6 (AT(4T 3 -300T 2 ) +AT 2 (6T 2 - 300T) + AT 3 (4T - 100) + AT 4 ] 
where R g is total bridge resistance, T is the temperature of one sensor 

assembly expressed in degrees Centigrade, AT is the temperature differential 

between sensor assemblies in degrees Centigrade, and ky ky ky ky k s , and 

k s are constants. Similarly, bridge resistance R^ can be related to temperature 

T and temperature differential AT by the simplified expression 

Rj - + B 2 T + BjT 2 + B 4 AT + BjATCAT + 2T) + B fi (T - 1U0)T 3 
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where B., B., B. t B , 3_, and B, are constants. Equations (6) and (7) are 
1 4 J 4 o o 

transcendental in T and AT. To find absolute and differential temperatures, 

an iterative simultaneous solution of both calibration equations is required 

The reduction of voltage differences and bridge resistances to temperatures 

and temperature differences was accomplished differently at the Johnson Space 

Craft Center and at L-DGO. The calibration data consisting of bridge resistance 

and voltage ratio versus bridge average temperature and temperature difference 

were fitted to two dimensional functions of the following form by least squares. 

3 ' 3 1 k 

dT - E A (i,k) .81 • C y /V ) 

4-0 * © • 

k*o 

A(4,k) are 16 coefficients derived by least squares. 

3,3 

T - 2 B (2,y) • r; . (V /V r 

i-o ® 0 « 

k«o 

where B(4,k) represent 16 coefficients. 

THE ELECTRONIC MEASUREMENT SYSTEM 
System Operation 

It can be seen from equations (6) and (7) that the electrical measurements 
required to solve for T and AT are total bridge resistance R a (measured as a 
single element), excxtation voltage V s , and differential output voltage V 
The system by which these measurements are made is represented in the simplified 
block diagram in Figure 5, which incorporates the essential features of the 
bridge measurement method, although it Includes only one of the eight probe 
bridges and omits thermocouple circuits, power supplies, and the bulk of the 
logic and control circuits. 

Each of the bridges is selected through reed relays for excitation by 

direct current. The bridges produce a differential output-to-laput voltage 

-3 o 

racio v q/ v s of approximately 15.8x10 for a dynamic range of t2 K. With 8 


k * 1 (8) 


k * l (9) 
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voles applied to Che bridge excitation cables, the data-chain maximum input 
requirement is set ee ±34 millivolts. The gradient-bridge low-sensitivity 
range of ±20°K requires 0.8 vnle excitation for a similar output maximum. 
Measurements for total bridge resistance are mede at the 8-vole supply level. 
Excitation voltage 7 g is reduced to the maximtm level of output voltage by 
an attenuator in the excitaelon sense circuit. The series elements of the 
attenuator include cable resistance that does fluctuate slightly with the 
large surface-temperature variations, but the proportion of nble resistance 
in the tamper ature-matched attenuator is small, and very high attenuation- 
ratio stability is achieved. The output impedance Is a., Tanged to be the same 
as the bridge differential output ispedance. At the low-sensitivity excitation 
level of 0.8 volt, the attenuator output vL is one- tenth the bridge 

o 

unbalance voltage 7^ at ±20 K temperature extremes. To effectively norma * 
these readings and at the same time avoid attenuator switching or gain changes, 
excitation current Xg is sensed. This current is combined with the total bridge 
resistance value, determined from readings at 8 volts excitation mede within a 
short time of the current measurement, to calculate the low-level excitation 
supply voltage. Low sensitivity current is found from the potential difference 
V developed across a precision 18f2 ohm resistance in the supply line; high- 
sensitivity current is found from the potential difference v? across the 2-ohm 
resistor only. The latter value is used with the attenuated bridge-excitation- 
voltage measurement to determine total bridge reals ance R g , Differential 
output voltages are read directly. 

These various measuremsnts are presented sequentially through a low-level 
fleld-effect-transistor (FBI) multiplexer to the input of the conmon signal 
amplifier, which has a gain of 288 and full-scale output of ±10 volts. The 
amplified signal is converted to a 13-bit digital number and clocked into a 
shift register along with a 7-bit mode-identification and binary-measurement 



TABLE 3 . Sources of Error in Experiment Data 


Platinum-bridge error sources 

* initial calibration accuracy 

* element stability with age 

* dissimilar metal EMF’s at sens^.-wire connections 


Measurement-system error sources 

* initial calibration of electronics 

* system noise 

* spurious EMF's in multiplexer 

* excitation voltage 

stability 

sense-attenuation stability 
bridge beating 

* amplifier 

common mode rejection 
gain and offset stability 
linearity 
settliug time 

* analog-to-digital converter 

linearity 

reference-voltage stability 
quantizing magnitude 


Data-processing error sources 

* truncation errors 

* iterative solution accuracy of bridge equations 
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coda. The resulting 20-bit number Is serially shifted into the ALSEF Central 
Station for insertion as two 10-bit words into the ALSEP data stream for 
transmission to earth. On earth, the binary numbers are converted back 
to the sensor temperatures from which they originated by applying the cali- 
bration factors for each bridge of thermocouple and each measuring channel. 
The sources of error to which the reconverted temperature values may be 
subject are outlined in Table III. 


THE RATIO MEASUREMENT TECHNIQUE 


Errors due to long-term system gain and offset instability are circum- 
vented by using a ratio technique, which eliminates system gain and offset 
as factors in the reconstitution of the original bridge ratio, depending for 
its success only on gain and offset stability during the, brief measurement 
period. All measurements are made twice, at two different levels. To utilize 
the full system range and obtain maximum resolution, it 'is convenient for each 
bridge Input and output measurement to be made at reversed polarity, with equal 

bipolar excitation levels. The bridge is pulsed twice in 2.4 seconds for 

-3 


2.6 milliseconds at a maximum duty cycle (0.2x10 ), which limits self -heating 

to an acceptable 0.1 microwatt; two measurements are made during each excitation 

pulse. The sequence outlined in Figure 6 for a high-sensitivity differential 

measurement is typical of all measurements. Power is applied to the bridge 

from the excitation-pulse supply at bipolar 4-volt levels, giving a positive 

8-volt total excitation. After 1 millisecond to allow the system to settle, 

the attenuator output V^+corre spending to excitation level V^+ is converted 

to a digital number Output voltage V + is then selected for measurement, 

and, after a 2.3-millisecond positive pulse duration, conversion to is 

executed. The entire process is repeated 2.4 seconds from the start of the 

sequence, with the pulse-excitation supply output reversed to -8 volts; during 

application of this negative pulse, V- and V - are converted to digital 

E 0 
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Appendix C-Figure 6 


2.6 msec*-! . h* 2.6 msec 
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numbers and N 4 , respectively . 

If G represents amplifier and analog- to-digital-converter gain, 

QJ 2 «N 4 )/(N x -N 3 ) * [(y 0 + + e 0 ♦ e)G - (y Q ~ + e Q +e)Gj/[(V' + + e x + e)G - 

<V + h • <V - V o' 1/CV Z + - 7 e‘ > 

- Q>* / TJ + )[1-C» 0 '/V 0 + )l/tl - <v-’ / Vj + )) 


where e is the amplifier and analog- co-digital-converter offset, e Q is the 
output-sense-line offset, and is the excitation-sense-line offset. Since 
the temperature of a bridge In the lunar regolith will not change during the 
2.4-second measurement period. 



/ - 7- /V' - . C 


(9) 


where C is a constant regardless of the excitation magnitudes of 
Rearranging equation (9) to read 


and V . 
E 



and combining equations (8) and (10) results in 


( 10 ) 


(N 2 - n 4 )/(n x - n 3 ) - V 0 + /v^ + (11) 

+ + 

where V /V' is the ratio between maasured bridge signals. Independent of 

0 e 

gain and offset. The required bridge ratio 7 /V is obtained by applying 

0 £ 

the excitation-sense measurement attenuation factor Z, the accuracy and 

stability of which clearly affect the result . Since Z ■ V'/V_, 

** £ 
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V 0 7 V E * V 0 /ZV E 

OTHER DESIGN FEATURES 

The use of excitation -polarity reversal, with one data converter for 

all measurements , eliminates some major sources of error and obviates 

need for separate periodic calibrations on the moon and additional circuitry 

with which to perform them. Nevertheless, numerous sources of error remain. 

It is therefore pertinent to suamarlze here the more significant features of 

the principal components of the data chain, which account for demonstrated 

measurement accuracies of better than 0.02-percent-full-scale probable error 

o o 

over a temperature range of 0 C to 60 C. 

The pulsed-pover supply produces 4.0Q0-volt±l-oillivolt bipolar pulses, 
which are stable to within 0.002 percent during bridge excitation and output 
measurements. The supply operates from positive and negative reference levels, 
derived from constant-current-driven low-temperature-coefficient zener diodes, 
which are switched by the control logic for positive and negative bridge 
excitation outputs from two series-connected operational amplifiers, each 
having a push/pull output stage to supply the 20-milliampere bridge current. 

The multiplexer for the heat-flow instrument is a double-tiered N-channel 
field-effect-transistor (FET) commutator, with 32 differential input pairs. 

It Is divided into four sections, one for each probe set of gradient sensors 
and one each for the remote sensors and the thermocouples. Each FET section 
is powered only when necessary to restrict the effects of a single-channel 
failure. Field-effect-transistor temperature-dependent offset voltages and 
mismatched differential impedances, which act with circuit impedances and 
amplifier bias currents to introduce variable offset voltages for each channel 
and dif ferret offset's between channels, are cancelled by the ratio measurement 
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technique previously described. 

the moat basic design trade-offs In aaalog-to-digital converter design 
are speed and accuracy. Neither is of great consequence in this application, 
where extreme linearity, stability, and sensitivity are the only critical 
requirements. Digital conversion of the amplified sensor signals is achieved 
by 13-bit successive approximation to the sum of a generated +I0.000-volt±l- 
nd. Ill volt offset supply and the tlO-volt full-scale output from the data 
amplifier. The offset binary number produced is linear to within 0.0075 
percent full scale, with a resolution of 2.4414 millivolts per bit. Conver- 
sion speed Is 20 microseconds per bit. The device is a conventional successive- 
approxlaatlon analog- to-digital converter, with the two most significant bits 
of the ladder network trimmed to remove errors due to voltage drops across the 
switches. 

The signals are amplified by a differential-input slngle-ended-output 
two-stage amplifier, with differential and common-mode impedances exceeding 
50 and 20 megohms, respectively, at frequencies below 5 kilohertz. The common- 
mode rejection ratio is greater than 120 decibels over the same frequency range. 
The characteristics most important to the ratio technique are linearity, which 
Is within 0.005 percent full scale, and short-term stability, which is better 
than 0.001 percent for 5 seconds at a maximum rate of temperature change of 0.0 1°C 
per second. The common-mode rejection ratio becomes important when bridge current 
is sensed at 8-volt excitation levels. 

Constant-current sources are supplied to each base of the differential 
input pair to condensate for the bias current required by the matched semi- 
conductors, but the more usual constant-current-source common-emitter supply 
is replaced with a variable current source that maintains the combined collector 
currents constant. The differential output from the collectors connects to an 
operational amplifier, the output of which feeds back to the emitters of the 
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input stags to force the total of the collector currents to divide equally 
between the two halves of the differential pair. Stable operating points 
are thus established Independently of common-mode inputs over a wide tempera- 
ture range. The first-stage gain is 100. After filtering, the signal is 
applied to the input of the second-stage 'amplifier, which normally operates 
with a gain of 2.88. The total amplifier gain of 288 is used for all bridge 
measurements . 

Thermocouple measurements require an input range of ±10 millivolts rather 
than ±34 millivolts. During thermocouple sequences, therefore, the second- 
stage gain is changed from 2.3d to 10 by switching a field-effect transistor 
to introduce a potential divider in the feedback path of the second stage. 
THERMOCOUPLE MEASUREMENTS 

Eight cons tan tan and two chrbmel wires from the thermocouple junctions 
in the probe cables connect to Kovar leads at an isothermal block inside the 
electronics package. The Kovar leads convey the thermocouple voltages to the 
multiplexer. The isothermal block also contains a platlnum/Evanohm bridge 
thermometer with two constant-value resistance arms. A schematic representing 
one probe-cable set of thermocouple junctions and the reference-temperature 
bridge is presented in Figure 7. 

Thermocouple measurements of absolute temperature have an inherent 
accuracy an order of magnitude lower than platinum resistance measurements , 
in part because of thermocouple instability but primarily because of difficulties 
In measuring low-level voltage sources. Since the ratio technique is not 
applicable to thermocouple measurements, a calibration method is used to 
establish system gain and offset during the measurement sequence that deter- 
mines the isothermal block temperature with the bridge thermometer. 

Measurements on the reference-temperature bridge follow the same pattern 
as those on the probe bridges, except that the amplitude of excitation is 
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attenuated to give a ±lQ-millivolt full-scale output, which is the range of 
the thermocouple voltages. An amplifier nominal gain of 1000 is selected 
for reference-bridge and thermocouple measurements. To find system gain 
and error offset values, it is assumed that the attenuated positive and 
negative bridge-excitation voltages v are of known, equal, and opposite 
magnitudes since they are derived from the precisely controlled pulsed- 
power supply and are connected to the bridge by short leads within the 
electronics package. 

The sensed excitation signals v + and v are converted by the analog- 
to-digital converter to binary numbers N ^ and N 3 in accord with the relation- 
ships 

■ t(u . + + V o(f )G1/U 
H 3 * t(v e‘ + Vf )G l /u 

where is error offset voltage referred to the amplifier input, G is 

system gain, and u is the digit resolution unit of 2.4414 millivolts. Gain 
is found from 

G - - N 3 )/(v e + - v e ‘) 

when the offset voltages cancel. To find the error offset voltage the 

offset binary outputs and are converted to positive and negative numbers 
by subtracting offset count 4096. If N ^ • (N ^ - 4096) and - 4096), 

then v of£ - + N'J/7. 

The use of known balanced inputs in the reference-bridge measuring sequence 
thus establishes system gain and offset. The measured outputs are converted 
through electronics calibration data, and the bridge voltage ratio is then 
used to solve iteratively a third-order calibration equation for reference- 
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bridge temperature and hence for isothermal -block ?nd reference-thermocouple 
temperature. Measurements of thermocouple voltages made shortly thereafter 
(In the sequence N^, Nj, N^, as shown on the Figure 7 schematic) are 
processed, using calculated gain and offset values, to arrive at true thermo- 
couple voltage outputs. 

The output from the chromel/constantan reference junction in the isothermal 
block is measured relative to only one of the junctions in the cable - that 
inside the hollow of the gradient sensor at the top of the probe. The remain- 
ing three junctions in each probe cable are also measured relative to this 
top junction. The double-referencing arrangement is designed to center the 
mean of the thermocouple-output extremes close to 0 volts for maximum measure- 
ment resolution. The predicted temperature range for the thermocouple junction 
at the top of the probe is 200°K to 260°K; the full-scale thermocouple range is 
90°K to 3S0°K, and the isothermal-block temperature is controlled to between 
278°K and 328°K. 

The voltage/temperature characteristics of the thermocouples are described 
by calibration correction factors applied to standard tables of the National 
Bureau of Standards (NBS) (Adams and Davisson, 1963) . 

PACKAGE CONSTRUCTION AND THERMAL CONTROL 

The heat-flow instrument operates from a 29-volt d.c. supply and requires 
data- interlace and mode-control signals from the ALSEP Central Station. The 
unit is otherwise self-contained with respect to logic and power management 
for all the sensor measurements and for probe-heater control. 

In the deployed configuration, a ribbon cable made up of 40 flat copper 
conductors in a plastic film extends from the instrument package across 9 meters 
of the lunar surface to the ALSEP Central Station. This cable is connected to 
an astTomate connector, which the astronaut mates with the Central Station 
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outlet for power Co end eozmtuni cation with the heat-flow instrument. 

The two probe cables are each made up of 35 unshielded conductors, inter- 
woven for a uniform stress distribution so that the weight of the astronaut 
can be supported without degradation in cable performance. These cables are 
very flexible, exhibiting little residual torque when extended, and they are 
covered with a woven Teflon sleeve to provide a low coefficient of friction 
during deployment. Heat leaks from the cables on the lunar surface to the 
probes and the electronics package are small since the conductors are made 
from low- thermal-conductivity wire. 

The couplets surface package, shown in figure 8, weighs 7.0 pounds 
(3.2 kilogram ) and is 9.5x10.0x11.0 inches (24x25x28 centimeters) in size, 
including the feet. The astromate connector, ribbon cable, and spool weigh 
1.4 pounds (0.6 kilogram). 

Electronics Package 

The electronic circuitry shown in Figure 9, is made up of five multilayer 
printed-circuit boards, joined in a stack by Interlocking spacers bonded to 
each board. The stack is 2. 5x7. 6x5. 6 inches (6x19x14 centimeters) in size, 
the boards consisting of up to twelve layers, each 0.004 inch (0.1 millimeter) 
thick. The layers of circuit cracks are interconnected through plated holes 
in the composite board, and as many as 210 component subpackages are mounted 
on one board. The components are in physical contact with a heavy printed- 
circuit conductor, which forms a heat path to the interlocking rpecers. 

Column screws pass through the spacers to a metal plate, which serves both 
as a thermal control plate and as a support frv the board stack. A d.c./i.c. 
converter for the instrument power supplies is packaged on a board that is 
bonded to the thermal control plate. Heating elements are also mounted in 
good thermal contact with the plate. A short length of Manganin cable (at 
left in Figure 9) provides thermal insulation from the copper conductors in 
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the ribbon cabin. Bonded Co Cha thermal control plate at the otbar side of 
the unit la a 70-pin isothermal connector to which the two probe cables are 
soldered. When the unit is assembled, the board stack is enclosed by a metal 
cover with a compressible, electrically conducting gasket to provide a barrier 
to electromagnetic interference. 

The electronics unit weighs 3.3 pounds (1.5 kilograms). 

Thermal Control Design 

A thermal insulation bag, shaped as a container for the metal cover 
surrounding the board stack, is hooked by velcro pads to a low-thermal con- 
ductivity mounting ring fitted around the inside edge of the thermal plate. 

The bag is constructed of twelve layers of closely spaced high-reflectivity 
shields having very low transverse conductivity; bridal-silk netting separates 
the layers and covers the inside and outside of the bag. The infrared emir»sivlty 
of the 1/4-mil aluminized Mylar used for the shields is approximately 0.02. 

The electronics assembly la supported and protected by a thin fiberglass 
outer case, which is connected to the mounting ring by low-conductivity joints. 
When the unit is standing on the feet on this outar case, the well-insulated 
electronics compartment is situated beneath the exposed thermal control plate. 
Internally generated heat is conducted to the plate and radiated from a spectrally 
selective surface coating (S-13C) having a high infrared emittance (e^ - 0.9) 
and a low absorptance (o g - 0.2) at frequencies where solar power is most intense. 

The thermal plate is protected from direct solar radiation by a suashield 
fitted over the assembly as shown in Figure 8. The sunshield is an insulated 
box with one open side, which is placed to f-rce away from the equ- tor with its 
edge aligned in the east /west direction. The numbered marks on the sunshield 
are used as a shadowgraph with the shadow cast by the handling tool, which fits 
into the center socket. A specular reflector slope* from the top edge of the 


n 
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lunar night ia 10.5 watts. Should Che 2 9- vole operational supply be switched 
off under abnormal circumstances, a separate survival line can be activated 
Co a part of the thermal-plate heater for a power dissipation of 4 watts. 

The thermal control system is designed for Instrument deployments between 
lunar latitudes of ±45 degrees. 

CA1IBRAII0H AND TESTING 

To ensure chat Che heat-flow instrument will meet its performance require- 
ments, the calibration aid test program is extensive and thorough. Indeed, a 
substantial part of the large array of calibration end test apparatus used was 
developed specifically for this purpose. The sensors, probes, and electronics 
are subjected to worst -case session environments and calibrated as subunits 
before assembly as flight instruments. The subassembly calibrations are verified 
by tests on the complete unit, which then undergoes a further series of tests 
closely simulating all the conditions anticipated for travel to and operation 
on the moon. The stability of the Instrument following exposure to such critical 
environments as system v 'ration is carefully checked. 

Sensor Calibration 

The platinum resistance thermometers are calibrated as a bridge by a com- 
arisorr method. Intercomparlsons of NBS standards indicate that accuracies 
afforded by this method approach those obtained by absolute fixed-point 
calibration. Each sensor is immersed, along with a standard thermometer, in 
an isothermal bath of trichloroethylene. The baths have separate temperature 
controls so that temperatures can be independently set for each thermometer 
of the differential bridge. The standard thermometers are interchanged for 
several measuremet s to determine offset. The degree Kelvin is established 
from NBS-ca librated standard thermometers. Secondary absolute-value resistance 
standards of 1000 ohms, referred to an MBS standard, ?re maintained as a ref- 
erence *oi .he electrical measurements . 
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To check for random errors, each gradient bridge is calibrated at 42 

points, which are least-squares fitted to equations (6) and (7) or 8 and 9 

to yield the constants in these equations. The standard deviation of the 

least-squares fit to the data points in equation (6) does not exceed 
-3° 

0.48x10 K for any of the sensor assemblies tested, and more than 60 

such assemblies have beev produced. The gradient bridges are calibrated 

at differential temperatures of +20, +10, +2, 0, -1, -2, -10, and -20°K 

at absolute temperatures of 200, 212.5, 225, 237.5, and 250°K. The equivalent 

dif ferential-temperature drift of some randomly selected gradient sensors 

0 o 

tested periodically at 200 K and 270 K over a 3-year period is approximately 
-3o 

0.3x10 X. The average absolute-temperature drift for individual sensor 

-3o 

elements over the same 3-year period has been found to be 0.5x10 K. The 
ring sensors are calibrated at a minimum of 14 points to calculate 12 cali- 
bration constants. 

Thermal -plate reference-temperature bridges in the electronics package 

are calibrated at -20, 0, +25, +50, and +90°C to yield the constants X , X , 

1 2 

X^ and X^ in the relationship 

V V E ■ *1 + w + w 2 + vv 3 


where T is the reference-bridge temperature and V /V is the bridge voltage 
R w £ 

o 

ratio. Duri .g reference-bridge calibration at 0 C, all eight cable thermo- 
couple junctions are immersed in an isothermal bath for calibration at 90 , 


200, 250, and 350 X. 


Electronic Testing and Calibration 

Testing of the electronics data chain involves adjustment by resistance- 
value selection for zero off-set, common-mode rejection, and gain. A calibration 
factor is found for each measuring channel by calibrating the channel as it 
will be used - that is, by performing ratio measurements. A set of resistance 
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networks, calibrated to an accuracy of 0.002 percent, substitute for bridge 
and line resistances to simulate high-sensitivity and low-sensitivity ratios 
and bridge resistance. The various types of measurement are made through 
separate channels. Calibration is performed at nine different differential- 
ratio amplitudes and five levels of bridge resistance, and at least ten 
measurements are made for each of the input ratios. The mean m and standard 
deviation ? of the output ratios for a constant input are calculated, and a 
linear calibration factor is f'und from the mean values corresponding to ratios 
near the limits of channel range (n ■ 1 and n » 9 in Figure 10) . The ideal 
transfer function of the channel is the calibration factor. The error at any 
polnt.n is defined as the difference between the ideal ov ;put R^(n) and the 
mean output a (n), summed with the standard deviation a (a). To verify system 
linearity, these measurements are made at temperature intervals throughout 
the operating range of the heat-flow instrument. A maximum error of 0.0375 
percent full scale is specified, but the typical maximum for instruments tested 
to date is 0.02 percent. Calibration factors are modified for the actual 
bridge-lead resistances. Thermocouple channel accuracy is alsc checked by 
multiple measurements at nine points throughout the +10 millivolt input range. 
Probe Performance Testing 

Assembled probes are characterized over ambient- and differential-temperature 
operating ranges in a temperature-gradient apparatus or gradient chamber, which 
has an overall height of 11 feet (3.4 raters) and a diameter of 2 feet (0.6 
meter). The probes are Inserted in a gradient tube 2.5 centimeters in internal 
diameter, in which a positive or negative temperature field can be developed. 

A liquid-nitrogen bath surrounds the gradient-tube assembly to provide a con- 
stant-temperature heat sink for heat introduced in developing average and linear 
differential temperacures . A double-vacuum shell minimizes heat leaks and 
permits the probe to be inserted without unduly disturbing the thermal equili- 
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standard deviation 
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briua. Thermocouples and thermopiles In Che gradient cube, along with 

associated readout Instrumentation, provide temperature measurements to 

accuracies that approximate those of the heat- flow Instrument. Probe- 

measured differential temperature is compared with apparatus-measured 

differential temperature in the fora of a "shorting" ratio, which is 

approximately a linear function of the absolute temperature along the 

gradient, the probes being radiatively coupled to the gradient cube. 

This shorting ratio varies from bridge to bridge, depending on small 

differences in probe construction. The probes are tested at mean ambient 

o 

temperatures of 20S, 225, and 245 K, with linear gradients of 0, 2, and 
18°K across each probe half-section. At least five tests are performed 
with each flight-instrument probe, including two in which the probe is 
integrated with the electronics. Gradient-chamber tests at 0°K check the 
zero' offset of the gradient bridges and determine the offset values for the 
ring bridges, the cables of which are shortened during probe assembly. Each 
of the probe heaters is turned on at the 0.002-watt level by the electronics 
to characterize the "low-conductivity" performance of the probe in an in- 
finitely conducting environment, which the gradient-tube apparatus represents 
since it Is an almost perfect heat sink. 

Thermal conductivity testing is also performed In an apparatus that 
simulates the properties of the lunar regolith. The thermal conductivity 

test apparatus is a 10.5-by 3-foot (3.2-by 0.9-meter) Dewar flask, with 

o o 

inner-vessel-wall cooling to the temperature range 200 K to 250 K. A probe 
bore tube extends Into a bed of glass microbeads in the inner container. The 
outer cavity and the boretube are evacuated, and the thermal conductivity of 
the glass-bead bad is varied over the anticipated lunar range by controlling 
the gas-filling pressure. Two such thermal conductivity cast units have 
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been constructed to permit the simultaneous testing of the two probes of a 
heat-flow Instrument. 

The Heat-Flow Instrument Testing 

For functional testing of the probes and the heat-flaw instrument as 
a whole, thermal slmllators are used to provide a stable temperature environ- 
ment for the probes in the lunar-temperature operating range. Each simulator 
Is an insulated container 8 feet (2.S meters) deep, filled with solid carbon 
dioxide and ethanol. In which a heavy-metal double boretube Is vertically 
situated. The ethanol fluid stabilizes at the temperature of the subliming 
dry ice that it surrounds. Since the temperature of sublimation is a function 
of pressure, the fluid temperature increases with depth to impose a gradient 
along the boretubes. Thermocouples and thermopiles within the metal block 
measure the temperature profile along the tubes. The laboratory thermal- 
simulator apparatus has a compartment for the electronics package so that 
the whole instrument can be operated in a vacuum. When inserted into the bore- 
tubes, the two probes of the instrument are in a similar, uncontrolled but 
slowly changing temperature field, and comparisons can be made of processed 
measurements on the probe bridges at given heights in the bath. 

A similar simulator has been built into an adjunct to a 27- by 20-foot 

(8.3- by 6.2 meter) thermal vacuum chamber. During ALSEF system testing, the 

heat-flow electronics package is deployed inside the main chamber, with the 

probes situated In the stabilized thermal-simulator boretubes. The entire 

ALSEF system. Including the heat-flow instrument, is exposed to the temperature 

-8 

conditions of lunar dawn, lunar noon, and lunar night in a vacuum of 5 x 10 
Torr. Cold space is simulated in the chamber by pumping liquid nitrogen under 
pressure through black, optically tight panels. The sun is represents 1 by 
carbon arc '.amps and dust degradation on instrument surfaces is simulated by 
subjecting them to additional infrared radiation. 


TABLE 4 . Apollo 16 Heat-Flow Instrument Test Results In Lunar-Night Environment 
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Table IV presents a comparison of thermal-* simulator temperature measure- 
ments made in a changing temperature environment by the heat-flow-instrument 
probes and electronics package during simulated lunar-night testing of the 
Apollo 16 ALSCF system. The means of the differences between the temperatures 
measured by adjacent bridges, corrected for bridge shorting ratio, and the 
standard deviations from the means are tabulated. The very small differences 
noted may have been caused by differences in the transient responses of the 
two probes, by different probe responses to nonlinear gradients, or by slightly 
different temperature fields, the result either of inherent differences between 
the apparatus bore tubes or of different times of measurement. Outside of these 
factors, the probe-measurement differences reflect the consistency of total 
Instrument calibration, stability, and signal processing. 

The low-sensitivity measurements on the upper gradient bridges illustrate 

an effect of data averaging that is generally seen in heat-flow-instrument test 

o 

results. Although resolution in the +20 K range is ten times less than in the 
high-sensitivity range, a fact reflected by the standard deviations, the averages 
are about the same in the two ranges. The number of measurements that do not 
follow this pattern decreases as the number of data cycles increases. Thus, 
system resolution is effectively increased by multiple-measurement averaging 
in both the low and the high sensitivity ranges. Use of this technique can 
result in an effective resolution better than analog-to-digital-converslon 
quantization would allow for a single measurement. 

During the simulated lunar-night test in which these results were obtained, 
the probe cables were lying in an uncontrolled temperature zone. There is little 
value, therefore, in comparing thermocouple readings 2,„ 3, and 4. However, the 
mean of the measurement differences for the thermocouples situated at the top 
of the probes (0.073°K) relative to the single-measurement resolution of 0.17°K 


i 



again demonstrates the data-averaging effect. The readings also provide a 
good comparative check of the thermocouples. 

INSTRUMENT PACKAGING 

All the components of the heat-flow instrument are pictured in Figure II. 
These include a packing container for the probes and an emplacement tool used 
by the astronaut to insert the probes into the borestems. Overall package 
design is strongly influenced by the requirement for deployability by a 
suited astronaut, as evidenced by the various pull rings, handles, and carry 
straps seen in the photograph. 

The probes are folded, with two molded packing pieces secured by nylon 
cloth and velcro pads holding the sections slightly apart. They are stored 
for transportation in two thin aluminum containers, which are carried to separate 
deployment sites on the lunar surface, each with one probe inside. The probes 
are held within the containers in nylon bags by soft foam bulkheads, so spaced 
along the folded probe as to suppress the first free-bending mode of the probe. 

The cables are coiled around the inside of the probe containers in troughs 
formed by an inner wall on each side. 

One probe container has provision for stowing the collapsed emplacement 
tool - four telescoping fiberglass tubes that lock when extended to give a 
tool length of 224 centimeters. A slotted piece of fiberglass at the end of 
the smallest-diamater tube fits over the probe cable to implement probe deploy- 
ment. The tool is white, with markings at 2-centimeter intervals and a bright 
orange band to indicate the depth to which the probe should be pushed. 

The two probe containers fit together to form a single package, strengthened 
by the emplacement cool inside. The package is covered with white thermal -control 
paint and secured by velcro straps with pull rings. A complete probe-container 
assembly is 3.4x4.5x25.5 inches (3.6x11.4x64.3 centimeters) in size, excluding 
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handles, and weighs 3.5 pounds (l.L kilograms). 

Carried to the moon In an equipment bay of the Lunar Module , the heat-flow 
instrument Is mounted via quick-release fasteners to a subpallet, which stands 
vertically on a pallet as shown in Figure 12, supported by energy-absorbing 
mounts. The electronics package has shock-absorbing washers fitted to the 
four lugs protruding from the mounting ring. The probe package is secured 
through tubes at the four corners of the box. The astromate connector and cable 
(coiled inside the reel) are seen in the figure beneath the electronics package. 
The nylon cover on the electronics package, which serves to protect the thermal 
control surfaces from lunar dust, is removed before final leveling and align- 
ment of the instrument on the moon. 
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